FR 3004635 A1 20141024 FR 1353696 A 20130423 - 1 - « Procédé de correction/calibration d'un paramètre de caractérisation d'un tissu osseux » Domaine technique La présente invention concerne un procédé de correction/calibration d'un paramètre de caractérisation d'un tissu osseux. Elle concerne également un système mettant en oeuvre un tel procédé. Le domaine de l'invention est le domaine médical et plus précisément le domaine de l'étude des tissus osseux avec des appareils d'imagerie comprenant un émetteur de rayons de particules électriquement neutres et mettant en oeuvre un capteur de rayons. L'invention s'applique plus particulièrement au domaine de l'imagerie médicale de haute précision et de l'étude de la caractérisation des os. Etat de ma technique Actuellement, il existe de nombreux appareils d'imagerie de tissus osseux. Ces appareils comprennent une source de particules électriquement neutres prévue pour émettre un rayonnement et éventuellement un récepteur doté d'un capteur numérique ou analogique. Le tissu osseux est positionné entre l'émetteur et le récepteur. Le procédé mis en oeuvre dans ces appareils comprend une émission d'un rayonnement de particules électriquement neutres qui viennent frapper le tissu osseux. Les particules électriquement neutres traversent plus ou moins le milieu étudié en fonction de sa densité. L'atténuation est plus importante dans un milieu osseux que dans les tissus mous. L'image obtenue par ces procédés est une image en niveaux de gris et le niveau de gris en un point de l'image obtenue dépend de la quantité de particules frappant le capteur en ce point. L'image en niveaux de gris est ensuite analysée et traitée pour déterminer la valeur d'un ou plusieurs paramètres relatifs au milieu osseux, aussi appelé « paramètre de caractérisation » dans la présente demande. Un tel paramètre de - 2 - caractérisation peut par exemple être le paramètre de texture H, aussi connu sous le nom de paramètre de Hurst. La valeur mesurée d'un paramètre de caractérisation d'un tissu osseux dépend fortement des conditions de mesure, et plus particulièrement du capteur utilisé pour générer l'image en niveaux de gris et de la dose de particules émise pour l'acquisition de l'image en niveaux de gris. Par conséquent il est difficile, voire impossible, de comparer, pour un même paramètre de caractérisation, deux valeurs mesurées dans deux configurations de mesures différentes. Actuellement, pour réaliser une telle comparaison, l'une des valeurs mesurées est calibrée/corrigée en utilisant une régression linéaire de type Mesure (configuration 1)= constante* Mesure (Configuration 2) s'appliquant directement aux valeurs mesurées du paramètre de caractérisation. Cependant, une telle calibration/correction donne des valeurs très approximatives, voire aberrantes, et ne peut donc être utilisée. Un but de l'invention est de pallier les inconvénients cités ci-dessus. Un autre but de l'invention est de proposer un procédé et système de calibration/correction réalisant une calibration/correction de la valeur mesurée d'un paramètre de caractérisation d'un tissu osseux en fonction de la configuration de mesure utilisée, fournissant de meilleurs résultats quelle que soit la configuration de mesure. Un autre but de la présente invention est de proposer un procédé et système de calibration/correction réalisant une calibration/correction de la valeur mesurée d'un paramètre de caractérisation d'un tissu osseux en fonction de la configuration de mesure utilisée, fournissant de meilleurs résultats quel(le) que soit le capteur utilisé ou la dose de particules utilisées. Enfin, un autre but de l'invention est de proposer un procédé et 30 système de mesure d'un paramètre de caractérisation d'un tissu osseux mettant en oeuvre une telle calibration/correction. - 3 - Exposé de l'invention L'invention propose d'atteindre au moins l'un de ces buts par un procédé de correction/calibration d'une valeur, dite mesurée, d'un paramètre de caractérisation d'un tissu osseux déterminée par mesure fractale sur une image en niveaux de gris réalisée dans une première configuration, dite de mesure, par rapport à une deuxième configuration, dite de référence, se distinguant de ladite configuration de mesure par une différence affectant la valeur mesurée dudit paramètre de caractérisation, ledit procédé étant caractérisé en ce qu'il comprend les étapes suivantes : - détermination, d'une dimension fractale, dite mesurée, de ladite image en niveaux de gris réalisée dans la configuration de mesure, - correction de ladite dimension fractale mesurée en référence à ladite configuration de référence, ladite correction fournissant une valeur dite corrigée de ladite dimension fractale, - détermination d'une valeur, dite corrigée, dudit paramètre de caractérisation à partir de ladite valeur corrigée de ladite dimension fractale. Le procédé selon l'invention propose donc de réaliser une correction d'une valeur mesurée d'un paramètre de caractérisation d'un tissu osseux à partir de la dimension fractale de l'image en niveau de gris utilisée pour déterminer ladite valeur mesurée. Le procédé selon l'invention permet donc de réaliser la correction/calibration du paramètre de caractérisation de manière indirecte, c'est-à-dire sans utiliser une relation qui s'applique directement sur la valeur mesurée du paramètre de caractérisation, alors que les procédés de l'état de la technique réalisent une correction se basant directement sur la valeur mesurée du paramètre de caractérisation. L'inventeur de la présente invention a découvert que la dimension fractale de l'image en niveaux de gris constitue une base permettant de réaliser une meilleure correction/calibration prenant mieux en compte les différences dans les configurations de mesure comparée à une correction/calibration basée sur la valeur mesurée du paramètre de caractérisation. - 4 - La correction/calibration selon l'invention permet de fournir des résultats répétables et confirmés par des essais cliniques, permettant ainsi une comparaison correcte de valeurs mesurées dans différentes configurations de mesure. Selon tous les aspects de l'invention, la dimension fractale d'une image en niveaux de gris peut être déterminée : - par analyse de l'image en niveaux de gris, et/ou - suivant une relation prédéterminée prenant en compte la valeur mesurée du paramètre de caractérisation. Selon l'invention, l'étape de correction de la dimension fractale mesurée, par rapport à la configuration de référence, peut être réalisée suivant une relation prédéterminée de manière empirique avec plusieurs échantillons de mesure pour lesquels les valeurs du paramètre de caractérisation, et/ou les valeurs de la dimension fractale, sont connues dans la configuration de mesure et dans la configuration de référence. Une telle relation prédéterminée peut être une régression linéaire prenant en compte la dimension fractale de l'image en niveaux de gris obtenue dans la configuration de mesure et la dimension fractale de l'image en niveaux de gris obtenue dans la configuration de référence. Un exemple de relation prédéterminée peut être la suivante : Dm= A Dr avec Dm la dimension fractale de l'image en niveaux de gris obtenue dans la configuration de mesure, D'f la dimension fractale de l'image en niveaux de gris obtenue dans la configuration de référence, et A une constante déterminée de manière empirique, et qui correspond à la pente de la régression linéaire passant par l'origine réalisée entre Dm et Dr. Selon un premier exemple de réalisation, la configuration de mesure peut différer de la configuration de référence par le détecteur ou capteur utilisé pour réaliser l'image en niveau de gris. Chaque détecteur/capteur peut être numérique ou analogique. - 5 - Dans le cas où le détecteur est un détecteur analogique, le procédé selon l'invention peut en outre comprendre une numérisation de l'image en niveaux de gris. Selon un autre exemple de réalisation, la configuration de mesure peut différer de la configuration de référence par la quantité et/ l'énergie des particules utilisées pour réaliser l'image en niveau de gris. Bien entendu, la configuration de mesure peut différer de la configuration de référence par au moins deux paramètres à la fois, par exemple à la fois par le détecteur utilisé et par la dose de particules utilisées. Dans ce cas, la correction/calibration réalisée prend en compte la combinaison de chacune de ces différences. Selon un autre aspect de l'invention, il est proposé un procédé de détermination/mesure d'un paramètre de caractérisation d'un tissu osseux, ledit procédé comprenant les étapes suivantes : - acquisition d'une image en niveaux de gris dudit tissu osseux, et - détermination d'une valeur, dite corrigée/calibrée, dudit paramètre de caractérisation suivant par le procédé selon l'invention. La valeur mesurée du paramètre de caractérisation peut en particulier être réalisée dans une région, dite d'intérêt, dudit tissu osseux préalablement identifiée par l'intermédiaire de repères anatomique. Une telle région d'intérêt peut être définie et identifiée sur une image en niveaux de gris préalablement acquise. Dans ce cas, le procédé de détermination/mesure selon l'invention comprend une étape préalable d'acquisition d'une première image en niveaux de gris, puis une identification de la région d'intérêt par des marqueurs anatomiques puis une acquisition d'une deuxième image en niveaux de gris, centrée sur la région d'intérêt. La détermination de la valeur mesurée du paramètre de - 6 - caractérisation est alors réalisée sur cette deuxième image en niveaux de gris centrée sur la région d'intérêt. Avantageusement, le paramètre de caractérisation du tissu osseux déterminé peut être le paramètre de texture osseux, également appelé le paramètre de Hurst, la densité osseuse, ou tout autre paramètre de caractérisation d'un tissu osseux dont la valeur est déterminée par mesure fractale sur une image en niveaux de gris. Selon encore un autre aspect de l'invention, il est proposé un dispositif de correction/calibration d'une valeur, dite mesurée, d'un paramètre de caractérisation d'un tissu osseux déterminée par mesure fractale sur une image en niveaux de gris réalisée dans une première configuration, dite de mesure, par rapport à une deuxième configuration, dite de référence, se distinguant de ladite configuration de mesure par une différence affectant la valeur mesurée dudit paramètre de caractérisation, ledit dispositif caractérisé en ce qu'il comprend des moyens agencés et configurés pour : - déterminer une dimension fractale de ladite image en niveaux de gris, - corriger ladite dimension fractale en référence à ladite configuration de référence, ladite correction fournissant une valeur dite corrigée de ladite dimension fractale, et - déterminer une valeur, dite corrigée, dudit paramètre de caractérisation à partir de ladite valeur corrigée de ladite dimension fractale. Selon encore un autre aspect de l'invention il est proposé un système 30 de détermination/mesure d'un paramètre de caractérisation d'un tissu osseux, ledit système comprenant : - des moyens d'imagerie pour l'acquisition d'au moins une image en niveaux de gris dudit tissu osseux, et - un dispositif de correction/calibration selon l'invention. -7 - Les moyens d'imagerie peuvent comprendre au moins un détecteur qui peut être un détecteur numérique ou analogique. Dans le cas où, le détecteur est un détecteur analogique, le système 5 selon l'invention peut en outre comprendre des moyens de numérisation de l'image analogique fournie par le détecteur analogique. Le dispositif, respectivement le système, selon l'invention peut être mis en oeuvre, pour la correction/calibration, respectivement la 10 détermination/mesure, d'au moins un paramètre de caractérisation d'un tissu osseux, par exemple la texture osseuse, définie par un paramètre dérivé de la mesure d'une dimension fractale de l'image, par exemple, le paramètre de Hurst. 15 Description des figures et modes de réalisation D'autres avantages et caractéristiques apparaîtront à l'examen de la description détaillée d'un mode de réalisation nullement limitatif, et des dessins annexés sur lesquels 20 la figure 1 est une représentation sous la forme d'un diagramme des étapes d'un exemple de procédé de correction/calibration selon l'invention de la valeur mesurée d'un paramètre de caractérisation d'un tissu osseux ; la figure 2 est une représentation sous la forme d'un diagramme des 25 étapes d'un exemple de procédé de détermination/mesure selon l'invention de la valeur d'un paramètre de caractérisation d'un tissu osseux ; la figure 3 est une représentation schématique d'un dispositif selon l'invention ; 30 - la figure 4 est une représentation schématique d'un système selon l'invention, et - les figures 5 et 6 correspondent chacune à une représentation schématique d'un exemple de résultats obtenus selon l'invention. - 8 - II est bien entendu que les modes de réalisation qui seront décrits dans la suite ne sont nullement limitatifs. On pourra notamment imaginer des variantes de l'invention ne comprenant qu'une sélection de caractéristiques décrites par la suite isolées des autres caractéristiques décrites, si cette sélection de caractéristiques est suffisante pour conférer un avantage technique ou pour différencier l'invention par rapport à l'état de la technique antérieur. Cette sélection comprend au moins une caractéristique de préférence fonctionnelle sans détails structurels, ou avec seulement une partie des détails structurels si c'est cette partie qui est uniquement suffisante pour conférer un avantage technique ou pour différencier l'invention par rapport à l'état de la technique antérieur. La FIGURE 1 est une représentation sous la forme d'un diagramme des étapes d'un exemple non limitatif d'un procédé de correction/calibration selon l'invention de la valeur mesurée d'un paramètre de caractérisation d'un tissu osseux. Dans l'exemple donné en FIGURE 1, la configuration de mesure diffère d'une configuration de référence prédéterminée soit par le 20 capteur/détecteur utilisé, soit par la dose utilisée en émission pour l'acquisition de l'image en niveau de gris, ou encore les deux. Le procédé 100 représenté sur la FIGURE 1 comprend une étape 102 réalisant une mesure d'une valeur, dite mesurée, d'un paramètre de caractérisation, par analyse d'une image en niveaux de gris d'un tissu 25 osseux. Dans l'exemple présent le paramètre de caractérisation est, à titre d'exemple uniquement et de manière nullement limitative, le paramètre de Hurst H et la valeur mesurée est notée I-11. Lors d'une étape 104, une valeur D1, dite mesurée, de la dimension fractale D de l'image en niveaux de gris est déterminée. La dimension 30 fractale D est déterminée soit par analyse de l'image en niveaux de gris, soit à partir de la valeur mesurée H1 du paramètre de Hurst, selon une relation prédéterminée. Dans l'exemple présent, la valeur D1 est déterminée à partir de la valeur H1 selon la relation prédéterminée : D1=2-H1 (1) - 9 - Lors d'une étape 106, la dimension fractale mesurée D1 est corrigée par rapport à une configuration de référence. La valeur mesurée D1 est corrigée suivant une relation linéaire préalablement déterminée de manière empirique et reliant la valeur de la dimension fractale dans la configuration de référence et la configuration de mesure, et en particulier en fonction de la pente de la régression linéaire passant par l'origine réalisée entre les valeurs de la dimension fractale dans la configuration de référence et dans la configuration de mesure. Cette relation linéaire est définie en mesurant, pour des échantillons connus, la valeur de la dimension fractale dans la configuration de référence et la configuration de mesure. Dans l'exemple présent, la valeur corrigée Dc' de la dimension fractale mesurée est obtenue suivant la relation : Dcor= D1*1,17 (2) Lors d'une étape 108 la valeur corrigée/calibrée du paramètre de caractérisation Hcor est déterminée à partie de la valeur corrigée/calibrée Dcor de la dimension fractale, par exemple en utilisant la relation (1) ci-dessus, à savoir : Fic0r=2-Dcor (3) La FIGURE 2 est une représentation sous la forme d'un diagramme des étapes d'un exemple non limitatif d'un procédé de détermination/mesure de la valeur d'un paramètre de caractérisation d'un tissu osseux selon l'invention. Le procédé 200 représenté sur la FIGURE 2 comprend une étape 202 d'acquisition d'une image en niveaux de gris d'un tissu osseux. Cette étape 202, peut être suivie d'une étape optionnelle 204 lors de laquelle une région d'intérêt est définie sur l'image en niveau de gris acquise lors de l'étape 202, par des repères anatomiques. Cette étape optionnelle 204 est alors suivie d'une étape 206 d'acquisition d'une nouvelle image en niveaux de gris centrée autour/sur la - 10 - région d'intérêt. Cette étape 206 n'est réalisée que si l'étape optionnelle 204 de définition d'une région d'intérêt est réalisée. L'étape 202, le cas échéant l'étape 206, est suivie des étapes du procédé 100 représenté en figure 1, déterminant une valeur corrigée/calibrée d'une valeur d'un paramètre de caractérisation d'un tissu osseux. La FIGURE 3 est une représentation schématique d'un exemple de dispositif selon l'invention. Le dispositif 300, représenté sur le FIGURE 3, comprend : - un module 302 de mesure d'une valeur, dite mesurée, d'un paramètre de caractérisation, à partir d'une image en niveaux de gris, - un module 304 de mesure d'une valeur, dite mesurée, de la dimension fractale de l'image en niveau de gris à partir soit de l'image en niveaux de gris soit à partir de la valeur mesurée du paramètre de caractérisation, - un module 306 de détermination par calcul d'une valeur corrigée de la dimension fractale par rapport à une configuration de référence, par exemple suivant une relation prédéterminée de manière empirique prenant en compte la valeur mesurée de la dimension fractale, et - un module 308 de détermination par calcul d'une valeur corrigée du paramètre de caractérisation prenant en compte la valeur corrigée de la dimension fractale. Chaque module peut être un moyen de calcul, tel qu'un processeur, agencé pour exécuter un programme informatique ou des lignes de commandes dédiées à l'étape ou l'opération à réaliser. Bien que représentés séparément, au moins deux des modules, en particulier tous les modules, peuvent être intégrés dans un même processeur. - 11 - La FIGURE 4 est une représentation schématique d'un exemple de système selon l'invention. Le système 400 représenté sur la FIGURE 4 comprend un émetteur 402 pour émettre une dose de particules neutres et un récepteur 404, placé en regard de l'émetteur 402, pour l'acquisition d'une image en niveaux de gris. Le tissu osseux à caractériser est disposé entre l'émetteur 402 et le récepteur 404, de sorte que les particules émises par l'émetteur 402 traversent le tissu osseux avant d'arriver sur le récepteur 404. Le système 400 comprend en outre un dispositif 300 tel que représenté sur la FIGURE 3 relié au récepteur 404. Les FIGURES 5 et 6 sont chacune une représentation schématique d'un exemple de résultats obtenus selon l'invention. Les résultats représentés sur les FIGURES 5 et 6 correspondent à la mesure du paramètre de texture H, pour plusieurs échantillons connus appelés PHANTOM, réalisée dans une configuration de mesure qui diffère d'une configuration de référence uniquement par le capteur utilisé. Dans la configuration de mesure le capteur utilisé est un capteur de marque HAMAMTSU et dans la configuration de référence le capteur utilisé est un capteur de marque DEXELA. Le tableau 502 donne, pour chaque PHANTOM, la valeur corrigée/calibrée du paramètre H selon les procédés de l'état de la technique basée directement sur la valeur mesurée. Le tableau 504 donne, pour chaque PHANTOM, la valeur corrigée/calibrée du paramètre H selon l'invention. Le graphique 602 montre la valeur calibrée du paramètre de texture H obtenue avec les techniques de calibration de l'état de l'art en fonction de la valeur mesurée. Le graphique 604 montre la valeur calibrée du paramètre de texture H obtenue avec l'invention en fonction de la valeur mesurée. Les tableaux 502 et 504 montrent que la correction/calibration selon l'invention est meilleure que la correction/calibration selon les procédés de - 12 - l'état de la technique, et fournit des valeurs corrigées plus proches de la configuration de référence. Ceci est également illustré par la figure 6. De plus, en ce qui concerne la valeur SD de déviation standard des niveaux de gris et la valeur MEAN de niveau de gris moyen, les valeurs corrigées obtenues par la correction/calibration selon l'invention sont en adéquation avec les valeurs de SD et de MEAN obtenues dans la configuration de référence. Pour les procédés de l'état de la technique, les valeurs corrigées de SD et de MEAN sont très différentes des valeurs obtenues dans la configuration de référence. Dans les exemples donnés le paramètre de caractérisation considéré est le paramètre de texture H. Bien entendu, l'invention n'est pas limitée au 15 paramètre de texture H et peut être appliquée à tout paramètre mesuré par analyse fractale, tel que par exemple la densité osseuse. L'invention n'est bien sûr pas limitée aux exemples qui viennent 20 d'être décrits. REVENDICATIONS1) Procédé (100) de correction/calibration d'une valeur, dite mesurée, d'un paramètre de caractérisation d'un tissu osseux déterminée par mesure fractale sur une image en niveaux de gris réalisée dans une première configuration, dite de mesure, par rapport à une deuxième configuration, dite de référence, se distinguant de ladite configuration de mesure par une différence affectant la valeur mesurée dudit paramètre de caractérisation, ledit procédé étant caractérisé en ce qu'il comprend les étapes suivantes : - détermination (104), d'une dimension fractale, dite mesurée, de ladite image en niveaux de gris réalisée dans la configuration de mesure, - correction (106) de ladite dimension fractale mesurée en référence à ladite configuration de référence, ladite correction fournissant une valeur dite corrigée de ladite dimension fractale, - détermination (108) d'une valeur, dite corrigée, dudit paramètre de caractérisation à partir de ladite valeur corrigée de ladite dimension fractale. 2) Procédé (100) selon la revendication 1, caractérisé en ce que l'étape de correction de la dimension fractale mesurée, par rapport à la configuration de référence, est réalisée suivant une relation prédéterminée de manière empirique avec plusieurs échantillons de mesure pour lesquels les valeurs du paramètre de caractérisation et/ou les valeurs de la dimension fractale sont connues dans la configuration de mesure et dans la configuration de référence. 3) Procédé (100) selon l'une quelconque des revendications précédentes, caractérisé en ce que la configuration de mesure diffère de la configuration de référence par le détecteur utilisé pour réaliser l'image en niveau de gris. 4) Procédé (100) selon l'une quelconque des revendications précédentes, caractérisé en ce que la configuration de mesure diffère de la configuration- 14 - de référence par la quantité et/ou l'énergie des particules utilisées pour réaliser l'image en niveau de gris. 5) Procédé (200) de détermination d'un paramètre de caractérisation d'un tissu osseux, ledit procédé (200) comprenant les étapes suivantes : acquisition (202-206) d'une image en niveaux de gris dudit tissu osseux, et détermination (100) d'une valeur, dite corrigée/calibrée, dudit paramètre de caractérisation suivant le procédé (100) selon l'une quelconque des revendications précédentes. 6) Procédé (200) selon l'une quelconque des revendications précédentes, caractérisé en ce que le paramètre de caractérisation est la texture osseuse, ou la densité osseuse. 7) Dispositif (300) de correction/calibration d'une valeur, dite mesurée, d'un paramètre de caractérisation d'un tissu osseux déterminée par mesure fractale sur une image en niveaux de gris réalisée dans une première configuration, dite de mesure, par rapport à une deuxième configuration, dite de référence, se distinguant de ladite configuration de mesure par une différence affectant la valeur mesurée dudit paramètre de caractérisation, ledit dispositif (300) étant caractérisé en ce qu'il comprend des moyens (302-308) agencés et configurés pour : - déterminer une dimension fractale de ladite image en niveaux de gris, - corriger ladite dimension fractale en référence à ladite configuration de référence, ladite correction fournissant une valeur dite corrigée de ladite dimension fractale, et - déterminer une valeur, dite corrigée, dudit paramètre de caractérisation à partir de ladite valeur corrigée de ladite dimension fractale. 8) Système (400) de détermination d'un paramètre de caractérisation d'un tissu osseux, ledit système (400) comprenant :- 15 - des moyens (402, 404) d'imagerie pour l'acquisition d'une image en niveaux de gris dudit tissu osseux, et un dispositif (300) de correction/calibration selon la revendication précédente. 9) Système (400) selon la revendication précédente, caractérisé en ce que les moyens d'imagerie (402, 404) comprennent un détecteur (404) numérique ou analogique. 10) Utilisation du dispositif (300) selon la revendication 7, respectivement du système selon l'une quelconque des revendications 8 ou 9, pour la correction/calibration, respectivement la détermination, de la texture osseuse, ou la densité osseuse d'un tissu osseux.15