FR 3115975 A1 20220513 FR 2011409 A 20201106 Procédé non-invasif de mesure d’une grandeur physique représentative de l’élasticité de la peau, comportant : - la détermination (86) des vitesses de phase du mode fondamental des différentes composantes d'une onde de cisaillement générée à la surface de la peau, le long d'un axe de mesure, l'ensemble des paires, formées chacune par une fréquence f i et la vitesse de phase V i du mode fondamental déterminée pour cette fréquence f i , formant une courbe de dispersion du mode fondamental dans une direction de mesure parallèle à l'axe de mesure, où l'indice "i" est un numéro d'ordre de la fréquence f i et de la vitesse de phase V i , et - la transformation (100) de cette courbe de dispersion en un profil de vitesses de phase en fonction de la profondeur, la vitesse de phase à une profondeur donnée étant une grandeur physique représentative de l'élasticité de la peau à cette profondeur. Fig. 3 Procédé non-invasif de mesure d’une grandeur physique représentative de l’élasticité de la peau L'invention concerne un procédé et un appareil non-invasif de mesure d’une grandeur physique représentative de l’élasticité de la peau. La mesure des propriétés mécaniques de la peau telles que sont l’élasticité est utile, par exemple, pour aider aux diagnostics de certaines maladies de la peau ou, pour mesurer l’efficacité et les conséquences de certains traitements cosmétiques sur la peau, ou pour mesurer l’élasticité de la peau in-vitro reconstruite par ingénierie tissulaire (bio-printing), ou pour mesurer l’effet d’une stimulation mécanique de la peau in-vitro reconstruite par ingénierie tissulaire (bio-printing). Un tel procédé connu est par exemple décrit dans l’article suivant : M. Ayadh et Al : " Methods for characterizing the anisotropic behavior of the human skin’s relief and its mechanical properties in vivo linked to age effects ", IOP Publishing, Surf. Topogr.: Metrol. Prop. 8 (2020) 014002, 19/03/2020. Par la suite, cet article est désigné par le terme "Ayadh2020". Ce procédé connu est avantageux en ce qu’il est non-invasif. Il est donc simple à mettre en œuvre sur la peau d’un patient. Ce procédé connu permet de mesurer la vitesse de Rayleigh à la surface de la peau dans différentes directions. La vitesse de Rayleigh est proportionnelle à une vitesse de phase représentative de l'élasticité de la peau. La peau ressemble, par certains aspects, à un substrat multicouche. Le procédé décrit dans Ayadh2020 permet uniquement d’obtenir des informations sur les propriétés mécaniques de la couche superficielle de la peau et notamment l’épiderme. Par-contre, aucune information sur les propriétés mécaniques des couches de la peau situées sous cette couche superficielle ne peut être obtenue à l’aide de ce procédé connu. En particulier, ce procédé connu ne permet pas de mesurer les propriétés mécaniques de sous-couches telles que le derme et l’hypoderme. L’invention vise donc à proposer un procédé de mesure d’une grandeur représentative de l’élasticité de la peau qui permet en plus de mesurer une grandeur physique représentative de l’élasticité des sous-couches de la peau. Elle a donc pour objet un tel procédé non-invasif de mesure d'une grandeur physique représentative de l'élasticité de la peau, ce procédé comportant les étapes suivantes : a) la déformation de la peau au niveau d'un point d'impact à l'aide d'un stimulateur pour générer une onde de cisaillement comportant des composantes à des fréquences différentes, ces différentes composantes se propageant à la surface de la peau et provoquant un déplacement de la surface de la peau, b) la mesure, à l'aide d'un dispositif de mesure, du déplacement de la surface de la peau au cours du temps en au moins trois points de mesure alignés les uns derrière les autres le long d'un axe de mesure, dans lequel le procédé comporte également les étapes suivantes : c) la détermination des vitesses de phase du mode fondamental des différentes composantes de l'onde de cisaillement générée, le long de l'axe de mesure, à partir des mesures du dispositif de mesure, l'ensemble des paires, formées chacune par une fréquence f i et la vitesse de phase V i du mode fondamental déterminée pour cette fréquence f i , formant une courbe de dispersion du mode fondamental dans une direction de mesure parallèle à l'axe de mesure, où l'indice "i" est un numéro d'ordre de la fréquence f i et de la vitesse de phase V i , et d) la transformation de cette courbe de dispersion en un profil de vitesses de phase en fonction de la profondeur, la vitesse de phase à une profondeur donnée étant une grandeur physique représentative de l'élasticité de la peau à cette profondeur. Les modes de réalisation de ce procédé peuvent comporter une ou plusieurs des caractéristiques suivantes : 1) la transformation de la courbe de dispersion en un profil de vitesses de phase dans la direction de mesure comporte les opérations suivantes : o1) convertir chaque fréquence f i de la courbe de dispersion, en une longueur d'onde λ i correspondante à l'aide de la relation suivante : λ i = V i /f i , puis o2) calculer la valeur d'un coefficient α pour cette direction de mesure à l'aide de la relation suivante : α = Z max / λ max , où : - Z max est égal à la moitié de la distance qui sépare les deux points de mesure les plus éloignés l'un de l'autre le long de l'axe de mesure, et - λ max est la plus grande des longueurs d'onde λ i obtenue à l'issue de l'exécution de l'opération o1) pour cette direction de mesure, o3) convertir chaque longueur d'onde λ i obtenue à l'issue de l'opération o1) en une profondeur p i correspondante à l'aide de la relation suivante : p i = αλ i , où α est le coefficient calculé lors de l'opération o2). 2) le procédé comporte l'exécution des étapes a), b), c) et d) pour au moins un premier et un second axes de mesure décalés angulairement l'un par rapport l'autre, ces premier et second axes de mesure passant par le même point d'impact. 3) l'étape c) de détermination des vitesses de phase du mode fondamental comporte la réitération des opérations suivantes pour plusieurs fréquences f i différentes, : - pour chaque point de mesure : - le filtrage, à l'aide d'un filtre passe-bande centré sur la fréquence f i et dont la bande passante à -3dB est comprise entre les fréquences f i-1 et f i+1 , du signal u(x,t), mesuré au niveau du point de mesure de coordonnée x le long de l'axe de mesure par le dispositif de mesure, pour obtenir un signal filtré u i (x,t), - l'identification d'un instant t i,m (x) auquel le signal filtré u i (x,t) passe par son minimum absolu, puis - le calcul de la vitesse V i de propagation de ce minimum le long de l'axe de mesure à partir des instants t i,m (x) et des positions x auxquels ce minimum se produit, la vitesse V i ainsi calculée étant la vitesse de phase du mode fondamental à cette fréquence f i . 4) l'étape c) comporte : - une opération d'identification automatique, parmi l'ensemble des fréquences f i , d'une fréquence minimale f min en dessous de laquelle les fréquences f i ne vérifient plus la condition (1) suivante : où : - µ i,1 et µ i,0 sont les coefficients de la droite, déterminée par la méthode des moindres carrés, approchant au mieux les points de coordonnées (x; t i,m (x)), - x p est égal à la position x du p-ième point de mesure compté à partir du premier point de mesure le plus proche du point d'impact, - P max est égal au nombre de points de mesure répartis le long de l'axe de mesure, - err max est une constante prédéterminée, - ensuite seule les fréquences f i supérieures ou égales à cette fréquence f min sont retenues pour former la courbe de dispersion. 5) l'étape c) comporte : - une opération d'identification automatique, parmi l'ensemble des fréquences f i , d'une fréquence maximale f max au-dessus de laquelle les fréquences f i ne vérifient plus la condition (1), cette opération d'identification automatique de la fréquence f max étant réalisée en testant cette condition (1) pour plusieurs fréquences f i , ensuite - seule les fréquences f i inférieures ou égales à cette fréquence f max sont retenues pour former la courbe de dispersion. 6) toutes les fréquences f i sont comprises entre 1 Hz et 3000 Hz. L’invention a également pour objet un appareil non-invasif de mesure d’une grandeur physique représentative de l’élasticité de la peau, pour la mise en œuvre du procédé ci-dessus, cet appareil comportant : - un stimulateur apte à déformer la peau au niveau d'un point d'impact pour générer une onde de cisaillement comportant des composantes à des fréquences différentes, ces différentes composantes se propageant à la surface de la peau et provoquant un déplacement de la surface de la peau, - un dispositif de mesure apte à mesurer le déplacement de la surface de la peau au cours du temps en au moins trois points de mesure alignés les uns derrière les autres le long d'un axe de mesure, dans lequel l'appareil comporte une unité de traitement apte : - à déterminer les vitesses de phase du mode fondamental des différentes composantes de l'onde de cisaillement générée, le long de l'axe de mesure, à partir des mesures du dispositif de mesure, l'ensemble des paires, formées chacune par une fréquence f i et la vitesse de phase V i du mode fondamental déterminée pour cette fréquence f i , formant une courbe de dispersion du mode fondamental dans une direction de mesure parallèle à l'axe de mesure, où l'indice "i" est un numéro d'ordre de la fréquence f i et de la vitesse de phase V i , et - à transformer cette courbe de dispersion en un profil de vitesses de phase en fonction de la profondeur, la vitesse de phase à une profondeur donnée étant une grandeur physique représentative de l'élasticité de la peau à cette profondeur. Les modes de réalisation de cet appareil peuvent comporter une ou plusieurs des caractéristiques suivantes : 1) - le dispositif de mesure comporte un réseau de capteurs apte chacun à mesurer l'amplitude de la déformation de la surface de la peau en un point de mesure respectif, ce réseau comportant au moins trois capteurs qui mesurent chacun le déplacement de la surface de la peau en trois points de mesure respectifs alignés les uns derrière les autres le long d'un axe de mesure, et - l'appareil comporte un bras articulé sur lequel est monté le réseau de capteurs, ce bras articulé étant apte à faire tourner le réseau de capteurs d'un angle prédéterminé autour d'un axe de rotation, pour aligner l'axe de mesure du réseau de capteur sur un premier axe de mesure et, en alternance, sur un second axe de mesure décalé angulairement par rapport au premier axe de mesure. 2) - le réseau de capteurs comporte : - une ligne de capteurs optiques qui captent la lumière réfléchie par chaque point de mesure, et - un microprocesseur configuré pour déterminer le déplacement de la surface de la peau au niveau de chaque point de mesure éclairé à partir de la lumière réfléchie captée par la ligne de capteurs optiques, - le dispositif de mesure comporte un émetteur d'un faisceau lumineux qui éclaire chaque point de mesure aligné le long de l'axe de mesure. 3) le stimulateur est apte à projeter sur la surface de la peau un jet de fluide qui provoque la déformation de la peau au niveau du point d'impact. L’invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d’exemple non limitatif et faite en se référant aux dessins sur lesquels : - la est une illustration schématique de l’architecture d’un appareil de mesure d’une grandeur physique représentative de l’élasticité de la peau ; - la est une illustration schématique d’un stimulateur et d’un dispositif de mesure de l’appareil de la ; - la est un organigramme d’un procédé de mesure d’une grandeur physique représentative de l’élasticité de la peau à l’aide de l’appareil de la ; - la est un graphe en trois dimensions illustrant des mesures de déplacement acquises par l’appareil de : - la est un graphe illustrant le calcul de la vitesse de déplacement d’une onde de cisaillement ; - la est un graphe en trois dimensions illustrant, pour une composante fréquentielle, des mesures de déplacement acquises par l’appareil de la ; - la est un graphe illustrant la détermination d’une vitesse de phase de la composante fréquentielle de la ; - la est un graphe illustrant l’amortissement de la composante fréquentielle de la ; - la est un graphe représentant une courbe de dispersion déterminée à l’aide de l’appareil de la ; - la est un graphe illustrant un profil de vitesses de phase déterminé à l’aide de l’appareil de la ; - la est un graphe illustrant une tomographie construite par l’appareil de la ; - la est un graphe illustrant l’amortissement de la composante fréquentielle de la en fonction de la profondeur ; - la est un graphe illustrant, sur une même représentation graphique, les vitesses de phase mesurées à l’aide du procédé de la et celles déterminées à l’aide d’un autre procédé connu de détermination de vitesses de phase. Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes éléments. Dans la suite de cette description, les caractéristiques et fonctions bien connues de l’homme du métier ne sont pas décrites en détail. Dans cette description, un exemple détaillé de mode de réalisation est d’abord décrit dans un chapitre I en référence aux figures. Ensuite, dans un chapitre II, des variantes de ce mode de réalisation sont présentées. Enfin, les avantages des différents modes de réalisation sont présentés dans un chapitre III. Chapitre I : Exemple de mode de réalisation La représente un appareil 2 non-invasif de mesure d’une grande physique représentative de l'élasticité de la peau. L’appareil 2 comporte : - un bras articulé 4, - un stimulateur 8 et un dispositif 10 de mesure, fixés sur une extrémité distale du bras 4, et - une unité 12 de traitement informatique raccordée au stimulateur 8 et au dispositif 10 de mesure. Une extrémité proximale du bras 4 est fixée sans aucun degré de liberté sur un support fixe 14. Le bras 4 comporte plusieurs articulations 20 qui permettent de déplacer à la fois le stimulateur 8 et le dispositif 10 de mesure simultanément en rotation autour d’un axe 22 de projection. Sur la , l’axe 22 est vertical. Ici, les articulations 20 permettent de déplacer le stimulateur 8 et le dispositif 10 de mesure par rapport au support 14 avec six degrés de liberté. Une fois que le bras 4 a été déformé pour placer le stimulateur 8 et le dispositif 10 de mesure dans la position souhaitée, le bras 4 maintient ce stimulateur 8 et ce dispositif 10 immobiles dans cette position. Les articulations 20 sont, par exemple, actionnées manuellement par un utilisateur ou actionnées par des moteurs électriques. Le bras 4 sert aussi de support pour des conducteurs électriques qui raccordent le stimulateur 8 et le dispositif 10 de mesure à l’unité 12 de traitement. Le stimulateur 8 déforme la peau d’un être humain au niveau d’un point d’impact lorsqu'il est actionné. Le point d'impact est situé à l’intersection de l’axe 22 de projection et de la surface de la peau. La déformation de la peau produite par le stimulateur 8 est telle qu’elle génère une onde de cisaillement qui se propage ensuite le long de la surface de la peau. Cette onde de cisaillement comporte des composantes fréquentielles à plusieurs fréquences différentes. Typiquement, dans le cas de la peau humaine, les fréquences de ces composantes se situent entre 1 Hz et 3000 Hz et, habituellement, entre 1 Hz et 1000 Hz. Le dispositif 10 mesure la déformation de la surface de la peau, causée par l’onde de cisaillement, en plusieurs points de mesure alignés les uns derrière les autres le long d’un axe de mesure. Cet axe de mesure s'étend parallèlement à une direction appelée "direction de mesure" dans ce texte. Typiquement, le dispositif 10 comporte plus de trois, dix ou cent points de mesure. Ici, le dispositif 10 comporte 400 points de mesure. Par la suite, la position d’un point de mesure le long de l’axe de mesure est repérée par son abscisse x, par exemple exprimée en mm ou en µm, mesurée à partir d’une origine O. Ici, l’axe de mesure et l’axe 22 de projection se coupent sensiblement à angle droit. L’origine O est prise égale au point d’intersection de ces deux axes. La distance entre le point de mesure le plus proche du point d’impact et le point de mesure le plus éloigné du point d’impact est notée L max . Ici, à titre d’illustration, la distance L max est égale à 7 mm. Dans ce mode de réalisation, les points de mesure sont uniformément repartis le long de l’axe de mesure. La distance entre deux points de mesure consécutifs est donc égale ici à 17,5 µm. L’unité 12 de traitement est raccordée au dispositif 10 pour acquérir les mesures de ce dispositif. Plus précisément, l’unité 12 acquiert au cours du temps, avec une fréquence d’échantillonnage f e , le déplacement mesuré en chacun des points de mesure. Par la suite, le déplacement mesuré au niveau d’un point de mesure d’abscisse x à l’instant t est noté u(x, t). Par exemple, ici, la fréquence f e est égale à 8 kHz. L’unité 12 est capable de traiter les signaux u(x, t) ainsi acquis pour en extraire une grandeur physique représentative de l'élasticité de la peau à différentes profondeurs. À cet effet, l’unité 12 comporte une unité centrale 30 et une interface machine 32. L’unité centrale 30 comporte : - une mémoire 34 comportant des instructions pour exécuter le procédé de la , et - un microprocesseur 36 programmable apte à exécuter les instructions enregistrées dans la mémoire 34. L’interface 32 permet d’afficher la grandeur physique mesurée représentative de l'élasticité de la peau. Typiquement, à cet effet, l’interface 32 comporte un écran 38. Ici, l’interface 32 comporte aussi un clavier 40 pour, par exemple, acquérir une commande de déclenchement de l’exécution du procédé de mesure de la . La représente plus en détail le stimulateur 8 et le dispositif 10 de mesure. Sur la , la peau est schématiquement représentée sous la référence 46 et la surface de la peau 46 porte la référence 48. Sur cette figure, la surface 48 est représentée sous une forme déformée après avoir subi un impact appliqué par le stimulateur 8. Le point d’impact sur la surface de la peau porte la référence 49. Pour créer l’onde de cisaillement qui se propage à la surface de la peau, le stimulateur 8 utilise, dans ce mode de réalisation, un jet d’air qui vient percuter la peau au niveau du point d’impact 49. À cet effet, le stimulateur 8 comporte : - un réservoir 50 d’air sous pression, - un détendeur 52 fluidiquement raccordé au réservoir 50, - une électrovanne commandable 54, et - une buse 56 fluidiquement raccordée à une sortie de l’électrovanne 54. Par exemple, la pression d’air contenue dans le réservoir 50 est supérieure à 0,6 ou 0,8 MPa. Ici, le détendeur 52 abaisse la pression d’air. Par exemple, la pression de l’air en sortie du détendeur 52 est comprise entre 0,1 MPa et 0,6 MPa ou entre 0,1 et 0,4 MPa. L’électrovanne 56 est déplaçable entre une position ouverte et une position fermée et vice versa sous la commande de l’unité 12. Dans la position fermée, l’électrovanne 54 empêche l’air de s’échapper du réservoir 50. Dans la position ouverte, au contraire, l’électrovanne 54 autorise l’air à s’échapper du réservoir 50. L’air qui s’échappe du réservoir 50 est alors guidé par la buse 56 pour former un jet d’air le long de l’axe 22 qui vient percuter la surface 48 de la peau 46 au niveau du point d’impact 49. L’électrovanne 54 permet de régler la durée du jet d’air projeté sur la surface 48. Typiquement, la durée du jet d’air est inférieure à 20 ms ou 10 ms. Ici, la durée du jet d’air est comprise entre 5 ms et 10 ms. Au moins une partie de la buse 56 s’étend le long de l’axe 22 pour diriger le jet d’air le long de cet axe. L’extrémité de la buse 56, en vis-à-vis de la surface 48, est mécaniquement séparée de cette surface, de sorte qu’il n’existe aucun contact mécanique direct entre l’appareil 2 et la surface 48 lors de son utilisation. Le dispositif 10 de mesure est un dispositif optique de mesure. Dans ce mode de réalisation, il comporte à cet effet : - un émetteur 60 d’un faisceau lumineux qui éclaire chaque point de mesure, - une ligne 62 de capteurs optiques qui captent la lumière réfléchie par chacun des points de mesure, et - un microprocesseur 64 programmé pour déterminer le déplacement de la surface 48 au niveau de chacun des points de mesure à partir de la lumière réfléchie captée par chacun des capteurs optiques. Ici, le dispositif 10 est positionné par rapport au stimulateur 8 de sorte que le point de mesure le plus proche du point d’impact 49 soit séparé de ce point d’impact par une distance supérieure à 0,5 mm ou 0,8 mm et, typiquement, inférieure à 5 mm. De préférence, la distance entre le point de mesure le plus proche du point d'impact et ce point d'impact est comprise entre 0,7 mm et 1,3 mm ou entre 0,9 mm et 1,1 mm. Ici, cette distance est égale à 1 mm. De tels dispositifs de mesure sont connus et commercialisés. Par exemple, dans ce mode de réalisation, le dispositif 10 est celui commercialisé par la société KEYENCE ® sous la référence LJ-V 7020. Dans ce cas, l’émetteur 60 est une source laser qui émet un faisceau lumineux monochromatique et collimaté. Ainsi, le dispositif 10 n’est pas décrit plus en détail par la suite. Le fonctionnement de l’appareil 2 va maintenant être décrit à l’aide de la et en référence aux graphes des figures 4 à 13. Lors d’une étape initiale 68, le bras 4 est déformé pour placer le stimulateur 8 et le dispositif 10 de mesure à proximité d’une partie du corps humain recouverte de la peau à étudier. L’appareil 2 permet d’étudier n’importe quelle partie du corps humain. Par exemple, les résultats expérimentaux présentés sur les figures 4 à 13 ont été obtenus sur l’avant-bras d’un être humain. Typiquement, le stimulateur 8 et le dispositif 10 sont placés par rapport à la surface 48 de la peau de manière à ce que l’axe 22 de projection fasse un angle compris entre 75° et 110°, et de préférence entre 80° et 100°, par rapport à la direction normale à la peau au niveau du point d’impact 49. L’extrémité inférieure de la buse 56, tournée vers la surface 48, est séparée du point d’impact 49 par une distance supérieure à 1 mm ou 2 mm ou 5 mm, et généralement inférieure à 20 mm. Une fois l’appareil 2 correctement positionné par rapport à la peau 46, il est procédé à une phase 70 d’acquisition des signaux u(x, t). Plus précisément, une fois que le stimulateur 8 et le dispositif 10 sont correctement positionnés par rapport à la surface 48, lors d’une étape 72, l’unité 12 commande le stimulateur 8 pour provoquer l’émission d’un jet d’air et déformer la peau. Ici, l’unité 12 commande l’électrovanne 54 pour générer ce jet d’air. Ce jet d’air percute alors la peau au niveau du point d’impact 49. Cela provoque une brève déformation de la surface 48 au niveau de ce point d’impact 49. Cette déformation de la surface 48 génère à son tour une onde de cisaillement qui se propage le long de la surface 48 dans toutes les directions, et donc notamment le long de l’axe de mesure du dispositif 10. Lors d’une étape 74, le dispositif 10 mesure, en chacun des points de mesure, le déplacement de la surface 48 causé par l’onde de cisaillement qui se propage à la surface de la peau. En parallèle, lors d’une étape 76, à chaque instant d’échantillonnage et pour chaque point de mesure, l’unité 12 acquiert la mesure du dispositif 10. Ainsi, l’unité 12 acquiert chacun des signaux u(x, t). La représente sur un graphe en trois dimensions un exemple des signaux u(x, t) acquis. Dans ce graphe : - l’axe horizontal représente le temps en millisecondes, - l’axe vertical représente l’amplitude, en millimètres, du déplacement de la surface 48, et -l’axe des profondeurs représente, en µm, la position x du point de mesure le long de l’axe de mesure. Une fois les différents signaux u(x, t) acquis, la phase 70 d’acquisition se termine et une phase 80 de traitement des signaux u(x, t) débute. La phase 80 est réalisée par l’unité 12. Lors d’une étape 82, pour chaque signal u(x, t), l’unité 12 recherche et identifie l’instant t min (x) auquel le signal u(x, t) passe pas son minimum absolu. Ensuite, lors d’une étape 84, l’unité 12 détermine l’équation de la droite D qui approche au mieux le nuage de points formé par les points de coordonnées (x ; t min (x)). L’équation de la droite est la suivante : t = µ 1 x + µ 0 . Les coefficients µ 1 et µ 0 sont ceux obtenus en mettant en œuvre la méthode des moindres carrés. La représente le nuage de points formé par les points de coordonnées (x ; tmin(x)) et la droite D obtenue par la méthode des moindres carrés. L’onde de cisaillement qui se propage dans la peau 46 est une onde dispersive, c’est-à-dire que sa vitesse de phase dépend de la fréquence. La vitesse de phase n'est donc pas la même pour chaque composante fréquentielle de cette onde de cisaillement. Par vitesse de phase, on désigne ici, en absence d’indication contraire, la vitesse de phase du mode fondamental de l’onde de cisaillement. Le mode fondamental correspond au mode pour lequel l’amplitude du déplacement de la surface 48 est maximale. Cette vitesse de phase est notée V i pour la fréquence f i , où l’indice i est un numéro d’ordre identifiant la fréquence f i . Ici, les fréquences sont classées de la plus petite fréquence, notée f 1 , à la plus grande fréquence notée f Pmax . Dans le cas de la peau, il a été observé que les différentes composantes fréquentielles qui se propagent dans la peau 46 sont typiquement comprises entre 1 Hz et 3000 Hz et, le plus souvent, entre 1 Hz et 1000 Hz ou entre 1 Hz et 500 Hz ou entre 1 Hz et 400 Hz. A titre d'illustration, ici, les fréquences f i sont choisies dans cet intervalle allant de 1 Hz à 1000 Hz. La plus petite fréquence, notée f min , pour laquelle il existe une vitesse de phase a une importance particulière. Il a été observé que cette plus petite fréquence f min est généralement comprise entre 1 Hz et 10 Hz. Ainsi, dans ce procédé, le pas d’échantillonnage des fréquences f i dans l’intervalle [1 Hz ; 10 Hz] est choisi petit, c’est-à-dire ici inférieur à 2 Hz ou 1 Hz. À l’inverse, dans l’intervalle [10 Hz, 1000 Hz], le pas d’échantillonnage des fréquences f i est choisi plus grand. Par exemple, dans cet intervalle [10 Hz ; 1000 Hz], le pas d’échantillon est supérieur à 5 Hz ou 10 Hz ou 20 Hz. Ainsi, les fréquences f i sont séparées les unes des autres par un pas de 1 Hz dans l’intervalle [1 Hz ; 10 Hz] alors qu’elles sont séparées les unes des autres par un pas supérieur à 5 Hz ou 10 Hz dans l’intervalle [10 Hz ; 1000 Hz]. Lors d’une étape 86, l’unité 12 détermine, pour chaque fréquence f i choisie, si elle existe, la vitesse de phase V i correspondante. Pour cela, lors d’une opération 88, l’unité 12 filtre, à l’aide d’un filtre passe-bande centré sur la fréquence f i , chaque signale u(x, t). Le signal u(x, t) filtré à la fréquence f i est noté par la suite u i (x, t). La bande passante à - 3 dB de ce filtre passe-bande est comprise entre les fréquences f i-1 et f i+1 . Typiquement, cette bande passante est inférieure à 20 Hz ou 10 Hz pour les fréquences f i comprises dans l'intervalle [10 Hz; 1000 Hz] et inférieure à 2 Hz pour les fréquences f i comprises dans l'intervalle [1 Hz; 10 Hz]. La représente différents signaux ui(x, t) obtenus en filtrant, à la fréquence de 20 Hz, les signaux u(x, t) de la . Ce graphe est identique au graphe de la sauf qu’il représente les signaux ui(x, t) et non pas les signaux u(x, t). Lors d’une opération 90, l’unité 12 recherche et identifie, pour chaque signal u i (x, t), l’instant t i,m (x) où ce signal passe par son minimum absolu. Pour cela, si l’instant t i-1 , m (x) ou t i+1 , m (x) a déjà été identifié pour, respectivement, la fréquence f i-1 ou la fréquence f i+1 , l’unité 12 recherche l’instant t i , m (x), en priorité, dans un intervalle de temps centré sur cet instant t i-1 , m (x) ou t i+1,m (x). Si aucun instant t i-1,m (x) ou t i+1,m (x) n’a précédemment été identifié, alors l’instant t i,m (x) est recherché dans un intervalle de temps centré sur l’instant t min (x) identifié lors de l’étape 82. Lors d’une opération 92, l’unité 12 calcule la vitesse de phase V i pour la fréquence f i à partir des signaux u i (x, t). Pour cela, l’unité 12 détermine l’équation de la droite D i qui approche au mieux le nuage de points formé par les points de coordonnées (x ; t i,m (x)). L’équation de la droite D i est la suivante : t i,e (x) = µ i,1 x + µ i,0 . Les coefficients µ i,1 et µ i,0 sont obtenus en mettant en œuvre la méthode des moindres carrés comme décrit pour l’étape 84. La représente le nuage de points de coordonnées (x, ti,m(x)) ainsi que la droite Di. Une fois que l’équation de la droite D i est déterminée, l’unité 12 estime aussi l’erreur d’approximation, c’est-à-dire l’écart qui existe entre les points de coordonnées (x ; t i,m (x)) et les points de coordonnées (x ; t i,e (x)) qui ont les mêmes abscisses et qui se trouvent sur la droite D i . Si cette erreur dépasse un seuil prédéterminé err max , alors il est considéré qu’il n’existe pas de vitesse de phase pour la fréquence f i car l’énergie de la composante fréquentielle de fréquence f i de l’onde de cisaillement est négligeable. Ici, par exemple, l’erreur d’approximation est estimée à l’aide de la relation (1) suivante : où : - µ i,1 et µ i,0 sont les coefficients de la droite D i , - x p est égal à la position x du p-ième point de mesure compté à partir du premier point de mesure le plus proche du point d'impact, - P max est égal au nombre de points de mesure répartis le long de l'axe de mesure, - err max est une constante prédéterminée. Si cette erreur d’approximation est inférieure ou égale au seuil err max , alors la vitesse de phase V i est prise égale à 1/µ i,1 . La vitesse de phase V i ainsi obtenue est la vitesse de phase du mode fondamental à la fréquence f i . En effet, elle est seulement obtenue à partir des minimums des signaux u i (x, t), c’est-à-dire à partir des points où l’amplitude du déplacement est maximale. Lors d’une opération 94, l’unité 12 détermine aussi l’atténuation A i (x) de la composante fréquentielle à la fréquence f i . Ici, l’atténuation A i (x) est prise égale à l’amplitude du minimum du signal u i (x, t) identifié lors de l’opération 90. Ainsi, dans ce mode de réalisation, l’atténuation A i (x) est égale à u i (x ; t i , m (x)). La représente l’évolution de l’atténuation Ai(x) en fonction de la position x pour les signaux ui(x,t) de la . Les opérations 88 à 94 sont réitérées pour chacune des fréquences f i . Par la suite, on note f min et f max , respectivement, la plus petite et la plus grande des fréquences f i pour lesquelles une vitesse de phase V i a été déterminée. La courbe formée par l’ensemble des points de coordonnées (V i ; f i ) est appelée courbe de dispersion. La représente la courbe de dispersion obtenue à partir des signaux u(x,t) de la . Plus la longueur d’onde λ i d’une composante fréquentielle de l’onde de cisaillement est grande, plus cette composante fréquentielle se propage en profondeur sous la surface 48 de la peau 46. Ainsi, la vitesse V i d’une composante fréquentielle de longueur d’onde λ i représente les propriétés mécaniques de la peau à une profondeur p i . À ce stade, il est rappelé que la vitesse de phase V i est reliée aux propriétés mécaniques de la peau 46 par la relation suivante : V i =(E i /(2ρ i (1+ν i ))) 0,5 , où : - E i est le module de Young de la peau à la profondeur p i , - ρ i est la densité volumique de la peau à la profondeur p i , et - ν i est le coefficient de poisson de la peau à la profondeur p i . Ainsi, si l’on considère que la densité volumique ρ i et le coefficient ν i sont des constantes connues, la vitesse V i est directement représentative du module de Young E i à la profondeur p i . Ainsi, la vitesse V i est représentative de l’élasticité de la peau à la profondeur p i . Dès lors, le profil des vitesses V i en fonction de la profondeur p i correspond à une vue en coupe, le long de l’axe de mesure, des propriétés mécaniques de la peau. Lors d’une étape 100, l’unité 12 transforme la courbe de dispersion en un profil de vitesses en fonction de la profondeur p i . Pour cela, lors d’une opération 102, l’unité 12 convertit chaque fréquence f i en une longueur d’onde λ i correspondante. A cet effet, l’unité 12 utilise la relation (2) suivante : λ i = V i /f i . Ensuite, lors d’une opération 104, l’unité 12 calcule la valeur d’un coefficient α qui permet de convertir chaque longueur d’onde λ i en une profondeur p i correspondante. Ici, la profondeur p i est la distance qui sépare un point, enfoui sous la peau, de la surface 48 de la peau. Ce coefficient α relie chaque longueur d’onde λ i à la profondeur p i correspondante selon la relation (3) suivante : p i = αλ i . Ce coefficient α est une constante pour un axe de mesure donné. Par contre, dans le cas de matériaux viscoélastiques anisotropes comme la peau, le coefficient α varie en fonction de la direction de mesure. Autrement dit, le coefficient α dépend de la direction dans laquelle les mesures sont faites. Ici, le coefficient α est calculé par l’unité 12 à l’aide de la relation (4) suivante : α = Z max /λ max , où : - Z max est égal la moitié de la distance L max , et - λ max est la plus grande des longueurs d’onde λ i obtenues à l’issue de l’opération 102. Une fois la valeur du coefficient α calculée, lors d'une opération 106, l’unité 12 convertit chaque longueur d’onde λ i en une profondeur p i à l’aide de la relation (3) ci-dessus. Ainsi, une fréquence f i correspond à une longueur d'onde λ i qui elle-même correspond à une profondeur p i . En remplaçant dans chaque point de coordonnées (V i ; f i ), la fréquence f i par la profondeur correspondante p i , le profil de vitesses mesuré par l’appareil 2 est obtenu. Un exemple d’un tel profil de vitesses est représenté sur la . Ici, l’étape 68 et les phases 70 et 80 sont réitérées plusieurs fois en faisant tourner à chaque fois l’axe de mesure d’un angle prédéterminé autour de l’axe 22. Pour cela, lors de chaque nouvelle itération de l’étape 68, le bras 4 est déformé de manière à faire tourner sur eux-mêmes le stimulateur 8 et le dispositif 10 de mesure. Cette rotation ne modifie pas la position de l’axe 22 de projection, et donc la position du point d’impact 49. Par exemple, lors de chaque nouvelle itération de l’étape 68, l’axe de mesure est décalé angulairement par rapport à sa précédente position d’au moins 1° ou 5°, et par exemple de 10° ou 20°. Lors de chaque nouvelle itération de la phase 80 de traitement des données acquises, l’étape 100 est de nouveau exécutée. En effet, comme expliqué précédemment, la valeur du coefficient α dépend fortement de la direction de l’axe de mesure dans le cas de la peau humaine. Enfin, lors d’une étape 110, les propriétés mécaniques de la peau en fonction de la profondeur pi sont affichées sur l’écran 38. Différentes représentations graphiques sont possibles. Par exemple, le profil de vitesses, tel que celui représenté sur la , est affiché sur l’écran 38. Toutefois, de préférence, les différents profils de vitesses obtenus pour différentes directions de mesure sont simultanément affichés sur un même graphe pour former une tomographie de la peau 46 au niveau du point d’impact 49. Une telle tomographie est représentée sur la . Dans cette tomographie, l’axe des ordonnées représente la profondeur pi. L’axe 22 correspond à l’axe de projection de l’appareil 2. Différents plans de mesures Pl1 à Pl6 sont représentés. Ces plans Pl1 à Pl6 sont décalés angulairement les uns par rapport aux autres. Ici, chacun de ces plans de mesure contient l’axe 22 et s’étend parallèlement à une direction de mesure respective. Chaque plan de mesure contient le profil de vitesses mesuré le long de l’axe de mesure parallèle à cette direction de mesure respective. Les axes des abscisses représentent donc les coordonnées Vix et Viy de la vitesse Vi mesurée. La représente l’atténuation Ai(x) en fonction de la profondeur. Plus précisément, l’axe des abscisses représente la position x du point de mesure. L’axe des ordonnées représente la profondeur pi en millimètre. La couleur de chaque point de coordonnées (x ; p) code l’atténuation de la vitesse Vi pour cette profondeur p et cette abscisse x. Pour cela, la profondeur est transformée en une longueur d’onde λi correspondante à l’aide de la relation (3), puis la longueur d’onde λi ainsi obtenue est transformée en une fréquence fi correspondante à l’aide de la relation (2). Lors de l’étape 94, l’atténuation Ai(x) en fonction de la position x a été relevée pour toutes les fréquences fi, et donc pour la fréquence fi particulière correspondant à cette profondeur p. Ceci a été illustré sur le graphe de la . Sur le graphe de la , l’atténuation Ai(x) correspondant à l’abscisse x est relevée. C’est la valeur de cette atténuation Ai(x) relevée qui est codée par une couleur particulière au niveau du point de coordonnées (x ; p) dans le graphe de la . Les différentes vitesses de phase V i peuvent être déterminées, à partir des signaux u(x, t), par d’autres méthodes, comme par exemple la méthode connue sous l’acronyme MASW (Multichannel Analysis of Surface Waves). La méthode MASW permet de déterminer, pour chaque fréquence f i , la vitesse V i du mode fondamental ainsi que les vitesses de phase des modes d’ordre supérieur à l’ordre du mode fondamental. La représente sur un même graphique la vitesse de phase Vi déterminée par le procédé de la et les différentes vitesses de phase déterminées par la méthode MASW. L’axe des abscisses représente la fréquence fi, l’axe des ordonnées représente l’amplitude de la vitesse de phase déterminée. Les vitesses Vi déterminées par le procédé de la sont représentées par une courbe en pointillés 120. Les vitesses de phase pour une fréquence f i donnée telle que déterminée par la méthode MASW sont codées par une couleur d'autant plus sombre que l’amplitude de cette vitesse de phase est grande. Pour une même fréquence f i , la méthode MASW détermine plusieurs vitesses de phase correspondant, respectivement, au mode fondamental et aux autres modes d’ordre supérieur. Parmi ces différentes vitesses de phase déterminées, celle qui a l’amplitude la plus importante correspond à la vitesse de phase du mode fondamental. Comme le montre ce graphe, dans la plupart des cas, la vitesse V i déterminée par le procédé de la passe à l’emplacement le plus sombre du graphe de la . Cela indique que les vitesses de phase du mode fondamental déterminées par les deux méthodes différentes coïncident. Par contre, il existe des zones, entourées par des ellipses sur le graphe de la ou ce n’est pas le cas. Cette situation se rencontre en particulier pour les fréquences élevées, c’est-à-dire ici les fréquences supérieures à 240 Hz, mais aussi pour quelques fréquences plus basses. Dans ces zones, la vitesse de phase du mode fondamental, telle que déterminée par la méthode MASW, chute brusquement avant de remonter tout aussi brusquement. Par la suite, ces brusques chutes de la vitesse de phase du mode fondamental sont appelées "sauts de phase". Un tel saut de phase est entouré par un cercle au voisinage de 90 Hz sur la . À l’inverse, le procédé de la ne produit pas de tels sauts de phase. De plus, le procédé de la permet de déterminer la vitesse de phase du mode fondamental pour des fréquences bien plus élevées que ce qui est possible en mettant en œuvre la méthode MASW. À cause de cela, le procédé de la est considéré comme étant plus précis que les méthodes connues. Chapitre II : Variantes Variantes du dispositif de mesure : Dans une variante particulière, le dispositif est apte à mesurer, simultanément et long de plusieurs axes de mesure décalés angulairement les uns par rapport aux autres, le déplacement de la surface de la peau. Ainsi, avec un tel dispositif de mesure, il n’est pas nécessaire de le faire tourner autour de l’axe 22 de projection ou le nombre de rotations à réaliser est plus petit. Par exemple, pour cela, un tel dispositif de mesure comporte un réseau de capteurs optiques pour chaque axe de mesure. Dans un autre mode de réalisation, l’axe de mesure ne passe pas par le point d’impact 49 mais à côté de ce point d’impact. Le dispositif 10 peut aussi être réalisé à l'aide d'une caméra qui acquière des images de la surface 48 à une fréquence élevée. D’autres modes de réalisation du stimulateur 8 sont possibles. Par exemple, dans une variante, le jet d’air est remplacé par un jet d’un autre gaz, comme par exemple du dioxyde de carbone, ou par un jet de liquide comme de l’eau. Le stimulateur n’émet pas nécessairement un jet de fluide pour générer l’onde de cisaillement à la surface de la peau. Une telle onde de cisaillement peut également être générée par un stimulateur qui vient directement toucher la peau à l’aide d’un ustensile. Par exemple, le stimulateur peut être un marteau qui vient heurter la peau au niveau du point d’impact 49. Il peut aussi s’agir d’un projectile tel qu’une bille en caoutchouc ou similaire qui est projetée le long de l’axe 22 pour venir percuter la peau au niveau du point d’impact 49. Le dispositif de mesure n’est pas nécessairement un dispositif de mesure optique. Ceci est particulièrement vrai lorsque l'appareil 2 est appliqué sur un substrat de grande surface où les contraintes d'encombrement sont relâchées. Par exemple, les signaux u(x, t) peuvent aussi être mesurés en disposant directement sur chacun des points de mesure un capteur de déplacement. Il existe de très nombreux capteurs de déplacement connus susceptibles de convenir pour une telle application. Par exemple, ce capteur de déplacement peut être un accéléromètre. Dans un autre mode de réalisation, il est disposé sur chacun des points de mesure un élément rayonnant des ondes électromagnétiques ou réfléchissant des ondes électromagnétiques à une longueur d’onde particulière. Les capteurs du dispositif de mesure mesurent le déplacement de la surface en chacun de ces points de mesure à partir du rayonnement électromagnétique émis ou réfléchi au niveau du point de mesure. Variantes du procédé : D’autres méthodes sont possibles pour construire la vitesse de phase du mode fondamental en fonction de la profondeur. Par exemple, comme indiqué en regard de la , la méthode connue sous l’acronyme MASW peut être appliquée même si celle-ci est aujourd’hui considérée comme étant moins précise. Dans un autre mode de réalisation, c’est une variante de la méthode MASW qui est appliquée, cette variante ayant été modifiée pour réduire le problème des sauts de phase observés sur la . Enfin, d’autres méthodes que la méthode MASW ont été développées dans d’autres domaines techniques comme par exemple le domaine de la géophysique et peuvent être transposées ici à partir du moment où ces méthodes déterminent la vitesse de phase du mode fondamental. En variante, le pas d’échantillonnage des fréquences f i est le même sur tout l’intervalle d'analyse. Par exemple, dans l'exemple de mode de réalisation précédent, le pas d'échantillonnage est le même sur tout l'intervalle allant de 1 Hz à 1000 Hz. D’autres méthodes pour rechercher et identifier l’instant t i,m (x) sont possibles. Par exemple, les instants t i,m (x) sont recherchés sans limiter la recherche à un intervalle de temps prédéterminé. Dans ce cas, la recherche est réalisée sans tenir compte des instants t i-1,m (x) ou t i+1,m (x) ou t min (x). Ainsi, les étapes 82 et 84 peuvent être omises. L’erreur d’approximation peut être estimée différemment. En particulier, de nombreuses autres relations sont possibles pour calculer cette erreur d’approximation. Par exemple, la relation (5) suivante est utilisable à la place de la relation (1) : . Dans un autre mode de réalisation, le coefficient α n’est pas déterminé en fonction des signaux u(x, t) mesurés par le dispositif 10. Par exemple, dans un mode de réalisation simplifié, le coefficient α, utilisé pour une direction de mesure donnée, est fourni par l’utilisateur du dispositif 12. Au lieu de faire tourner sur lui-même le stimulateur 8 pour faire des mesures dans différentes directions de mesure, il est aussi possible de maintenir constante la direction de mesure et de déplacer le point d'impact le long d'une droite pour scanner progressivement une portion de la peau. Les mesures le long d'un même axe de mesure peuvent être répétées à différents instants pour voir l'évolution dans le temps de ces mesures. Par exemple, cela peut être appliqué pour mesurer l'évolution dans le temps des propriétés mécaniques de la peau suite à l'application d'un produit hydratant. Autres variantes : L’appareil 2 décrit ici peut être appliqué à d’autres matériaux viscoélastiques et anisotropes que la peau. Par exemple, il peut être appliqué à tous matériaux viscoélastiques similaires comme de la peau artificielle. Il peut aussi être appliqué à la peau des animaux, y compris la peau des animaux marins telle la peau des poissons. L’appareil 2 peut aussi être appliqué à d’autres matériaux viscoélastiques comme par exemple la peau des légumes ou des fruits. Lorsque l’appareil 2 est appliqué à d’autres matériaux viscoélastiques que la peau, l’intervalle dans lequel se situent les fréquences f i pour lesquelles l’unité 12 calcule la vitesse V i peut être différent de l’intervalle [1 Hz ; 1000 Hz]. De même, si le matériau viscoélastique sur lequel est appliqué l’appareil 2 n’est pas anisotrope, le coefficient α peut être calculé pour une première direction de mesure et la même valeur de ce coefficient α est ensuite utilisé pour d’autres directions de mesure décalées angulairement par rapport à la première direction de mesure. Dans ce cas, pour ces autres directions de mesure, l’opération 104 n’est pas réitérée. D’autres propriétés mécaniques que le module de Young de la peau à différentes profondeurs peuvent être déduites à partir de la vitesse de phase V i . Par exemple, la viscosité peut aussi être estimée à partir de la vitesse V i . Le procédé de la décrit pour construire la courbe de dispersion peut être mis en œuvre dans toutes applications non-invasives de mesure d’une grandeur physique représentative de propriétés mécaniques d’un substrat. En effet, comme expliqué précédemment, ce procédé permet d’obtenir une vitesse de phase du mode fondamental plus précise. Par exemple, le procédé de la peut aussi être mis en œuvre dans un appareil non-invasif de mesure de propriétés mécaniques en profondeur d’un substrat tel qu’un revêtement routier ou une structure multicouche quelconque. Dans le cas de tels substrats, il n’est pas nécessaire de mesurer le profil de vitesses dans plusieurs directions différentes angulairement décalées les unes par rapport aux autres. Il n’est pas non plus nécessaire de déterminer la valeur du coefficient α pour chacune des directions de mesure. Enfin, dans le cas d’autres substrats que la peau, le stimulateur 8 et le dispositif 10 de mesure sont adaptés à ce substrat. Par exemple, dans le cas d’un revêtement routier, le stimulateur 8 est formé d’une masse qui vient percuter la surface du revêtement routier. Le dispositif 10 mesure, dans le cas d'un revêtement routier, l’onde de cisaillement sur une distance typiquement plus grande que 7 mm. La détermination de la valeur du coefficient α en fonction des mesures du dispositif de mesure peut aussi être mise en œuvre dans tout autre appareil de mesure des propriétés mécaniques d’un substrat à partir d’un profil de vitesses de phase du mode fondamental. En effet, le calcul de la valeur du coefficient α comme décrit précédemment permet d’augmenter la précision de la conversion de la courbe de dispersion en un profil de vitesses. Chapitre III : Avantages des modes de réalisation décrits Le procédé de mesure décrit ici permet de mesurer la vitesse de phase du mode fondamental de l’onde de cisaillement à différentes profondeurs et non pas seulement en surface. Il permet donc de révéler des propriétés mécaniques de la peau à différentes profondeurs sous la peau et non pas seulement en surface. De plus, l’appareil 2 permet de révéler les propriétés mécaniques à différentes profondeurs tout en restant non-invasif, c’est-à-dire sans nécessité d’incision de la peau. Le fait de calculer la valeur du coefficient α à partir des signaux u(x, t) mesurés et plus précisément à partir de la courbe de dispersion, permet d’adapter automatiquement la valeur du coefficient α à la peau sur laquelle la mesure est faite et à la direction de mesure choisie. En effet, contrairement à d’autres substrats, il a été observé que la valeur du coefficient α varie de façon importante d’un être humain à un autre et varie également de façon importante en fonction de la direction de mesure choisie. Ainsi, le calcul automatique de la valeur du coefficient α permet donc d’améliorer la précision du profil de vitesses de phase. On obtient ainsi une meilleure précision sur la profondeur d'observation. Le fait de construire un profil de vitesses de phase pour plusieurs directions décalées angulairement les unes par rapport aux autres permet de générer une tomographie d’une propriété mécanique de la peau. Une telle tomographie permet notamment d’observer l’anisotropie de cette propriété mécanique en fonction de la direction de mesure. Le calcul de la vitesse de phase du mode fondamental à partir des instants t i , m (x) et des positions x où se produit le minimum du signal u i (x, t) permet de déterminer précisément la vitesse de phase du mode fondamental à la fréquence f i . Cela accroît donc la précision du profil de vitesses de phase construit, et donc la précision sur la mesure des propriétés mécaniques de la peau. Le fait d’identifier automatiquement la fréquence f min permet d’augmenter la reproductibilité du procédé de mesure car cette fréquence minimum n’est pas déterminée manuellement par l’utilisateur. De plus, lorsque la condition (1) est utilisée, cela permet de déterminer plus précisément la longueur d’onde λ max , ce qui améliore également la précision de la détermination de la valeur du coefficient α, et donc en fin de compte la précision du profil de vitesses de phase. Le fait d’identifier automatiquement la fréquence f max permet d’éviter de déterminer des vitesses de phase pour des composantes fréquentielles de l’onde de cisaillement qui n’existent pas ou qui sont négligeables. Cela améliore donc la précision du profil de vitesses de phase construit par l’appareil 2. Procédé non-invasif de mesure d'une grandeur physique représentative de l'élasticité de la peau, ce procédé comportant les étapes suivantes : a) la déformation (72) de la peau au niveau d'un point d'impact à l'aide d'un stimulateur pour générer une onde de cisaillement comportant des composantes à des fréquences différentes, ces différentes composantes se propageant à la surface de la peau et provoquant un déplacement de la surface de la peau, b) la mesure (74), à l'aide d'un dispositif de mesure, du déplacement de la surface de la peau au cours du temps en au moins trois points de mesure alignés les uns derrière les autres le long d'un axe de mesure, caractérisé en ce que le procédé comporte également les étapes suivantes : c) la détermination (86) des vitesses de phase du mode fondamental des différentes composantes de l'onde de cisaillement générée, le long de l'axe de mesure, à partir des mesures du dispositif de mesure, l'ensemble des paires, formées chacune par une fréquence f i et la vitesse de phase V i du mode fondamental déterminée pour cette fréquence f i , formant une courbe de dispersion du mode fondamental dans une direction de mesure parallèle à l'axe de mesure, où l'indice "i" est un numéro d'ordre de la fréquence f i et de la vitesse de phase V i , et d) la transformation (100) de cette courbe de dispersion en un profil de vitesses de phase en fonction de la profondeur, la vitesse de phase à une profondeur donnée étant une grandeur physique représentative de l'élasticité de la peau à cette profondeur. Procédé selon la revendication 1, dans lequel la transformation (100) de la courbe de dispersion en un profil de vitesses de phase dans la direction de mesure comporte les opérations suivantes : 1) convertir (102) chaque fréquence f i de la courbe de dispersion, en une longueur d'onde λ i correspondante à l'aide de la relation suivante : λ i = V i /f i , puis 2) calculer (104) la valeur d'un coefficient α pour cette direction de mesure à l'aide de la relation suivante : α = Z max / λ max , où : - Z max est égal à la moitié de la distance qui sépare les deux points de mesure les plus éloignés l'un de l'autre le long de l'axe de mesure, et - λ max est la plus grande des longueurs d'onde λ i obtenue à l'issue de l'exécution de l'opération 1) pour cette direction de mesure, 3) convertir (106) chaque longueur d'onde λ i obtenue à l'issue de l'opération 1) en une profondeur p i correspondante à l'aide de la relation suivante : p i = αλ i , où α est le coefficient calculé lors de l'opération 2). Procédé selon la revendication 2, dans lequel le procédé comporte l'exécution des étapes a), b), c) et d) pour au moins un premier et un second axes de mesure décalés angulairement l'un par rapport l'autre, ces premier et second axes de mesure passant par le même point d'impact. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape c) de détermination des vitesses de phase du mode fondamental comporte la réitération des opérations suivantes pour plusieurs fréquences f i différentes, : - pour chaque point de mesure : - le filtrage (88), à l'aide d'un filtre passe-bande centré sur la fréquence f i et dont la bande passante à -3dB est comprise entre les fréquences f i-1 et f i+1 , du signal u(x,t), mesuré au niveau du point de mesure de coordonnée x le long de l'axe de mesure par le dispositif de mesure, pour obtenir un signal filtré u i (x,t), - l'identification (90) d'un instant t i,m (x) auquel le signal filtré u i (x,t) passe par son minimum absolu, puis - le calcul (92) de la vitesse V i de propagation de ce minimum le long de l'axe de mesure à partir des instants t i,m (x) et des positions x auxquels ce minimum se produit, la vitesse V i ainsi calculée étant la vitesse de phase du mode fondamental à cette fréquence f i . Procédé selon la revendication 4, dans lequel l'étape c) comporte : - une opération d'identification automatique, parmi l'ensemble des fréquences f i , d'une fréquence minimale f min en dessous de laquelle les fréquences f i ne vérifient plus la condition (1) suivante : où : - µ i,1 et µ i,0 sont les coefficients de la droite, déterminée par la méthode des moindres carrés, approchant au mieux les points de coordonnées (x; t i,m (x)), - x p est égal à la position x du p-ième point de mesure compté à partir du premier point de mesure le plus proche du point d'impact, - P max est égal au nombre de points de mesure répartis le long de l'axe de mesure, - err max est une constante prédéterminée, - ensuite seule les fréquences f i supérieures ou égales à cette fréquence f min sont retenues pour former la courbe de dispersion. Procédé selon la revendication 5, dans lequel l'étape c) comporte : - une opération d'identification automatique, parmi l'ensemble des fréquences f i , d'une fréquence maximale f max au-dessus de laquelle les fréquences f i ne vérifient plus la condition (1), cette opération d'identification automatique de la fréquence f max étant réalisée en testant cette condition (1) pour plusieurs fréquences f i , ensuite - seule les fréquences f i inférieures ou égales à cette fréquence f max sont retenues pour former la courbe de dispersion. Procédé selon l'une quelconque des revendications précédentes, dans lequel toutes les fréquences f i sont comprises entre 1 Hz et 3000 Hz. Appareil non-invasif de mesure d'une grandeur physique représentative de l'élasticité de la peau, pour la mise en œuvre d'un procédé conforme à l'une quelconque des revendications précédentes, cet appareil comportant : - un stimulateur (8) apte à déformer la peau au niveau d'un point d'impact pour générer une onde de cisaillement comportant des composantes à des fréquences différentes, ces différentes composantes se propageant à la surface de la peau et provoquant un déplacement de la surface de la peau, - un dispositif (10) de mesure apte à mesurer le déplacement de la surface de la peau au cours du temps en au moins trois points de mesure alignés les uns derrière les autres le long d'un axe de mesure, caractérisé en ce que l'appareil comporte une unité (12) de traitement apte : - à déterminer les vitesses de phase du mode fondamental des différentes composantes de l'onde de cisaillement générée, le long de l'axe de mesure, à partir des mesures du dispositif de mesure, l'ensemble des paires, formées chacune par une fréquence f i et la vitesse de phase V i du mode fondamental déterminée pour cette fréquence f i , formant une courbe de dispersion du mode fondamental dans une direction de mesure parallèle à l'axe de mesure, où l'indice "i" est un numéro d'ordre de la fréquence f i et de la vitesse de phase V i , et - à transformer cette courbe de dispersion en un profil de vitesses de phase en fonction de la profondeur, la vitesse de phase à une profondeur donnée étant une grandeur physique représentative de l'élasticité de la peau à cette profondeur. Appareil selon la revendication 8, dans lequel l'unité (12) de traitement est configurée pour exécuter les opérations suivantes pour transformer la courbe de dispersion en un profil de vitesses de phase dans la direction de mesure: 1) convertir chaque fréquence f i de la courbe de dispersion, en une longueur d'onde λ i correspondante à l'aide de la relation suivante : λ i = V i /f i , puis 2) calculer la valeur d'un coefficient α pour cette direction de mesure à l'aide de la relation suivante : α = Z max / λ max , où : - Z max est égal à la moitié de la distance qui sépare les deux points de mesure les plus éloignés l'un de l'autre le long de l'axe de mesure parallèle, et - λ max est la plus grande des longueurs d'onde λ i obtenue à l'issue de l'exécution de l'opération 1) pour cette direction de mesure, 3) convertir chaque longueur d'onde λ i obtenue à l'issue de l'opération 1) en une profondeur p i correspondante à l'aide de la relation suivante : p i = α * λ i , où α est le coefficient calculé lors de l'opération 2). Appareil selon la revendication 8 ou 9, dans lequel : - le dispositif de mesure comporte un réseau (62) de capteurs apte chacun à mesurer l'amplitude de la déformation de la surface de la peau en un point de mesure respectif, ce réseau comportant au moins trois capteurs qui mesurent chacun le déplacement de la surface de la peau en trois points de mesure respectifs alignés les uns derrière les autres le long d'un axe de mesure, et - l'appareil comporte un bras (4) articulé sur lequel est monté le réseau de capteurs, ce bras articulé étant apte à faire tourner le réseau de capteurs d'un angle prédéterminé autour d'un axe de rotation, pour aligner l'axe de mesure du réseau de capteur sur un premier axe de mesure et, en alternance, sur un second axe de mesure décalé angulairement par rapport au premier axe de mesure. Appareil selon la revendication 10, dans lequel : - le réseau de capteurs comporte : - une ligne (62) de capteurs optiques qui captent la lumière réfléchie par chaque point de mesure, et - un microprocesseur (64) configuré pour déterminer le déplacement de la surface de la peau au niveau de chaque point de mesure éclairé à partir de la lumière réfléchie captée par la ligne de capteurs optiques, - le dispositif de mesure comporte un émetteur (60) d'un faisceau lumineux qui éclaire chaque point de mesure aligné le long de l'axe de mesure. Appareil selon l’une quelconque des revendications précédentes, dans lequel le stimulateur (8) est apte à projeter sur la surface de la peau un jet de fluide qui provoque la déformation de la peau au niveau du point d'impact.