Brevet FR3069150A1 — publié le 20190125

DOMAINE DE L'INVENTION La présente invention concerne un procédé de caractérisation d'un os. ETAT DE LA TECHNIQUE L'utilisation d'ondes pour caractériser un corps vivant est largement connue. L'imagerie par tomographie à rayons X est par exemple une technique d'imagerie médicale utilisant des rayons X traversant le corps à imager. Cette technique présente toutefois l'inconvénient d'exposer le corps à un rayonnement ionisant potentiellement dangereux. Encore aujourd'hui, on évite d'exposer le corps de jeunes enfants à de tels rayons X sauf absolue nécessité. D'autres techniques connues ayant pour but de caractériser un corps vivant utilisent des ondes ultrasonores, qui sont des ondes moins dangereuses que les rayons X. Des ondes ultrasonores sont conventionnellement émises par un réseau d'émetteurs récepteurs, et leurs échos sur un corps à caractériser sont reçus par le même réseau ou un autre réseau après une certaine durée de propagation entre l'émission et la réception. Des images montrant une coupe du corps dans lequel les ondes ultrasonores se sont propagées peuvent ensuite être construites sur la base des signaux d'écho reçus par le (ou les) réseau(x) d'émetteurs-récepteurs. Au cours de la construction de telles images, il est fait l'hypothèse que le corps est un milieu homogène, et que, par conséquent, la vitesse du son est uniforme dans le corps étudié. La vitesse du son choisie est généralement une vitesse du son moyenne dans du tissu biologique non osseux (peau ou muscle par exemple), qui est généralement de l'ordre de 1540 mètres par seconde à une erreur près de 5% à 10%. Les images obtenues sur la base de cette hypothèse présente ainsi une qualité satisfaisante dans des régions d'intérêt montrant du tissu biologique non osseux. Or, la vitesse du son dans un os (généralement comprise entre 2800 mètres par secondes et 4200 mètres par secondes pour une onde de compression dans l'os cortical) est très supérieure à la vitesse du son dans un tissu biologique non osseux. Des images obtenues sur la base d'une hypothèse de vitesse uniforme dont la valeur correspond à une vitesse du son dans du tissu biologique non-osseux présentent une mauvaise qualité dans des régions d'intérêt montrant l'os. Cette mauvaise qualité se traduit typiquement par une faible intensité et un rendu flou au niveau de l'os. Pour cette raison, une idée reçue largement répandue dans le milieu médical est que les ondes ultrasonores ne - pénètrent pas facilement » dans de l'os. 2 La différence entre la vitesse du son dans un os et dans les tissus biologiques non osseux entourant l'os crée un effet de réfraction. Il n'est donc pas possible de reconstruire des images de bonne qualité en se fondant sur une valeur prédéterminée de vitesse de son dans l'os, comme on a déjà pu le faire avec une vitesse de son dans un tissu biologique non 5 osseux. La vitesse du son dans un os dépend de plusieurs facteurs. Premièrement, la vitesse du son dans un os varie d'un individu à un autre. C'est d'ailleurs la raison pour laquelle elle constitue une information très utile pour diagnostiquer ultérieurement des troubles tels que l'ostéoporose. 10 Deuxièmement, la vitesse du son dans un os n'est pas la même dans toutes les directions. L'os est en effet un milieu anisotrope. Ceci est dû notamment au fait que l'os dit cortical comprend des canaux s'étendant parallèlement à l'axe longitudinal d'un os long (tibia par exemple), pour héberger des vaisseaux sanguins. Ainsi, une onde sonore de compression se propage dans un os suivant une direction parallèle à l'axe longitudinal d'un 15 os long plus rapidement que suivant une autre direction. En définitive, déterminer précisément la vitesse de son constitue un enjeu majeur pour caractériser un os. A cet égard, il a déjà été proposé un procédé de détermination de la vitesse du son dans un os parallèlement à l'axe longitudinal de l'os. Au cours de ce procédé, des ondes 20 ultrasonores sont émises par un réseau d'émetteurs-récepteurs d'ondes ultrasonores alignés parallèlement à l'axe longitudinal d'un os long. Sur la base de signaux d'écho reçus par les récepteurs, on détermine une vitesse d'une onde guidée par l'os, dite onde de tête, se propageant le long de la surface extérieure de l'os. Cette vitesse est déterminée assez facilement en faisant l'hypothèse que la relation entre l'instant de réception d'une onde 25 ultrasonore par un des récepteurs et la distance séparant ce récepteur de l'émetteur de la même onde est une fonction linéaire. Cependant, ce procédé ne peut pas être utilisé pour déterminer une vitesse du son dans l'os suivant une autre direction qu'une direction parallèle à l'axe longitudinal d'un os long. En conséquence, il ne permet de caractériser un os que de façon partielle. 30 EXPOSE DE L'INVENTION Un but de l'invention est de proposer un procédé qui puisse être utilisé pour caractériser de façon plus complète un os. Il est dès lors proposé un procédé de caractérisation d'un os, le procédé comprenant 35 des étapes de : 3 - émission de premières ondes ultrasonores vers un corps comprenant l'os et du tissu biologique non osseux entourant l'os, - réception de premiers signaux d'écho des premières ondes ultrasonores émises, - détermination d'une vitesse du son dans le tissu biologique non osseux suivant une 5 première direction, la détermination comprenant des sous-étapes de : o pour plusieurs premières valeurs candidates prédéterminées, construction d'une première image montrant le tissu biologique non osseux et le périoste de l'os, à partir des premiers signaux d'écho et sous l'hypothèse que la vitesse du son dans le tissu biologique non osseux suivant la première 10 direction est égale à la première valeur candidate, o pour chaque première image, calcul d'au moins une première métrique indicative d'une qualité de mise au point du périoste et/ou du tissu biologique non osseux entourant l'os dans la première image, o sélection d'une des premières valeurs candidates en tant que vitesse du 15 son dans le tissu biologique non osseux suivant la première direction, en fonction des premières métriques, - localisation d'une première courbe de démarcation entre le tissu biologique non osseux et le périoste dans l'une des premières images, - détermination d'une vitesse du son dans l'os suivant la première direction, la 20 détermination comprenant des sous-étapes de : o pour plusieurs deuxièmes valeurs candidates prédéterminées, construction d'une deuxième image montrant du tissu osseux cortical de l'os et l'endoste de l'os, à partir des premiers signaux d'écho, de la vitesse du son dans le tissu biologique non osseux suivant la première direction déterminée, et de 25 la première courbe de démarcation, et sous l'hypothèse que la vitesse du son dans l'os suivant la première direction est égale à la deuxième valeur candidate, o pour chaque deuxième image, calcul d'au moins une deuxième métrique indicative d'une qualité de mise au point du tissu osseux cortical de l'os 30 et/ou de l'endoste de l'os dans la deuxième image, o sélection d'une des deuxièmes valeurs candidates en tant que vitesse du son dans l'os suivant la première direction, en fonction des deuxièmes métriques. 35 4 Le procédé proposé peut également comprendre les caractéristiques ou étapes optionnelles suivantes, prises seules ou en combinaison lorsque cela est techniquement possible. La construction d'une première image à l'aide d'une première valeur candidate peut 5 comprendre les étapes suivantes mises en oeuvre pour au moins un point du tissu biologique non osseux : - estimation de premières trajectoires de premières ondes ultrasonores émises par des émetteurs, puis passées par le point du tissu biologique non osseux, puis reçues par des récepteurs, à partir des premiers signaux d'écho, sous l'hypothèse que la 10 vitesse du son dans le tissu biologique non osseux suivant la première direction est égale à la première valeur candidate, - calcul de durées de propagation des ondes ultrasonores via les premières trajectoires estimées, - calcul d'une intensité d'un pixel de la première image au point du tissu biologique 15 non osseux, à partir des durées de propagation, des premiers signaux d'écho et de positions des émetteurs et des récepteurs. La construction d'une deuxième image à l'aide d'une deuxième valeur candidate comprend les étapes suivantes mises en oeuvre pour au moins un point de l'os : - estimation de deuxièmes trajectoires de premières ondes ultrasonores émises par 20 des émetteurs, puis passées par le point de l'os, puis reçues par des récepteurs, à partir des signaux d'écho, à partir de la vitesse du son dans le tissu biologique non osseux suivant la première direction déterminée, de la première courbe de démarcation, et sous l'hypothèse que la vitesse du son dans l'os suivant la première direction est égale à la deuxième valeur candidate, 25 - calcul de durées de propagation des ondes ultrasonores via les deuxièmes trajectoires estimées, - calcul d'une intensité d'un pixel de la deuxième image au point de l'os, à partir des durées de propagation, des premiers signaux d'écho et de positions des émetteurs et des récepteurs. 30 La localisation de la première courbe de démarcation peut être mise en oeuvre dans la première image qui a été construite à l'aide de la première valeur sélectionnée en tant que vitesse du son dans le tissu biologique non osseux suivant la première direction. Les premières ondes ultrasonores peuvent être des ondes émises par des émetteurs et les signaux d'échos reçus par des récepteurs alignés le long d'un axe perpendiculaire à un 5 axe longitudinal de l'os, et dans lequel la première direction est une direction perpendiculaire à l'axe longitudinal de l'os. Le procédé peut comprendre un affichage de la première image construite à l'aide de la première valeur sélectionnée, et/ou de la deuxième image construite à l'aide de la 5 deuxième valeur sélectionnée. Le procédé peut en outre comprendre des étapes de : - émission de deuxièmes ondes ultrasonores vers le corps, - réception de signaux d'écho des deuxièmes ondes ultrasonores émises, dits deuxièmes signaux d'écho, 10 - détermination d'une vitesse du son dans le tissu biologique non osseux suivant une deuxième direction différente de la première direction, comprenant des sous-étapes de : o pour plusieurs troisièmes valeurs candidates, construction d'une troisième image montrant le tissu biologique non osseux et le périoste de l'os, à partir 15 des deuxièmes signaux d'écho et sous l'hypothèse que la vitesse du son dans le tissu biologique non osseux suivant la deuxième direction est égale à la deuxième valeur candidate, o pour chaque troisième image, calcul d'au moins une troisième métrique indicative d'une qualité de mise au point du périoste dans la troisième 20 image, o sélection d'une des troisièmes valeurs candidates en tant que vitesse du son dans le tissu biologique non osseux suivant la deuxième direction, en fonction des troisièmes métriques, - détermination d'une vitesse du son dans l'os suivant la deuxième direction, à l'aide 25 des deuxièmes signaux d'écho, - localisation d'une deuxième courbe de démarcation entre le tissu biologique non osseux et le périoste dans l'une des troisièmes images, - détermination d'un paramètre d'anisotropie de l'os susceptible d'être utilisé par une fonction prédéterminée en combinaison avec la vitesse du son dans l'os suivant 30 la première direction et avec la vitesse du son dans l'os suivant la deuxième direction pour calculer une vitesse du son dans l'os suivant une direction quelconque, la détermination du paramètre d'anisotropie comprenant des sous-étapes de : o pour plusieurs quatrièmes valeurs candidates prédéterminées, construction 35 d'une quatrième image montrant du tissu osseux cortical et l'endoste de 6 l'os, à partir des deuxièmes signaux d'écho, de la vitesse de son dans l'os suivant la première direction, de la deuxième courbe de démarcation, de la vitesse de son dans l'os suivant la deuxième direction, de la fonction prédéterminée, et sous l'hypothèse que le paramètre d'anisotropie de l'os 5 est égal à la quatrième valeur candidate, o pour chaque quatrième image, calcul d'une quatrième métrique indicative d'une qualité de mise au point de l'endoste et/ou du tissu osseux cortical de l'os dans la quatrième image, o sélection d'une des quatrièmes valeurs candidates en tant que paramètre 10 d'anisotropie de l'os, en fonction des quatrièmes métriques. La construction d'une quatrième image peut comprendre les étapes suivantes mises en oeuvre pour au moins un point de l'os : - estimation de troisièmes trajectoires de deuxièmes ondes ultrasonores émises par des émetteurs, puis passées par le point de l'os, puis reçues par des récepteurs, à 15 partir des deuxièmes signaux d'écho, de la vitesse de son dans l'os suivant la première direction déterminée, de la vitesse du son dans l'os suivant la deuxième direction déterminée, de la deuxième courbe de démarcation, de la fonction prédéterminée et sous l'hypothèse que le paramètre d'anisotropie de l'os est égale à la quatrième valeur candidate, 20 - calcul de durées de propagation des ondes ultrasonores via les troisièmes trajectoires estimées, - calcul d'une intensité d'un pixel de la quatrième image au point de l'os, à partir des durées de propagation et des deuxièmes signaux d'écho et de positions des émetteurs et des récepteurs. 25 La localisation de la deuxième courbe de démarcation peut être mise en oeuvre dans la troisième image ayant été construite à l'aide de la troisième valeur sélectionnée en tant que vitesse du son dans le tissu biologique non osseux suivant la deuxième direction. Le procédé peut comprendre un affichage de la troisième image construite à l'aide de la troisième valeur sélectionnée, et/ou de la quatrième image construite à l'aide de la 30 quatrième valeur sélectionnée et de la vitesse du son dans l'os suivant la deuxième direction. Les deuxièmes ondes ultrasonores peuvent être des ondes émises par des émetteurs et les deuxièmes signaux d'écho sont reçus par des récepteurs alignés le long d'un axe compris dans un plan parallèle à un axe longitudinal de l'os, et dans lequel la deuxième 35 direction est une direction parallèle à l'axe longitudinal de l'os. 7 Chaque vitesse de son dans l'os déterminée peut être une vitesse de propagation d'ondes de compression ou une vitesse de propagation d'ondes de cisaillement. Chaque vitesse de son dans l'os déterminée est une vitesse de propagation d'ondes de 5 compression, et la fonction prédéterminée être une fonction V de la forme : V(0) = V axial_os (Vaxial_os Vradial_os) Lig sin2(0) cos2(0) + cos4(0)] dans lequel - V radial_os est la vitesse de propagation d'ondes de compression dans l'os suivant la première direction déterminée, 10 - V axial_os est la vitesse de propagation d'ondes de compression dans l'os suivant la deuxième direction déterminée, - est le paramètre d'anisotropie de l'os, - 0 est un angle, - V(0) est une vitesse de propagation d'ondes de compression dans l'os suivant une 15 direction formant l'angle 0 avec la première direction. La vitesse de son dans l'os suivant la première direction déterminée peut être une vitesse de propagation d'ondes de cisaillement, et la fonction prédéterminée être une fonction V de la forme : V(0) = Vcis_axial_os [1 + sin2(0) cos2(0)] 20 dans lequel - V cis_axial_os est la vitesse de propagation d'ondes de cisaillement dans l'os suivant la première direction déterminée, - est le paramètre d'anisotropie de l'os, - 0 est un angle, 25 - V(0) est une vitesse de propagation d'ondes de cisaillement dans l'os suivant une direction formant l'angle 0 avec la première direction. Il peut être déterminé, conjointement à la vitesse du son dans l'os suivant la première direction, un paramètre d'anisotropie de l'os susceptible d'être utilisé par une fonction prédéterminée en combinaison avec la vitesse du son dans l'os suivant la première 30 direction déterminée pour calculer une vitesse du son dans l'os suivant une direction quelconque, ladite détermination conjointe comprenant des sous-étapes de : - pour plusieurs paires de deuxièmes et troisièmes valeurs candidates prédéterminées, construction d'une deuxième image montrant du tissu osseux cortical de l'os et l'endoste de l'os, à partir des premiers signaux d'écho, de la 35 vitesse du son dans le tissu biologique suivant la première direction déterminée, et 8 de la première courbe de démarcation, sous l'hypothèse que la vitesse du son dans l'os suivant la première direction est égale à la deuxième valeur candidate, et sous l'hypothèse que le paramètre d'anisotropie de l'os est égal à la troisième valeur candidate, 5 - pour chaque deuxième image, calcul d'au moins une deuxième métrique indicative d'une qualité de mise au point du tissu osseux cortical de l'os et/ou de l'endoste de l'os dans la deuxième image, - sélection conjointe d'une deuxième valeur candidate en tant que vitesse du son dans l'os suivant la première direction, et d'une troisième valeur candidate en tant 10 que paramètre d'anisotropie de l'os, la sélection étant mise en oeuvre en fonction des deuxièmes métriques. DESCRIPTION DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de la 15 description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : - La figure 1 est une vue schématique en perspective d'un corps comprenant de l'os. - La figure 2 est une représentation schématique d'un système de caractérisation d'un os, selon un mode de réalisation de l'invention. 20 - La figure 3 montre certaines étapes d'un procédé de caractérisation d'un os, selon un premier mode de réalisation de l'invention. - Les figures 4 et 6 détaillent des sous-étapes d'étapes du procédé selon le premier mode de réalisation de l'invention. - Les figures 5a et 5b sont deux exemples d'images construites au cours de la mise en 25 oeuvre du procédé selon le premier mode de réalisation de l'invention. - La figure 7 est un plan de coupe longitudinal d'un corps comprenant de l'os, montrant également le trajet d'une onde ultrasonore dans ce corps. - La figure 8 est un ensemble de courbes, chaque courbe montrant l'évolution d'une métrique calculée au cours de la mise en oeuvre du procédé selon le premier mode 30 de réalisation de l'invention, en fonction d'une valeur candidate de vitesse de propagation d'ondes de compression radiales dans un os. - La figure 9 est montre d'autres étapes du procédé de caractérisation d'un os selon le premier mode de réalisation de l'invention. - Les figures 10 et 11 détaillent des sous-étapes d'étapes représentées en figure 9. 9 - La figure 12 est un ensemble de courbes, chaque courbe montrant l'évolution d'une métrique calculée lors de la mise en oeuvre du procédé selon le premier mode de réalisation de l'invention, en fonction d'une valeur candidate pour un paramètre d'anisotropie d'un os. 5 - La figure 13 est un organigramme d'étapes d'un procédé de caractérisation d'un os, selon un deuxième mode de réalisation de l'invention. - La figure 14 détaille des sous-étapes d'une étape du procédé selon le deuxième mode de réalisation de l'invention. Sur l'ensemble des figures, les éléments similaires portent des références identiques. 10 DESCRIPTION DETAILLEE DE MODES DE REALISATION DE L'INVENTION On a illustré de manière schématique en figure 1 un exemple de partie d'un corps C comprenant un os B et du tissu biologique non osseux T. L'os B s'étend le long d'un axe longitudinal X. L'os B comprend notamment de la 15 moelle, un endoste E s'étendant autour de la moelle (représenté en pointillés sur la figure 1), et un périoste PE qui entoure l'endoste E. L'os B comprend par ailleurs du tissu osseux cortical se trouvant entre l'endoste E et le périoste PE. Le tissu biologique non osseux T s'étend autour de l'os, et plus spécifiquement autour du périoste PE avec lequel il est en contact. 20 Le tissu biologique non osseux T comprend de la chair voire de la peau entourant la chair. De façon bien connue, des ondes ultrasonores peuvent se propager dans un tel corps C, A cet égard, on a représenté en figure 1 plusieurs vitesses propagation d'ondes sur la 25 figure 1, qui sont les suivantes : - Vradial_os Vitesse de propagation d'une onde de compression dans l'os B suivant une direction perpendiculaire à l'axe longitudinal X de l'os B (on parlera alors de direction radiale », et d'ondes de compression radiales » ou transverses ») - Vaxial_os Vitesse de propagation d'une onde de compression dans l'os B suivant une 30 direction parallèle à l'axe longitudinal X (on parlera alors de direction axiale » et d'ondes de compression axiales ») - V radial_tissu : vitesse de propagation d'une onde de compression dans le tissu biologique non osseux T suivant une direction perpendiculaire à un axe longitudinal d'un os, 10 - V - axial_tissu : vitesse de propagation d'une onde de compression dans le tissu biologique non osseux T suivant une direction parallèle à l'axe longitudinal X de l'os B. La propagation d'ondes de compression dans un os suivant une direction quelconque 5 peut être calculée au moyen d'une fonction prédéterminée V combinant un paramètre d'anisotrope )0 et des vitesses Vaxial_as et Vradial_os- La fonction prédéterminée V est alors typiquement de la forme suivante : V(e) = Vaxial_os (Vaxial_os Vradial _os) sin2(0) où 10 - 0 est un angle quelconque, - )0 est un paramètre d'anisotropie de l'os, - V(0) est une vitesse du son dans l'os B suivant une direction formant l'angle 0 avec la première direction. Cette fonction V se fonde sur l'hypothèse relativement réaliste que l'os B est isotrope 15 transverse, c'est-à-dire isotrope dans un plan perpendiculaire à son axe longitudinal X. En référence à la figure 2, un système 1 de caractérisation comprend une sonde 2 ultrasonore, un dispositif de traitement 4 de signaux d'écho acquis par la sonde 2, voire un écran d'affichage 12. La sonde 2 ultrasonore, connue en elle-même, comprend au moins un réseau 20 d'émetteurs-récepteurs 6 alignés le long d'un axe Y. De façon conventionnelle, la sonde 2 comprend une lentille en silicone (non illustrée) agencée devant la rangée d'émetteurs-récepteurs 6. Chaque émetteur-récepteur 6 est adapté pour émettre des ondes ultrasonores. Chaque émetteur-récepteur 6 est également adapté pour acquérir des signaux d'écho 25 d'ondes ultrasonores émises par un quelconque autre émetteur-récepteur 6. Les positions relatives des émetteurs-récepteurs 6 sont prédéterminées. Typiquement, les émetteurs-récepteurs 6 sont distants d'un pas constant le long de l'axe de la sonde 2. Par ailleurs, le dispositif de traitement 4 de signaux d'écho comprend 30 conventionnellement au moins un processeur 8 et une mémoire 10. Le processeur 8 est configure pour exécuter des calculs, et en particulier un algorithme de traitement d'images dont le fonctionnement sera détaillé dans la suite. La mémoire 10 stocke des données prédéterminées. Ces données ne sont pas propres au corps d'un individu mais sont des données génériques applicables à tout individu d'une 35 population. cos2(0) + cos4(0)] 11 Les données prédéterminées comprennent un ensemble de premières valeurs candidates prédéterminées pour la vitesse Vradial_tis su - Les premières valeurs candidates prédéterminées sont typiquement comprises entre 1400 et 1700 mètres par seconde. Les données prédéterminées comprennent par ailleurs un ensemble de deuxièmes 5 valeurs candidates de vitesse pour la vitesse Vradial_o s - Les deuxièmes valeurs candidates prédéterminées sont typiquement comprises entre 2600 et 3600 mètres par seconde. Les données prédéterminées comprennent un ensemble de troisièmes valeurs candidates prédéterminées pour la vitesse Vaxial_tissu - Les troisièmes valeurs candidates prédéterminées sont typiquement comprises entre 3400 et 4200 mètres par seconde. 10 Les données prédéterminées comprennent par ailleurs un ensemble de quatrièmes valeurs candidates pour le paramètre fi' d'anisotropie de l'os B. Comme on le verra dans la suite, l'algorithme mis en oeuvre par le processeur 8 fait la supposition que la vitesse de propagation dans un os B est régie par la fonction V décrite plus haut, laquelle dépend notamment du paramètre fi' d'anisotropie. 15 En référence à l'organigramme de la figure 3, un procédé de caractérisation de l'os B représenté en figure 1 au moyen du système de caractérisation 1 comprend les étapes suivantes. La sonde 2 est positionnée à proximité du corps C dans une première position. Dans la première position, les émetteurs-récepteurs 6 de la sonde 2 sont alignés 20 sensiblement perpendiculairement à l'axe longitudinal de l'os B. Autrement dit l'axe Y de la sonde 2 est perpendiculaire à l'axe X de l'os B dans la première position. La sonde 2 émet des premières ondes ultrasonores en direction du corps C suivant une première direction (étape 100). Les premières ondes sont par exemple des ondes de compression radiales. Autrement 25 dit, la première direction est une direction perpendiculaire à l'axe longitudinal de l'os B. Les premières ondes ultrasonores pénètrent dans le corps C et sont réfléchies par ce dernier à différents niveaux de profondeur dans le corps C. Des signaux d'écho de ces ondes, dits premiers signaux d'écho, sont ainsi acquis par les émetteurs récepteurs de la sonde 2 (étape 102). Bien entendu, une onde émise par un 30 émetteur-récepteur 6 d'indice i peut parfaitement donner naissance à un signal d'écho reçu par un autre émetteur-récepteur 6 d'indice j. Les premiers signaux d'écho sont numérisés, transmis au dispositif de traitement 4 et mémorisés dans la mémoire 10 sous une forme connue de l'état de la technique. 12 Le processeur 8 détermine la vitesse Vradial_tissu sur la base des premiers signaux d'écho reçus et sur la base des premières valeurs prédéterminées mémorisées dans la mémoire 10 (étape 104). En référence à la figure 4, la détermination 104 de la vitesse Vradial_tissu comprend 5 les sous-étapes suivantes. Le processeur 8 construit une image, dite première image, à partir des premiers signaux d'écho, sous l'hypothèse que la vitesse Vradial_tissu est égale à une première valeur candidate (étape 200). La première image est constituée d'une grille de pixels, chaque pixel étant défini par 10 une position dans la grille et par une intensité I, cette intensité étant typiquement représentative d'un niveau de gris. La première image représente par ailleurs une vue en coupe du corps C dans un plan dans lequel les premières ondes ultrasonores se sont propagées (ce plan de coupe étant parallèle à l'axe de la sonde 2). Chaque point de ce plan de coupe sera ainsi montré dans un pixel de la première image. 15 La construction de la première image est réalisée au moyen de la méthode de migration de Kirchhoff ou de la méthode dite - de mise au point totale » (- Total Focusing Method » en anglais). Ces méthodes sont connues en elles-mêmes, mais avantageusement appliquées aux premiers signaux ultrasonores que l'on a précédemment acquis. La construction 200 de la première image au moyen d'une de ces méthodes connues 20 comprend les sous-étapes suivantes. Pour un point P donné du corps C étudié, le processeur 8 estime les trajectoires de premières ondes ultrasonores émises par les émetteurs, puis passées par le point P, puis reçues par des récepteurs, à partir des premiers signaux d'écho, sous l'hypothèse que la vitesse Vradtal_os est égale à une des première valeur candidate. 25 Les ondes ultrasonores passées par le point P ont été chacune émise par un émetteur d'indice i, dont la position est connue le long de l'axe Y de la sonde 2, et reçue par un récepteur d'indice j, dont la position est également connue le long de l'axe Y de la sonde 2. On a donc au plus autant de premiers signaux d'écho d'ondes passées par le point P que de paires (i, j) d'indices d'émetteurs/récepteurs dans la sonde 2 (donc au plus M x N 30 signaux si M est le nombre d'émetteurs utilisés et N le nombre de récepteurs utilisés). L'estimation des trajectoires des ondes ultrasonores est mise en oeuvre en exploitant le principe de Fermat, selon lequel il est supposé qu'une onde se propage rectilignement dans un milieu homogène. Ici, le corps C comprenant du tissu biologique non osseux et de l'os B est considéré dans son ensemble comme un milieu homogène. Il est fait l'hypothèse 35 que la vitesse de propagation d'une onde est égale à la première valeur candidate 13 considérée, qui correspond à une valeur crédible de vitesse de propagation du son dans du tissu biologique non osseux T. Le processeur 8 calcule ensuite des durées de propagation des ondes passées par le point P via les trajectoires estimées. 5 Une durée de propagation se décompose en une durée de propagation tT(i, P) depuis l'émetteur d'indice i jusqu'au point P, et une durée de propagation tR(jP) depuis le point P jusqu'au récepteur d'indice j. Le processeur 8 calcule ensuite une intensité d'un pixel de la première image au point P considéré, à partir des durées de propagation estimées, des premiers signaux d'écho et 10 des positions des émetteurs et des récepteurs. L'intensité I du point P est typiquement calculée via la formule ci-dessous : M N 1(P) j)xD(t = tT(i,P)+ tR(j, j)) 1=1 j=1 dans laquelle : - D(t = tT(i,P)+ tR(j, j)) désigne une donnée représentative d'un signal d'écho 15 reçu à l'instant t par le récepteur d'indice j, l'écho provenant d'une onde initialement émise par l'émetteur d'indice i, - W(P,i,j) désigne un poids obtenu par application d'une fonction de pondération W prédéterminée. Typiquement, la fonction de pondération W est une fonction fenêtre d'observation 20 (également appelée fenêtre de pondération ou d'apodisation dans la littérature). On a W(P,i,j) = 1 si l'angle du segment de trajet retour de l'onde allant du point P au récepteur d'indice j, par rapport à une direction normale à un plan d'émission/réception de la sonde 2, est inférieur à un seuil angulaire prédéterminé, et l'on a W(P,i,j) = 0 sinon. Ce seuil angulaire est par exemple fixé à 50 degrés (cet angle correspondant à une perte de 25 sensibilité d'un récepteur de l'ordre de 50%). En répétant les sous étapes qui précèdent en plusieurs points P, on peut reconstruire une première image. L'étape de construction 200 d'une première image est alors terminée. L'étape de construction 200 d'une première image est répétée pour chaque première valeurs mémorisée dans la mémoire 10. 30 On a représenté en figure 5a un exemple de première image que l'on peut obtenir par la mise en oeuvre de l'étape de construction 200. Comme on peut le constater, cette première image est relativement riche en information au niveau du tissu biologique non osseux T entourant l'os. En revanche, cette première image est pauvre en information à l'intérieur de l'os. La première image présente notamment une zone blanchâtre floue au 14 niveau du périoste PE de l'os B. Ceci est dû au fait que le processeur 8 a pris comme vitesse de propagation de référence une première valeur candidate qui correspond à une valeur crédible de vitesse de propagation du son dans du tissu biologique non osseux T, mais qui ne correspond pas à une valeur crédible de vitesse de propagation du son dans un 5 os, comme cela a été évoqué dans l'introduction du présent texte. De retour à la figure 4, pour une première image construite au cours de l'étape 200, le processeur 8 calcule une première métrique représentative d'une qualité de mise au point dans une région d'intérêt de la première image (étape 202). La région d'intérêt choisie est typiquement une région montrant le périoste PE et/ou du tissu biologique non 10 osseux T entourant l'os B. La première métrique est de préférence fonction d'une intensité moyenne et/ou d'un contraste moyen dans la région d'intérêt considérée dans la première image. La première métrique est typiquement l'une ou une combinaison des métriques suivantes, connues de l'état de la technique : 15 - La métrique d'énergie spectrale latérale décrite dans le document - Sound Speed Correction in ultrasound image » par D. Napolitano, C.Chou, G. McLaughlin et al., publié en 2006, - la métrique dite de - critère de netteté par la méthode de Brenner » ou la métrique dite de - critère de netteté par la méthode de Tenenbaum » ou la métrique appelée 20 - critère de netteté par la méthode de la variance », toutes décrites dans le document - Automatic sound speed selection in photoacoustic image reconstruction using an autofocus approach », par B. Treeby, T. Varslot, E. Zhang et al., publié en 2011. Le processeur 8 répète l'étape 202 pour chaque première image construite au cours 25 de l'étape 200. Le processeur 8 génère ainsi autant de premières métriques que de premières images, et que de premières valeurs candidates prédéterminées. Le processeur 8 sélectionne ensuite en tant que valeur définitive pour la vitesse Vradtal _tissu une valeur optimale parmi les premières valeurs candidates utilisées pour produire les premières images (étape 204). Le processeur 8 se base pour cela sur les 30 premières métriques. La première valeur candidate sélectionnée en tant que vitesse Vradtai _tissu au cours de l'étape 204 est celle qui a servi de donnée d'entrée pour produire une première image dont la métrique associée est indicative d'une qualité de mise au point dans la région d'intérêt considérée qui est maximale parmi toutes les premières métriques calculées. Typiquement, 35 lorsque l'une des méthodes listées ci-dessus est utilisée pour calculer la première 15 métrique, on recherche la première métrique de valeur maximale parmi les premières métriques calculées. On a ainsi déterminé la vitesse Vradial_tissu réelle dans le tissu organique du corps C étudié. L'étape 104 est alors terminée. 5 De retour à la figure 3, le processeur 8 localise par ailleurs une première courbe de démarcation entre le tissu biologique non osseux et l'os B dans l'une des premières images construites (étape 106). Cette courbe de démarcation se situe en pratique au niveau d'une zone de l'image de forte intensité, correspondant au périoste PE. Très préférentiellement, la localisation 106 est mise en oeuvre dans la première 10 image construite à l'aide de la première valeur candidate sélectionnée en tant que vitesse Vradial_os- La localisation est alors beaucoup plus précise du fait de la haute qualité de mise au point de cette première image parmi toutes celles qui ont été construites par le processeur 8. La localisation 106 comprend conventionnement les étapes suivantes : 15 - la première image fait l'objet d'une segmentation, de sorte à identifier un groupe de pixels de la première image montrant la frontière entre le tissu biologique non osseux et l'os B (cette segmentation comprenant par exemple la mise en oeuvre d'un algorithme de Djikstra connu de l'état de la technique) ; - ce groupe de pixels est approximé en une courbe de démarcation définie par un 20 polynôme, par exemple une parabole. Le processeur 8 détermine ensuite la vitesse Vradial_os (étape 108) en exploitant avantageusement la première courbe de démarcation localisée dans l'étape 106, et la vitesse Vradial_tissu Par ailleurs déterminée au cours de l'étape 104. En référence à la figure 6, la détermination 108 de la vitesse Vradial_os comprend les 25 étapes suivantes. Pour chacune des deuxièmes valeurs candidates mémorisées dans la mémoire 10, le processeur 8 construit une deuxième image montrant du tissu osseux cortical de l'os B et l'endoste E de l'os, à partir des premiers signaux d'écho, de la vitesse Vradial_tissu déterminée, et de la première courbe de démarcation, et sous l'hypothèse que 30 la vitesse Vradial_os est égale à la deuxième valeur candidate (étape 300). La construction 300 d'une deuxième image à partir d'une deuxième valeur candidate comprend la mise en oeuvre de sous-étapes similaires à celles de l'étape de construction 200 d'une première image, à quelques différences près. Comme l'étape 200, l'étape 300 utilise la méthode de migration de Kirchhoff ou la 35 méthode dite - de mise au point totale » (- Total Focusing Method » en anglais). 16 Il est en particulier à nouveau mis en oeuvre une estimation de trajectoires d'ondes passées par un point P, en exploitant à nouveau le principe de Fermat, selon lequel il est supposé qu'une onde se propage rectilignement dans un milieu homogène. Toutefois, le corps C est cette fois considéré au cours de la mise en oeuvre de cette estimation de 5 trajectoires comme un milieu hétérogène : le tissu biologique est considéré comme un premier milieu homogène, dans lequel les premières ondes ultrasonores se sont propagées à la vitesse Vradtal_ttssu Préalablement déterminée. Par ailleurs, l'os B est considéré comme un deuxième milieu dans lequel ces premières ondes se sont propagées à une des deuxièmes valeurs candidates. La courbe de démarcation est ainsi une ligne d'interface 10 entre les deux milieux, en laquelle une onde subit une réfraction. En combinant ces hypothèses avec la courbe de démarcation préalablement localisée, il est ainsi possible d'estimer au cours de la construction 300 d'une deuxième image les trajectoires de premières ondes ultrasonores non seulement passées par un point P se trouvant dans le tissu biologique non osseux, mais encore les trajectoires de premières 15 ondes qui sont passées par un point P se trouvant dans l'os B. On a représenté en figure 7 un exemple de trajectoire d'une onde ultrasonore de compression depuis l'émetteur 6 d'indice i jusqu'au récepteur 6 d'indice j, passant par un point P de l'os (plus précisément au niveau du tissu cortical de l'os), qui a été estimée au cours de la mise en oeuvre de l'étape 300. Cette trajectoire estimée est continue par 20 morceaux et comprend les quatre segments successifs suivants : - un premier segment dans le tissu biologique non osseux T allant jusqu'à un premier point de la deuxième courbe de démarcation, - un deuxième segment prolongeant le premier segment à partir de la deuxième courbe de démarcation jusqu'au point P, le deuxième segment étant non parallèle 25 au premier segment du fait de la réfraction de l'onde lorsqu'elle change de milieu en le premier point de la deuxième courbe de démarcation, - un troisième segment allant du point P vers un autre point de la deuxième courbe de démarcation, et - un quatrième segment partant du deuxième point de la deuxième courbe de 30 démarcation et allant vers le récepteur d'indice j, le quatrième segment étant non parallèle au troisième segment du fait de la réfraction de l'onde lorsqu'elle change à nouveau de milieu en le deuxième point de la deuxième courbe de démarcation. Au cours de l'étape de construction 300 d'une deuxième image, le processeur 8 calcule ensuite des durées de propagation des ondes ultrasonores via les trajectoires 35 estimées et des intensités de pixels à partir des durées de propagation, des premiers 17 signaux d'écho et de positions des émetteurs et des récepteurs, de la même manière que celle décrite dans le cadre de la construction d'une des premières images. Ici, la durée de propagation tT(i,P) depuis l'émetteur d'indice i jusqu'au point P couvre le premier et le deuxième segment, et la durée de propagation tR(j; P), j,] depuis le 5 point P jusqu'au récepteur d'indice j couvre le troisième et le quatrième segment. En définitive, la construction 300 d'une deuxième image utilise globalement les mêmes principes que la construction 200 d'une première image, mais utilise davantage de données d'entrées, notamment la courbe de démarcation, permettant ainsi aux deuxièmes images de montrer des informations exploitables au niveau de l'os B. 10 On a représenté en figure 5b un exemple de deuxième image que l'on peut obtenir par la mise en oeuvre de l'étape de construction 300. Comme on peut le constater, cette deuxième image est relativement riche en information non seulement au niveau du tissu biologique non osseux T entourant l'os B, mais également à l'intérieur de l'os B, au niveau du périoste PE et au niveau de l'endoste E. 15 Le processeur 8 calcule ensuite, pour chaque deuxième image, au moins une deuxième métrique indicative d'une qualité de mise au point du tissu osseux cortical de l'os B et/ou de l'endoste E de l'os B dans la deuxième image (étape 302). Les troisièmes métriques sont par exemple du même type que les premières métriques. 20 Le processeur 8 sélectionne ensuite en tant que valeur définitive pour la vitesse Vradial_os une valeur optimale parmi les deuxièmes valeurs candidates utilisées pour produire les deuxièmes images (étape 304). Le processeur 8 se base pour cela sur les deuxièmes métriques, de la même manière qu'au cours de l'étape 204. On a représenté en figure 8 différentes courbes de deuxième métrique en fonction de 25 la deuxième valeur candidate utilisée. On constate que, quelle que soit la méthode utilisée pour calculer la deuxième métrique, la deuxième valeur sélectionnée en tant que vitesse Vradial_os est sensiblement la même (les maximas de ces courbes sont en effet très proches les uns des autres). On a ainsi déterminé la vitesse Vradial_os réelle dans l'os B du corps C étudié. L'étape 30 108 est alors terminée. Cette vitesse Vradial_os constitue en soi une donnée intéressante car peut servir de donnée d'entrée pour diagnostiquer une ostéoporose éventuelle de l'os B du corps C étudié. Il est toutefois avantageux de compléter cette information par d'autres vitesses afin 35 de caractériser de manière plus complète l'os B étudié. A cette fin, le procédé mis en 18 oeuvre par le système 1 de caractérisation d'os B comprend les étapes additionnelles suivantes, en référence à la figure 9. La sonde 2 est positionnée à proximité du corps C dans une deuxième position différente de la première position. Dans la deuxième position, les émetteurs-récepteurs 6 5 de la sonde 2 sont alignés dans une deuxième direction différente de la première direction. La sonde 2 émet des deuxièmes ondes ultrasonores en direction du corps C dans la deuxième direction (étape 400). Les deuxièmes ondes sont de préférence des ondes de compression axiales. Autrement dit, la deuxième direction est par exemple comprise dans un plan parallèle à l'axe 10 longitudinal de l'os B. Cette deuxième direction est idéalement parallèle à l'axe longitudinal de l'os, mais peut en pratique former un angle a de petite valeur avec l'axe longitudinal de l'os B. Les deuxièmes ondes ultrasonores pénètrent dans le corps C et sont réfléchies par ce derniers à différents niveaux de profondeur dans le corps C. 15 Des nouveaux signaux d'écho de ces deuxièmes ondes, dits deuxièmes signaux d'écho, sont ainsi acquis par les émetteurs-récepteurs 6 de la sonde 2 (étape 402). Les signaux d'écho sont numérisés, transmis au dispositif de traitement 4 et mémorisés dans la mémoire 10. Le processeur 8 détermine ensuite sur la base de deuxièmes signaux d'écho, la vitesse 20 Vaxtal_os (étape 404). L'étape 404 de détermination de la vitesse VaxtaLos peut être mise en oeuvre au moyen d'une méthode connue de l'état de la technique se fondant sur un repérage d'une onde de tête qui se propage le long de la surface extérieure de l'os B. Selon cette technique connue, on utilise deux deuxièmes signaux d'écho spécifiques 25 acquis au moyen de deux émetteurs-récepteurs 6 extrémaux de la sonde 2 (typiquement, celui d'indice 0 et celui d'indice maximal). L'une des deux émetteurs-récepteurs 6 reçoit les ondes émises sélectivement par l'autre émetteur-récepteur 6 et vice versa. La distance qui sépare les deux émetteurs-récepteurs 6 est connue, typiquement 28 cm. A cette échelle, il peut être fait l'hypothèse que la surface extérieure de l'os B est plane. Dès lors, 30 la vitesse de propagation d'une onde de tête le long de cette surface peut être déterminée facilement étant donné que la relation entre instant d'arrivée d'une onde par l'une des deux émetteurs-récepteurs 6 utilisés et la distance séparant les deux émetteurs-récepteurs 6 est une fonction linéaire. Sous cette hypothèse de linéarité, il est très facile de déterminer : 19 - la vitesse V1 d'une onde de tête qui s'est propagée le long de l'os B lorsque l'émetteur-récepteur 6 d'indice 0 était utilisé en tant qu'émetteur, et que l'émetteur-récepteur 6 d'indice maximal était utilisé en tant que récepteur, et - la vitesse V2 d'une onde de tête qui s'est propagée le long de l'os B lorsque 5 l'émetteur-récepteur 6 d'indice 0 était utilisé en tant qu'émetteur, et que l'émetteur-récepteur 6 d'indice maximal était utilisé en tant que récepteur, et La vitesse V,,,,,,Los est ensuite calculée par le processeur 8 au moyen de la formule suivante : 2V1V2 cos (a) Vaxial_os + V2 10 Par ailleurs, le processeur 8 détermine la vitesse Vaxial_tissu sur la base des deuxièmes signaux d'écho (étape 404). Cette étape utilise les mêmes principes fondamentaux que l'étape 104 ayant permis de déterminer la vitesse Vradial_tissu, à la différence près que l'étape 404 utilise les deuxièmes ondes qui ont été émises alors que la sonde 2 occupait la deuxième position et non la première position. 15 La détermination 404 de la vitesse Vaxial_tissu comprend plus précisément les étapes suivantes, en référence à la figure 10: - pour plusieurs troisièmes valeurs candidates, construction d'une troisième image montrant le tissu biologique et le périoste PE de l'os, à partir des deuxièmes signaux d'écho et sous l'hypothèse que la vitesse Vaxial_tissu est égale à la deuxième valeur 20 candidate (étape 500), - pour chaque troisième image, calcul d'une troisième métrique indicative d'une qualité de mise au point du périoste PE dans la troisième image (étape 502), - sélection d'une des troisièmes valeurs candidates en tant que vitesse Vaxial_tissu, en fonction des troisièmes métriques (étape 504). 25 Les troisièmes images sont des vues du corps C selon une coupe longitudinale parallèle à l'axe X de l'os B. Les troisièmes métriques sont par exemple du même type que les premières métriques et /ou des deuxièmes métriques. De retour à la figure 9, le processeur 8 localise une deuxième courbe de démarcation 30 entre le tissu biologique non osseux et l'os B dans l'une des troisièmes images construites (étape 408). Très préférentiellement, la localisation 408 de la deuxième courbe est mise en oeuvre dans la troisième image construite à l'aide de la troisième valeur candidate sélectionnée en tant que vitesse Vaxial_tissu- La localisation est alors beaucoup plus précise du fait de la 20 haute qualité de mise au point de cette troisième image parmi toutes celles qui ont été construites par le processeur 8 au cours de l'étape 500. Le processeur 8 détermine par ailleurs le paramètre d'anisotropie de l'os B étudié (étape 412). 5 Comme indiqué précédemment, ce paramètre est susceptible d'être utilisé par la fonction prédéterminée V en combinaison avec la vitesse Vrodtai_ os et avec la vitesse Vaxtal_ os pour calculer une vitesse de propagation V(0) d'une onde de compression dans l'os B suivant une direction quelconque formant un angle 6I avec l'axe X. Le paramètre d'anisotropie de l'os B est déterminé 412 au moyen des sous-étapes 10 suivantes, en référence à la figure 11. Pour chaque quatrième valeur candidate mémorisée dans la mémoire 10, le processeur 8 construit une quatrième image montrant du tissu osseux cortical et l'endoste E de l'os, à partir des deuxièmes signaux d'écho, de la vitesse Vrodtal os, de la vitesse Vaxtal os, de la fonction prédéterminée V, et sous l'hypothèse que le paramètre 15 d'anisotropie de l'os B est égal à la quatrième valeur candidate (étape 600). Cette étape 600 est similaire aux étapes 200, 300 et 500 en ce sens qu'elle utilise une des méthodes de migration de Kirchhoff ou dite de mise au point totale » (- Total Focusing Method » en anglais). Il convient de noter simplement que les trajectoires d'ondes estimées au cours de l'étape 600 se fondent sur l'hypothèse que la vitesse de propagation 20 d'une onde dans l'os B observé dans une direction quelconque respecte le modèle défini par la fonction V. On obtient donc à l'issue de l'étape 600 une pluralité de quatrièmes images, une pour chaque quatrième valeur candidate mémorisée initialement dans la mémoire 10. Pour chaque quatrième image, le processeur 8 calcule au moins une quatrième 25 métrique indicative d'une qualité de mise au point de l'endoste E et/ou du tissu biologique T entourant l'os B dans la quatrième image (étape 602). Les quatrièmes métriques sont par exemple du même type que les premières métriques. Le processeur 8 sélectionne ensuite une des quatrièmes valeurs candidates en tant que paramètre d'anisotropie de l'os B, fi', en fonction des quatrièmes métriques (étape 30 604). On a représenté en figure 12 différentes courbes de quatrièmes métriques en fonction de la quatrième valeur candidate utilisée. On constate que, quelle que soit la méthode utilisée pour calculer la quatrième métrique, la quatrième valeur sélectionnée en tant que paramètre est sensiblement la même (les maximas de ces courbes sont en effet 21 très proches les uns des autres), comme cela était également le cas pour la vitesse Vradial_os - Il est désormais possible pour le processeur 8 de calculer n'importe quelle vitesse de propagation d'une onde de compression dans l'os B en utilisation la fonction V paramétrée 5 avec les données Vradial_os Vaxial_os et /3 déterminées au cours du procédé. Dans le mode de réalisation de procédé qui vient d'être décrit, les différentes données caractéristiques de l'os B que sont Vradial_os, Vaxial os, et ont été déterminées dans des étapes distinctes, et non de manière conjointe. En particulier, les étapes de détermination 104, 108 et 412 n'ont utilisé chacune qu'un seul ensemble de valeurs 10 candidates (deuxièmes valeurs candidates pour l'étape 104, troisièmes valeurs candidates pour l'étape 108, quatrièmes valeurs candidates pour l'étape 412). Ceci permet en définitive de limiter le nombre de deuxièmes, troisièmes et quatrièmes images construites. Autrement dit les ressources matérielles (charge de calcul, mémoire 10) requises pour la mise en oeuvre de ces étapes sont relativement modérées. 15 Avantageusement, le processeur 8 commande un affichage l'écran d'affichage 12 au moins une des images suivantes : - la première image construite à l'aide de la première valeur sélectionnée au cours de l'étape 204, - la deuxième image construite à l'aide de la deuxième valeur sélectionnée au cours 20 de l'étape 304, - la troisième image construite à l'aide de la première valeur sélectionnée au cours de l'étape 504, - la deuxième image construite à l'aide de la deuxième valeur sélectionnée au cours de l'étape 604. 25 Ces images présentent en effet une qualité de mise au point élevée et sont par conséquent d'intérêt pour un praticien. L'invention ne se limite pas au seul mode de réalisation qui vient d'être décrit en relation avec les figures. L'invention est notamment applicable à l'étude d'ondes de cisaillement. 30 Dans ce cas, il est possible de modéliser la vitesse de propagation d'une onde de cisaillement dans l'os B au moyen de la fonction prédéterminée V de la forme suivante : V(0) = V15 axial os [1 + 10' sin2(0) cos2(0)] où - Vcis axial_os est une vitesse de propagation d'une onde de cisaillement dans l'os B 35 suivant une première direction parallèle à l'axe longitudinal X, 22 - B est un angle quelconque, - V(0) est une vitesse de propagation d'une onde de cisaillement dans l'os B suivant une direction formant l'angle B avec la première direction. On a par exemple représenté en figure 13 un organigramme d'étapes d'un procédé de 5 caractérisation de l'os selon un deuxième mode de réalisation. La sonde 2 est positionnée dans la deuxième position précédemment décrite, et des ondes ultrasonores de compression sont émises par la sonde 2 dans cette deuxième position (étape 600). De façon connue en elle-même, les ondes de compression émises entraînent 10 l'apparition d'ondes de cisaillement dans l'os B. Ces ondes de cisaillement sont générées lorsque les ondes de compression pénètrent dans le tissu osseux cortical de l'os B, situé entre le périoste PE et l'endoste E. Des signaux d'écho de ces ondes de cisaillement sont ensuite reçus par la sonde 2 (étape 602). 15 Le processeur 8 détermine ensuite une vitesse Vaxtal ttssu de propagation axiale d'ondes de compression dans le tissu biologique non osseux T (étape 704). Cette étape 604 est similaire à l'étape 404. Le processeur 8 localise une courbe de démarcation entre le tissu T et l'os B (étape 706). Cette étape est identique à l'étape 406. 20 Le processeur 8 détermine ensuite la vitesse Vcts oxtai os (étape 708) dans l'os B, et, conjointement, le paramètre d'anisotropie fi'. Cette détermination conjointe comprend les sous-étapes suivantes, en référence à la figure 14. On suppose qu'a été préalablement mémorisé dans la mémoire un jeu de cinquièmes valeurs candidates, pour la vitesse Vcts oxtai 0 s . Ces cinquièmes valeurs sont comprises entre 25 1300 et 2000 mètres par seconde. Par ailleurs, les quatrièmes valeurs candidates sont, dans ce mode de réalisation comprises entre -0,3 et 0.3. Pour une paire constituée d'une quatrième valeur candidate (pour fi') et d'une cinquième valeur candidate, le processeur construit une image à partir des signaux d'écho, 30 en faisant la double hypothèse que la vitesse Vcts oxtoLos est égale à la cinquième valeur candidate de la paire, et que le paramètre d'anisotropie est égal à la quatrième valeur de la paire (étape 800). Le processeur 8 répète cette étape 800 pour chaque paire disponible. 23 On constate ici que le nombre d'images construites est bien plus grand que le nombre d'images générées au cours d'une étape du procédé selon le premier mode de réalisation, car l'on travaille sur deux dimensions au lieu d'une seule. Pour chaque image construite, le processeur calcule une métrique de qualité de mise 5 au point (étape 802). Cette étape est similaire aux étapes 502, 602. Le processeur 8 sélectionne ensuite une des paires utilisées pour construire les images au cours de l'étape 800 (étape 804), sur la base des métriques calculées. Cette étape est similaire aux étapes 504, REVENDICATIONS1. Procédé de caractérisation d'un os (B), le procédé comprenant des étapes de : - émission (100, 700) de premières ondes ultrasonores vers un corps comprenant l'os (B) et du tissu biologique non osseux (T) entourant l'os (B), - réception (102, 702) de premiers signaux d'écho des premières ondes ultrasonores émises, - détermination (104, 704) d'une vitesse du son (V radial _tissu) dans le tissu biologique non osseux (T) suivant une première direction, la détermination (104, 704) comprenant des sous-étapes de : o pour plusieurs premières valeurs candidates prédéterminées, construction (200) d'une première image montrant le tissu biologique non osseux (T) et le périoste (PE) de l'os (B), à partir des premiers signaux d'écho et sous l'hypothèse que la vitesse (V radial_tissu) du son dans le tissu biologique non osseux (T) suivant la première direction est égale à la première valeur candidate, o pour chaque première image, calcul (202) d'au moins une première métrique indicative d'une qualité de mise au point du périoste et/ou du tissu biologique non osseux (T) entourant l'os dans la première image, o sélection (204) d'une des premières valeurs candidates en tant que vitesse du son (Vradial_tissu) dans le tissu biologique non osseux (T) suivant la première direction, en fonction des premières métriques, - localisation (106, 706) d'une première courbe de démarcation entre le tissu biologique non osseux (T) et le périoste dans l'une des premières images, - détermination (108, 708) d'une vitesse du son dans l'os (Vradial os) suivant la première direction, la détermination (108, 708) comprenant des sous-étapes de : o pour plusieurs deuxièmes valeurs candidates prédéterminées, construction (300, 800) d'une deuxième image montrant du tissu osseux cortical de l'os (B) et l'endoste (E) de l'os (B), à partir des premiers signaux d'écho, de la vitesse du son (V radial_tissu) dans le tissu biologique non osseux (T) suivant la première direction déterminée, et de la première courbe de démarcation, et sous l'hypothèse que la vitesse du son dans l'os (Vradlotos) suivant la première direction est égale à la deuxième valeur candidate, 25 o pour chaque deuxième image, calcul (302, 802) d'au moins une deuxième métrique indicative d'une qualité de mise au point du tissu osseux cortical de l'os et/ou de l'endoste de l'os dans la deuxième image, o sélection (304, 804) d'une des deuxièmes valeurs candidates en tant que 5 vitesse du son dans l'os (1/'.ad,,,/ ') suivant la première direction, en fonction des deuxièmes métriques. 2. Procédé selon la revendication 1, dans lequel la construction (200) d'une première image à l'aide d'une première valeur candidate comprend les étapes suivantes mises en 10 oeuvre pour au moins un point du tissu biologique non osseux (T) : - estimation de premières trajectoires de premières ondes ultrasonores émises par des émetteurs, puis passées par le point du tissu biologique non osseux (T), puis reçues par des récepteurs, à partir des premiers signaux d'écho, sous l'hypothèse que la vitesse du son dans le tissu biologique non osseux (T) suivant la première 15 direction est égale à la première valeur candidate, - calcul de durées de propagation des ondes ultrasonores via les premières trajectoires estimées, - calcul d'une intensité d'un pixel de la première image au point du tissu biologique non osseux (T), à partir des durées de propagation, des premiers signaux d'écho et 20 de positions des émetteurs et des récepteurs. 3. Procédé selon l'une des revendications précédentes, dans lequel la construction (300, 800) d'une deuxième image à l'aide d'une deuxième valeur candidate comprend les étapes suivantes mises en oeuvre pour au moins un point de l'os (B) : 25 - estimation de deuxièmes trajectoires de premières ondes ultrasonores émises par des émetteurs, puis passées par le point de l'os, puis reçues par des récepteurs, à partir des signaux d'écho, à partir de la vitesse du son dans le tissu biologique non osseux (T) (Vradial tissu) suivant la première direction déterminée, de la première courbe de démarcation, et sous l'hypothèse que la vitesse du son (Vradial_os) dans 30 l'os suivant la première direction est égale à la deuxième valeur candidate, - calcul de durées de propagation des ondes ultrasonores via les deuxièmes trajectoires estimées, - calcul d'une intensité d'un pixel de la deuxième image au point de l'os (B), à partir des durées de propagation, des premiers signaux d'écho et de positions des émetteurs et des récepteurs. 26 4. Procédé selon l'une des revendications précédentes, dans lequel la localisation de la première courbe de démarcation est mise en oeuvre dans la première image qui a été construite à l'aide de la première valeur sélectionnée (204) en tant que vitesse du son 5 (V radial_tissu) dans le tissu biologique non osseux (T) suivant la première direction. 5. Procédé selon l'une des revendications précédentes, dans lequel les premières ondes ultrasonores sont des ondes émises par des émetteurs et les signaux d'échos reçus par des récepteurs alignés le long d'un axe perpendiculaire à un axe longitudinal (X) de l'os (B), et 10 dans lequel la première direction est une direction perpendiculaire à l'axe longitudinal (X) de l'os (B). 6. Procédé selon l'une des revendications précédentes, comprenant en outre un affichage de la première image construite à l'aide de la première valeur sélectionnée, et/ou de la 15 deuxième image construite à l'aide de la deuxième valeur sélectionnée. 7. Procédé selon la revendication précédente, comprenant en outre des étapes de : - émission (400) de deuxièmes ondes ultrasonores vers le corps, - réception (402) de signaux d'écho des deuxièmes ondes ultrasonores émises, dits 20 deuxièmes signaux d'écho, - détermination (404) d'une vitesse du son (Vaxial_tissu) dans le tissu biologique non osseux (T) suivant une deuxième direction différente de la première direction, comprenant des sous-étapes de : o pour plusieurs troisièmes valeurs candidates, construction (500) d'une 25 troisième image montrant le tissu biologique non osseux (T) et le périoste de l'os, à partir des deuxièmes signaux d'écho et sous l'hypothèse que la vitesse du son dans le tissu biologique non osseux (T) suivant la deuxième direction est égale à la deuxième valeur candidate, o pour chaque troisième image, calcul (502) d'au moins une troisième 30 métrique indicative d'une qualité de mise au point du périoste dans la troisième image, o sélection (504) d'une des troisièmes valeurs candidates en tant que vitesse du son (Vaxial_tissu) dans le tissu biologique non osseux (T) suivant la deuxième direction, en fonction des troisièmes métriques, 27 - détermination (406) d'une vitesse du son (V,,,,,,L') dans l'os suivant la deuxième direction, à l'aide des deuxièmes signaux d'écho, - localisation (408) d'une deuxième courbe de démarcation entre le tissu biologique non osseux (T) et le périoste dans l'une des troisièmes images, 5 - détermination (410) d'un paramètre d'anisotropie de l'os (/g) susceptible d'être utilisé par une fonction prédéterminée (V) en combinaison avec la vitesse du son dans l'os suivant la première direction et avec la vitesse du son dans l'os suivant la deuxième direction pour calculer une vitesse du son dans l'os suivant une direction quelconque, la détermination du paramètre d'anisotropie (/g) comprenant des sous- 10 étapes de - o pour plusieurs quatrièmes valeurs candidates prédéterminées, construction (600) d'une quatrième image montrant du tissu osseux cortical et l'endoste de l'os, à partir des deuxièmes signaux d'écho, de la vitesse de son dans l'os suivant la première direction (Vradtal_os), de la deuxième courbe de 15 démarcation, de la vitesse de son dans l'os suivant la deuxième direction (Vaxtal_os), de la fonction prédéterminée (V), et sous l'hypothèse que le paramètre d'anisotropie de l'os (/g) est égal à la quatrième valeur candidate, o pour chaque quatrième image, calcul (602) d'une quatrième métrique indicative d'une qualité de mise au point de l'endoste et/ou du tissu osseux 20 cortical de l'os dans la quatrième image, o sélection (604) d'une des quatrièmes valeurs candidates en tant que paramètre d'anisotropie de l'os (/g), en fonction des quatrièmes métriques. 8. Procédé selon la revendication 7, dans lequel la construction d'une quatrième image 25 comprend les étapes suivantes mises en oeuvre pour au moins un point de l'os : - estimation de troisièmes trajectoires de deuxièmes ondes ultrasonores émises par des émetteurs, puis passées par le point de l'os, puis reçues par des récepteurs, à partir des deuxièmes signaux d'écho, de la vitesse de son dans l'os (V radial os) suivant la première direction déterminée, de la vitesse du son dans l'os (V axial os) 30 suivant la deuxième direction déterminée, de la deuxième courbe de démarcation, de la fonction prédéterminée (V) et sous l'hypothèse que le paramètre d'anisotropie de l'os ()g) est égale à la quatrième valeur candidate, - calcul de durées de propagation des ondes ultrasonores via les troisièmes trajectoires estimées, 28 - calcul d'une intensité d'un pixel de la quatrième image au point de l'os, à partir des durées de propagation et des deuxièmes signaux d'écho et de positions des émetteurs et des récepteurs. 5 9. Procédé selon l'une des revendication 7 à 8, dans lequel la localisation (408) de la deuxième courbe de démarcation est mise en oeuvre dans la troisième image ayant été construite à l'aide de la troisième valeur sélectionnée en tant que vitesse du son dans le tissu biologique non osseux (T) suivant la deuxième direction. 10 10. Procédé selon l'une des revendication 7 à 9, comprenant en outre un affichage de la troisième image construite à l'aide de la troisième valeur sélectionnée, et/ou de la quatrième image construite à l'aide de la quatrième valeur sélectionnée et de la vitesse du son dans l'os suivant la deuxième direction (V,,,,,a '). 15 11. Procédé selon l'une des revendications 7 à 10, dans lequel les deuxièmes ondes ultrasonores sont des ondes émises par des émetteurs et les deuxièmes signaux d'écho sont reçus par des récepteurs alignés le long d'un axe compris dans un plan parallèle à un axe longitudinal (X) de l'os (B), et dans lequel la deuxième direction est une direction parallèle à l'axe longitudinal (X) de l'os (B). 20 12. Procédé selon l'une des revendications précédentes, dans lequel chaque vitesse de son dans l'os (Vradial os, Vaxial os) déterminée est une vitesse de propagation d'ondes de compression ou une vitesse de propagation d'ondes de cisaillement. 25 13. Procédé selon l'une des revendications 7 à 12, dans lequel chaque vitesse de son dans l'os (V radial os, Vaxial os) déterminée est une vitesse de propagation d'ondes de compression, et dans lequel la fonction prédéterminée est une fonction V de la forme : V(e) = Vaxial_os (Vaxial_os Vradial_os)W sin2(0) cos2(0) + cos4(0)] dans lequel de compression dans l'os suivant la de compression dans l'os suivant la 30 - v - radial_os est la vitesse de propagation d'ondes première direction déterminée, - V axial_os est la vitesse de propagation d'ondes deuxième direction déterminée, - )0 est le paramètre d'anisotropie de l'os, 35 - 0 est un angle, 29 - V(0) est une vitesse de propagation d'ondes de compression dans l'os suivant une direction formant l'angle 0 avec la première direction. 5 14. Procédé selon l'une des revendications précédentes, dans lequel la vitesse de son dans l'os suivant la première direction déterminée est une vitesse de propagation d'ondes de cisaillement, et dans lequel la fonction prédéterminée est une fonction V de la forme : V(0) = Vcis_axial_os [1 + sin2(0) cos2(0)] dans lequel 10 - V cis_axial_os est la vitesse de propagation d'ondes de cisaillement dans l'os suivant la première direction déterminée, - )0 est le paramètre d'anisotropie de l'os, - 0 est un angle, - V(0) est une vitesse de propagation d'ondes de cisaillement dans l'os suivant une 15 direction formant l'angle 0 avec la première direction. 15. Procédé selon l'une des revendications 1 à 7, dans lequel il est déterminé (708), conjointement à la vitesse (Vcts_axtaLos) du son dans l'os suivant la première direction, un paramètre d'anisotropie de l'os (/0) susceptible d'être utilisé par une fonction 20 prédéterminée en combinaison avec la vitesse du son dans l'os suivant la première direction déterminée pour calculer une vitesse du son dans l'os suivant une direction quelconque, ladite détermination conjointe comprenant des sous-étapes de : - pour plusieurs paires de deuxièmes et troisièmes valeurs candidates prédéterminées, construction (800) d'une deuxième image montrant du tissu osseux 25 cortical de l'os et l'endoste de l'os, à partir des premiers signaux d'écho, de la vitesse du son dans le tissu biologique suivant la première direction déterminée, et de la première courbe de démarcation (périoste), sous l'hypothèse que la vitesse du son dans l'os suivant la première direction est égale à la deuxième valeur candidate, et sous l'hypothèse que le paramètre d'anisotropie de l'os ()0) est égal à 30 la troisième valeur candidate, - pour chaque deuxième image, calcul (802) d'au moins une deuxième métrique indicative d'une qualité de mise au point du tissu osseux cortical de l'os et/ou de l'endoste de l'os dans la deuxième image, - sélection (804) conjointe d'une deuxième valeur candidate en tant que vitesse du 35 son dans l'os suivant la première direction, et d'une troisième valeur candidate en 30 tant que paramètre d'anisotropie de l'os (/g), la sélection étant mise en oeuvre en fonction des deuxièmes métriques.