Brevet FR2957441A1 — publié le 20110916

Procédé d'extraction simultanée de la fonction d'entrée et des paramètres pharmacocinétiques d'un principe actif La présente invention concerne un procédé d'extraction simultanée de la fonction d'entrée et des paramètres pharmacocinétiques d'un principe actif dans une structure d'un organisme vivant à partir d'une image de mesure de la cinétique du principe actif dans la structure obtenue par un appareil d'imagerie, l'image comprenant une pluralité de régions d'intérêt à l'intérieur de chacune desquelles la cinétique du principe actif est mesurée. L'imagerie moléculaire in vivo permet de mesurer des paramètres pharmacocinétiques d'un principe actif dans des animaux et êtres humains intacts. L'imagerie moléculaire englobe des techniques nucléaires et des techniques optiques. Parmi les techniques nucléaires, la tomographie par émission de positons (TEP) ainsi que la tomographie par émission monophotonique (TEMP) sont les techniques les plus sensibles et les seules à offrir des mesures quantitatives. La quantification des mesures avec la TEP ou la TEMP pour l'étude d'un principe actif nécessite l'injection intraveineuse d'un traceur radioactif adapté à l'étude de ce principe actif. Une image de la distribution du traceur, et donc du principe actif, est alors reconstruite à partir d'une séquence temporelle d'images TEP ou TEMP de mesure de la concentration du traceur. Le post-traitement de l'image TEP ou TEMP requiert en général trois étapes pour accéder aux paramètres pharmacocinétiques. Tout d'abord, des régions d'intérêt délimitant les organes pharmacocinétiques, appelés pharmaco-organes, doivent être définis sur la séquence temporelle d'images TEP ou TEMP, ou sur une image morphologique recalée avec la séquence temporelle d'images TEP ou TEMP. Ensuite, les cinétiques des pharmaco-organes, c'est-à-dire les courbes d'activité temporelle du traceur correspondant aux variations au fil du temps de la concentration du traceur, doivent être extraites. Enfin, un modèle physiologique de cinétique peut être défini, sur la base des cinétiques et de la concentration du traceur dans le plasma, permettant le calcul des paramètres pharmacocinétiques intéressants pour le praticien. L'extraction des paramètres pharmacocinétiques en imagerie utilisant un traceur est donc l'ultime étape donnant accès à des informations sur le traitement du traceur par l'organisme, pouvant ainsi renseigner sur l'efficacité d'un médicament ou l'évolution d'une pathologie. On connaît des procédés permettant d'estimer les paramètres pharmacocinétiques et la concentration plasmatique du traceur non métabolisé en fonction du temps, appelée fonction d'entrée, à partir des cinétiques extraites. En effet, la fonction d'entrée est nécessaire à la connaissance des paramètres pharmacocinétiques mais elle n'est pas toujours visible directement dans l'image. Cependant, les paramètres pharmacocinétiques extraits par ces procédés souffrent d'un très fort biais lorsqu'ils ne bénéficient pas d'informations supplémentaires sous la forme de prélèvements veineux ou artériels. Par ailleurs, l'existence d'une telle restriction ne rend ces procédés applicables qu'à un nombre très réduit de traceurs. L'invention a pour but de proposer un procédé à caractère invasif minimal permettant d'extraire directement à partir de l'image, avec le moins de biais possible, les paramètres pharmacocinétiques simultanément avec la fonction d'entrée. A cet effet, l'invention a pour objet un procédé du type précité caractérisé en ce qu'il comprend les étapes suivantes : - détermination d'une pluralité de jeux de cinétiques à partir des cinétiques mesurées, les jeux de cinétiques déterminés étant différents les uns des autres et chaque jeu de cinétiques déterminé comportant plusieurs cinétiques mesurées ; - définition d'au moins un modèle de cinétique, des paramètres du modèle de cinétique défini étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques ; - estimation, pour chaque jeu de cinétiques déterminé, de la fonction d'entrée et des paramètres pharmacocinétiques par résolution du modèle de cinétique défini à partir du jeu de cinétiques déterminé ; et - déduction de la fonction d'entrée finale et des paramètres pharmacocinétiques finaux à partir des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés. Le procédé selon l'invention peut comporter une ou plusieurs des caractéristiques suivantes : - l'étape de détermination de la pluralité de jeux de cinétiques comprend les étapes suivantes : - constitution de tous les jeux de cinétiques possibles à partir des cinétiques mesurées ; et - sélection, parmi les jeux de cinétiques constitués, d'un nombre prédéterminé de jeux de cinétiques, cette étape comprenant les étapes suivantes : • évaluation, pour chaque jeu de cinétiques constitué, d'un critère de sélection de jeux de cinétiques permettant d'évaluer leur risque d'introduire des erreurs dans l'estimation de la fonction d'entrée et des paramètres pharmacocinétiques ; • classement des jeux de cinétiques constitués selon une règle du critère de sélection maximal ; et • choix du nombre prédéterminé des premiers jeux de cinétiques classés ; - le critère de sélection de jeux de cinétiques permet de sélectionner les jeux de cinétiques dont les cinétiques mesurées sont différentes les unes des autres en termes de forme et/ou d'amplitude, et/ou dont les cinétiques mesurées sont peu affectées par l'effet de volume partiel et/ou dont les cinétiques mesurées présentent un fort rapport signal sur bruit ; - le modèle de cinétique défini prend en compte la métabolisation du principe actif dans la structure de l'organisme vivant imagé ; l'étape d'estimation de la fonction d'entrée et des paramètres pharmacocinétiques est réalisée plusieurs fois pour chaque jeu de cinétiques déterminé de manière à trouver les différentes solutions du modèle de cinétique défini ; - les jeux de cinétiques déterminés comprennent des jeux de cinétiques comportant le même nombre de cinétiques mesurées ; - les jeux de cinétiques déterminés comprennent des jeux de cinétiques comportant un nombre différent de cinétiques mesurées ; - l'étape d'estimation de la fonction d'entrée et des paramètres pharmacocinétiques comprend la minimisation d'une fonction de coût permettant d'évaluer la dissemblance entre le jeu de cinétiques déterminé et le modèle de cinétique défini ; - l'étape d'estimation de la fonction d'entrée et des paramètres pharmacocinétiques est réalisée suivant une boucle itérative dans laquelle les fonctions d'entrée estimées et les paramètres pharmacocinétiques estimés lors d'une itération permettent de déterminer l'initialisation des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés à l'itération suivante et/ou d'apporter des contraintes sur les fonctions d'entrée estimées et les paramètres pharmacocinétiques estimés à l'itération suivante ; - l'étape de définition d'au moins un modèle de cinétique comprend la définition d'une pluralité de modèles de cinétique, des paramètres des modèles de cinétique définis étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques ; - l'étape d'estimation de la fonction d'entrée et des paramètres pharmacocinétiques comprend, pour chaque jeu de cinétiques déterminé, la résolution des modèles de cinétique définis à partir du jeu de cinétiques déterminé, et le procédé comprend une étape de choix, parmi la pluralité de modèles de cinétique, d'un modèle de cinétique optimal en termes de compromis entre l'adéquation du modèle de cinétique aux données et le faible nombre de paramètres de la fonction d'entrée et des paramètres pharmacocinétiques à extraire ; - l'étape de définition d'au moins un modèle de cinétique comprend, outre le choix de la forme du modèle de cinétique, l'introduction dans le modèle de cinétique de contraintes résultant de connaissances a priori non liées à des prélèvements artériels effectués sur l'organisme vivant imagé, les contraintes permettant de fixer au moins un paramètre du modèle de cinétique tel que les fractions vasculaires des régions d'intérêt ; - l'étape d'estimation de la fonction d'entrée et des paramètres pharmacocinétiques comprend l'introduction dans la résolution du modèle de cinétique de contraintes résultant de connaissances a priori non liées à des prélèvements artériels effectués sur l'organisme vivant imagé, les contraintes permettant d'ajouter des termes connus dans la résolution du modèle de cinétique tels que des termes prenant en compte des connaissances a priori sur les fractions vasculaires des régions d'intérêt et/ou sur le volume de distribution du principe actif sur une population donnée ; et - l'image de mesure de la cinétique du principe actif dans la structure est obtenue à partir d'une séquence temporelle d'images tridimensionnelles obtenues par tomographie par émission de positons ou par tomographie par émission monophotonique ; et - l'étape de déduction de la fonction d'entrée finale et des paramètres pharmacocinétiques finaux est réalisée en calculant la moyenne ou la médiane des paramètres des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés. L'invention a également pour objet un support de stockage d'information comprenant un code pour extraire simultanément la fonction d'entrée et les paramètres pharmacocinétiques d'un principe actif dans une structure d'un organisme vivant à partir d'une image de mesure de la cinétique du principe actif dans la structure obtenue par un appareil d'imagerie, l'image comprenant une pluralité de régions d'intérêt à l'intérieur de chacune desquelles la cinétique du principe actif est mesurée, caractérisé en ce que le code comprend des instructions pour : - déterminer une pluralité de jeux de cinétiques à partir des cinétiques mesurées, les jeux de cinétiques déterminés étant différents les uns des autres et chaque jeu de cinétiques déterminé comportant plusieurs cinétiques mesurées ; - définir au moins un modèle de cinétique, des paramètres du modèle de cinétique défini étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques ; - estimer, pour chaque jeu de cinétiques déterminé, la fonction d'entrée et les paramètres pharmacocinétiques par résolution du modèle de cinétique défini à partir du jeu de cinétiques déterminé ; et - déduire la fonction d'entrée finale et les paramètres pharmacocinétiques finaux à partir des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés. L'invention a en outre pour objet un dispositif destiné à extraire simultanément la fonction d'entrée et les paramètres pharmacocinétiques d'un principe actif dans une structure d'un organisme vivant à partir d'une image de mesure de la cinétique du principe actif dans la structure obtenue par un appareil d'imagerie, l'image comprenant une pluralité de régions d'intérêt à l'intérieur de chacune desquelles la cinétique du principe actif est mesurée, caractérisé en ce qu'il comprend : - un appareil d'imagerie; et - un système de traitement de données relié à l'appareil d'imagerie et comprenant : • des moyens pour déterminer une pluralité de jeux de cinétiques à partir des cinétiques mesurées, les jeux de cinétiques déterminés étant différents les uns des autres et chaque jeu de cinétiques déterminé comportant plusieurs cinétiques mesurées ; • des moyens pour définir au moins un modèle de cinétique, des paramètres du modèle de cinétique défini étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques ; • des moyens pour estimer, pour chaque jeu de cinétiques déterminé, la fonction d'entrée et les paramètres pharmacocinétiques par résolution du modèle de cinétique défini à partir du jeu de cinétiques déterminé ; et • des moyens pour déduire la fonction d'entrée finale et les paramètres pharmacocinétiques finaux à partir des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés, sur lesquels : - la Figure 1 est une vue schématique d'un dispositif mettant en oeuvre le procédé selon l'invention ; - la Figure 2 est un organigramme représentant les étapes principales du procédé selon l'invention ; - la Figure 3 est un organigramme représentant les étapes secondaires de la première étape principale de la Figure 2 ; - la Figure 4 est une vue schématique d'un exemple de modélisation compartimentale de processus physiologiques d'interaction entre un traceur et une région d'intérêt ; - la Figure 5 est un organigramme représentant des étapes supplémentaires du procédé selon l'invention ; - les Figures 6 et 7 sont des schémas présentant les résultats obtenus avec le procédé selon l'invention sur des données simulées pour respectivement la fonction d'entrée et les paramètres pharmacocinétiques ; et - les Figures 8 et 9 sont des schémas similaires à ceux respectivement des Figures 6 et 7 présentant les résultats obtenus sur des données expérimentales. Le procédé selon l'invention consiste à extraire simultanément la fonction d'entrée et les paramètres pharmacocinétiques d'un principe actif à partir d'une image de mesure de la cinétique du principe actif. Pour cela, un traceur radioactif, par exemple le FDG, correspondant au principe actif étudié est administré, par exemple par injection intraveineuse, et est véhiculé par le sang jusque dans une structure d'un organisme vivant, par exemple le cerveau d'un être humain. Toute structure d'un organisme est considérée comme un système qui agit sur le comportement du traceur. Le devenir du traceur dans la structure étudiée est ainsi modélisé par un modèle à compartiments, comme cela sera expliqué plus en détail ultérieurement. Dans ce modèle, les différents états du traceur sont modélisés par des compartiments, par exemple un compartiment représentant l'état intact du traceur dans le plasma, un compartiment représentant l'état métabolisé du traceur dans le plasma, un compartiment représentant l'état libre du traceur dans la structure, un compartiment représentant l'état lié non spécifique du traceur dans la structure, un compartiment représentant l'état lié spécifique du traceur dans la structure, etc., dans lesquels le traceur est présent à une certaine concentration. Le modèle est également caractérisé par des vitesses de transfert du traceur entre compartiments. Le nombre de compartiments et leur nature dépendent du type de traceur. La connaissance des concentrations du traceur dans chaque compartiment et des vitesses de transfert du traceur entre compartiments, ces vitesses de transfert formant les paramètres pharmacocinétiques du traceur, permet de remonter à des quantités physiologiques, comme par exemple l'affinité du traceur à sa cible, la consommation de glucose de la structure, etc. La mesure de la concentration radioactive du traceur n'est qu'une conséquence du traitement caché du traceur par la structure et est généralement donnée par la somme pondérée des concentrations du traceur dans les différents compartiments. Afin de remonter des mesures de concentration radioactive du traceur aux paramètres du modèle, il est nécessaire de connaître la concentration de traceur disponible aux structures, qui est la concentration de traceur non métabolisé dans le plasma sanguin en fonction du temps et appelée fonction d'entrée. La Figure 1 représente un dispositif 10 permettant d'extraire simultanément la fonction d'entrée FE et les paramètres pharmacocinétiques PP d'un principe actif au moyen d'un traceur radioactif correspondant au principe actif étudié. Le dispositif 10 comprend un appareil d'imagerie 12 et un système de traitement de données 14 relié à l'appareil d'imagerie 12. Le système de traitement de données 14 comprend des moyens 16 reliés à l'appareil d'imagerie 12 pour déterminer une pluralité de jeux de cinétiques à partir de cinétiques du traceur mesurées dans une structure d'un organisme vivant 18 grâce à l'appareil d'imagerie 12, et des moyens 20 pour définir au moins un modèle de cinétique. Un jeu de I cinétiques peut être défini comme un vecteur à I composantes dans lequel chaque composante est formée par une cinétique mesurée dans une séquence temporelle d'images de mesure au sein d'une région d'intérêt délinéant la structure considérée. Un modèle de cinétique pour une région d'intérêt i peut être défini comme une modélisation de la cinétique du traceur dans cette région, c'est-à-dire une modélisation de la variation en fonction du temps de la concentration du traceur dans la région, certains des paramètres du modèle étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques. Le système de traitement de données 14 comprend des moyens 22 reliés aux moyens de détermination 16 et aux moyens de définition 20 pour estimer la fonction d'entrée et les paramètres pharmacocinétiques, et des moyens 24 reliés aux moyens d'estimation 22 pour déduire la fonction d'entrée finale et les paramètres pharmacocinétiques finaux. Les différents moyens 16, 20, 22, et 24 du système de traitement de données 14 sont mis en oeuvre à l'aide d'un support de stockage d'information 26 comprenant un code 28 comportant des instructions, chaque instruction correspondant respectivement à chacun des moyens 16, 20, 22 et 24. Comme illustré sur la Figure 2, le procédé selon l'invention comprend quatre étapes principales : - une première étape 30 destinée à déterminer la pluralité de jeux de cinétiques et réalisée par les moyens de détermination 16 (Figure 1) ; - une deuxième étape 32 destinée à définir le modèle de cinétique et réalisée par les moyens de définition 20 (Figure 1) ; - une troisième étape 34 destinée à estimer la fonction d'entrée et les paramètres pharmacocinétiques par résolution du modèle de cinétique défini à l'étape 32 à partir des jeux de cinétiques déterminés à l'étape 30, et réalisée par les moyens d'estimation 22 (Figure 1) ; et - une quatrième étape 36 destinée à déduire la fonction d'entrée finale et les paramètres pharmacocinétiques finaux à partir des fonctions d'entrée et des paramètres pharmacocinétiques estimés à l'étape 34, et réalisée par les moyens de déduction 24 (Figure 1). Les quatre étapes principales 30, 32, 34 et 36 du procédé selon l'invention sont expliquées ci-après. En référence à la Figure 3, la première étape principale 30 de détermination d'une pluralité de jeux de cinétiques comprend elle-même quatre étapes secondaires. Initialement, un traceur radioactif correspondant au principe actif étudié est injecté au sein d'une structure d'un organisme vivant 18 et une image de mesure de la cinétique du traceur dans la structure est obtenue par l'appareil d'imagerie 12. L'image de mesure est obtenue à partir d'une séquence temporelle d'images tridimensionnelles de mesure de la concentration du traceur dans la structure, elles-mêmes obtenues par tomographie par émission de positons (TEP) ou par tomographie par émission monophotonique (TEMP). A l'aide d'un algorithme quelconque approprié de segmentation d'images, l'image obtenue est partitionnée en une pluralité O de régions d'intérêt. La cinétique du traceur est alors mesurée à l'intérieur de chacune des régions d'intérêt, par exemple par un algorithme quelconque approprié de correction de l'effet de volume partiel, conduisant ainsi à o cinétiques mesurées du traceur. Lors d'une première étape secondaire 40, tous les jeux de cinétiques possibles sont constitués à partir des O cinétiques mesurées de telle sorte que tous les jeux de cinétiques sont différents les uns des autres et que chaque jeu de cinétiques comporte plusieurs cinétiques mesurées, c'est-à-dire au moins deux cinétiques mesurées. 0 Soit I le nombre de cinétiques par jeu, avec 2 s 15O, alors au total 1Co jeux de I cinétiques sont constitués. Une deuxième étape secondaire 42 permet ensuite de sélectionner, parmi les jeux de cinétiques constitués, un nombre prédéterminé de jeux de cinétiques. En effet, soit un jeu r de I cinétiques. L'estimation de la fonction d'entrée et des paramètres pharmacocinétiques sera d'autant meilleure que : - les I cinétiques du jeu r sont différentes les unes des autres, tant d'un point de vue de la forme que de l'amplitude ; - les I cinétiques du jeu r sont peu affectées par l'effet de volume partiel ; et - les I cinétiques du jeu r présentent un fort rapport signal sur bruit. Par exemple, certaines des cinétiques du jeu r peuvent être issues du même ensemble fonctionnel, ou bien de deux ensembles différents mais aux fonctions semblables et ont donc de fortes chances d'être semblables. L'étape de sélection 42 empêche ainsi de sélectionner le jeu r comportant de telles cinétiques. L'étape de sélection 42 comprend elle-même trois sous-étapes. Une première sous-étape 44 consiste à évaluer, pour chaque jeu de cinétiques constitué à l'étape de constitution 40, un critère 4) de sélection de jeux de cinétiques permettant d'évaluer leur risque d'introduire des erreurs dans l'estimation de la fonction d'entrée et des paramètres pharmacocinétiques : i z CO 2 Y(t)ùYi(t) +CO 3 a / ~Y. ù Y. )2+CO41 T Y2(t) Y YI t-1 vi,t avec w; le coefficient de recouvrement de la région i, c'est-à-dire la proportion du signal provenant effectivement de la structure i qui contribue à la mesure moyenne de concentration d'activité dans la région i au temps t, rendant compte de l'affection de la cinétique Y; par l'effet de volume partiel, v;,, la variance de l'erreur commise du fait du bruit statistique dans l'image sur l'estimation de Y;(t), et col, ces, (03i (04, des poids reflétant l'importance relative accordée respectivement à l'affection de la cinétique Y; par l'effet de volume partiel, à la différence entre les formes des cinétiques, à la différence entre les normes des cinétiques et au rapport signal sur bruit des cinétiques. Une deuxième sous-étape 46 consiste ensuite à classer tous les jeux de cinétiques constitués à l'étape de constitution 40 selon une règle du critère de sélection maximal, c'est-à-dire que les jeux de cinétiques constitués sont triés par ordre de 4) décroissant. Une troisième sous-étape 48 consiste à choisir le nombre prédéterminé de jeux de cinétiques parmi les premiers jeux de cinétiques classés à la sous-étape 46. Lorsqu'un utilisateur choisit un nombre maximum de jeux de cinétiques à utiliser pour l'estimation de la fonction d'entrée et des paramètres pharmacocinétiques, par i 1 Iù1 {Y ) _ Iù1i=1 COI + 1 I 2 1 i I -1 j=i+1 exemple 100, alors la sous-étape 48 consiste à choisir les 100 meilleurs jeux de cinétiques selon Y, les autres jeux de cinétiques n'étant pas utilisés pour l'estimation. La deuxième étape principale 32 permet de définir au moins un modèle de cinétique, des paramètres du modèle de cinétique défini étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques. Comme mentionné précédemment, les processus physiologiques d'interaction entre le traceur et une région d'intérêt i sont modélisés par un modèle à compartiments. La Figure 4 représente un exemple d'un modèle à trois compartiments et quatre paramètres pharmacocinétiques. Le traceur traverse la barrière cellulaire avec une constante de vitesse K, et peut de la même façon repasser la barrière pour retourner dans le sang avec une constante de vitesse k2. Il est transformé ou se lie à un récepteur dans les cellules avec une constante de vitesse k3. Si cette transformation ou liaison est réversible, la transformation inverse ou dissociation est caractérisée par une constante de vitesse k4. Les trois compartiments étudiés sont le compartiment plasma-sérum Cp, le compartiment libre Cf; de la région i et le compartiment lié Obi de la région i. Soit Cp(t), Cf;(t) et Cb;(t) les concentrations du traceur dans chacun de ces trois compartiments. Alors les variations des concentrations du traceur dans les différents compartiments s'écrivent : dCf (t) dt KiCp(t) û (k2 + k3)Cf~ (t) + k4Cbi(t) dCbi(t) = k3Cf (t) û k4Cbi(t) dt Soit Ct;(t) la concentration du traceur dans le tissu de la région i : Ct;(t) = Cf;(t) + Cb;(t) La solution de Ct;(t) est donnée par : Ct;(t) = (Cp O h;)(t) où hi(t) est la fonction de réponse de la région i au traceur non métabolisé. Le procédé selon l'invention se propose de prendre en compte la métabolisation éventuelle du traceur, permettant ainsi d'en étendre l'application à un grand nombre de traceurs. En supposant que le traceur est métabolisé par l'organisme vivant imagé 18 en M métabolites radio-marqués, soit Cam(t) la cinétique du métabolite m du traceur dans l'artère, Cp, la concentration plasmatique du métabolite m du traceur et h;,m la fonction de réponse de la région i au métabolite m du traceur. Alors le modèle de cinétique C;(t) dans la région i s'écrit : Ct (t) = Vbti x Ca(t) + (1ù Vbti) x (Cp O hi )(t) + (Vbti x Cam (t) + (1ù Vbti) x (Cpm O hj m)(t )) m=1 avec Vb; la fraction de sang dans la région i et Ca(t) la cinétique du traceur non métabolisé dans l'artère. La troisième étape principale 34 permet d'estimer, pour chaque jeu de cinétiques déterminé à l'étape 30, la fonction d'entrée et les paramètres pharmacocinétiques par résolution du modèle de cinétique défini C;(t) à partir du jeu de cinétiques déterminé. Pour cela, le procédé met en oeuvre un algorithme quelconque approprié permettant de minimiser une fonction de coût d({Y;}1 Ils peuvent être introduits de multiples façons dans le processus d'estimation, en fixant ou en éliminant certains paramètres des équations, ou sous forme de termes supplémentaires dans la fonction de coût d, etc. Une deuxième stratégie consiste à multiplier les estimations : - avec le même jeu de cinétiques, c'est-à-dire que l'étape d'estimation 34 est réalisée plusieurs fois pour chaque jeu de cinétiques déterminé à l'étape de détermination 30, comme représenté par la boucle 52 sur la Figure 2. Cette méthode permet de réduire le biais mais de manière très limitée ; - avec des jeux de cinétiques différents comportant I cinétiques mesurées. En effet, l'étape 30 de détermination d'une pluralité de jeux de cinétiques permet de générer au maximum Co jeux de cinétiques comportant le même nombre I de cinétiques mesurées. Les biais affectant les estimations à partir de ces jeux de cinétiques se compensent ; et - avec des jeux de cinétiques différents comportant des nombres I différents de cinétiques mesurées. Comme vu précédemment, l'étape 30 de détermination d'une 0 pluralité de jeux de cinétiques permet de générer au maximum 1Co jeux de cinétiques. 1-z Ici encore, les biais affectant les estimations à partir de ces jeux de cinétiques se compensent. Une troisième stratégie consiste à procéder à une estimation itérative, illustrée par la boucle 54 sur la Figure 2. A la première itération, E fonctions d'entrée et E fonctions de réponse des régions d'intérêt sont estimées de la manière décrite précédemment. Lors des itérations q suivantes, les Eq_, fonctions d'entrée et les Eq_, fonctions de réponse des régions d'intérêt estimées à l'itération q-1 servent soit à déterminer l'initialisation de chacune des Eq estimations de l'itération q, soit à apporter des contraintes sur l'estimation au moyen de la première stratégie décrite précédemment. Les étapes 32 de définition du modèle de cinétique et 34 d'estimation de la fonction d'entrée et des paramètres pharmacocinétiques ont été décrites en considérant le modèle de cinétique connu. Dans le cas où ce modèle de cinétique n'est pas connu, qu'il s'agisse du comportement pharmacocinétique du traceur non métabolisé et/ou du comportement pharmacocinétique des éventuels métabolites du traceur, le procédé selon l'invention permet d'estimer la fonction d'entrée et les paramètres pharmacocinétiques du traceur en considérant une pluralité de modèles de cinétiques (Figure 5). Pour cela, le procédé comprend la définition d'une pluralité de modèles de cinétique, des paramètres des modèles de cinétique définis étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques, ainsi que la résolution, pour chaque jeu de cinétiques déterminé, des modèles de cinétique définis à partir du jeu de cinétiques déterminé. Plus précisément, le procédé comprend autant d'étapes supplémentaires 32x de définition d'un modèle de cinétique qu'il y a de modèles de cinétique considérés, x variant ainsi entre 1 et X-1, X étant le nombre total de modèles de cinétique considérés. Le procédé comprend également autant d'étapes supplémentaires 34x d'estimation de la fonction d'entrée et des paramètres pharmacocinétiques qu'il y a de modèles de cinétique considérés, x variant ainsi entre 1 et X-1 et chaque étape supplémentaire 34x permettant la résolution d'un modèle de cinétique considéré. Le procédé comprend une étape 60 de choix, parmi la pluralité de modèles de cinétique définis aux étapes 32 et 32x, d'un modèle de cinétique optimal en termes de compromis entre l'adéquation du modèle de cinétique aux données et le faible nombre de paramètres de la fonction d'entrée et des paramètres pharmacocinétiques à extraire. Cette étape de choix 60 est réalisée par des moyens 62 que comporte le dispositif 10 (Figure 1) et reliés aux moyens de déduction 24. Un exemple d'application du procédé selon l'invention dans le cadre d'une étude cérébrale chez l'être humain est décrit ci-après. Une image TEP du cerveau d'un être humain 18 est obtenue par l'appareil d'imagerie 12 et partitionnée en O régions d'intérêt au sein de chacune desquelles la cinétique du traceur est mesurée. L'étape de détermination 30 est effectuée de la manière décrite précédemment et conduit par exemple à 100 jeux de I cinétiques, avec 2 s 15O. Cette étape 30 permet notamment d'éviter de sélectionner un jeu de cinétiques comportant plusieurs cinétiques semblables, par exemple des cinétiques issues de l'hémisphère gauche et de l'hémisphère droit qui sont des ensembles différents mais aux fonctions semblables. L'étape 32 de définition d'un modèle de cinétique propose un modèle à trois compartiments : un compartiment vasculaire et un compartiment tissulaire comprenant lui- même deux compartiments, l'un où le traceur est libre dans le tissu et l'autre où il est lié de manière spécifique. On considère ici un traceur ayant un unique métabolite radio-marqué qui ne passe pas la barrière hémato-encéphalique, comme cela est souvent le cas dans les études pharmacocinétiques cérébrales. Le métabolite n'est alors présent que dans le compartiment vasculaire. On considère également le cas, fréquent, où la concentration plasmatique est égale à la concentration dans l'artère. On suppose que la fraction vasculaire Vb; de la région d'intérêt i a une densité de probabilité Gaussienne connue de moyenne ,1,71 et de variance 6,2b, , mesurées sur un échantillon de la population à laquelle le sujet imagé 18 appartient. On suppose également que le volume de distribution Vd du traceur a une densité de probabilité Gaussienne connue de moyenne ,uvd et de variance 6Ud . En incluant le métabolite dans le modèle de cinétique, ce modèle de cinétique dans la région i s'exprime comme : Cti (t) = Vbti x (Cp (t) + Cal (t» + (1 Vbti) x (Cp O hi )0 où Ca, est la concentration du métabolite dans l'artère. La fonction de réponse de la région i au traceur non métabolisé s'exprime comme : ùb; it ùb; zt h~(t)=ai,le +ai2e où a;,1i a;,2, b;,1 et b;,2 s'expriment comme des fonctions connues des paramètres pharmacocinétiques K1, k2, k3 et k4. On choisit pour fonction d'entrée une fonction paramétrique de la forme : Cp(t) = (AI (t û 2) û A2 û A3 )e-11(t-z) + Ale-12(t-z) + Aie-23 (t-z) où z est le délai entre l'injection du traceur et le début de la cinétique plasmatique dans le champ de vue imagé. La concentration artérielle du traceur métabolisé est donnée par : t t1 :a x (1û e-e-" ) ) Ici, on considère le modèle comme connu donc il n'y a pas d'étapes supplémentaires 32x à effectuer. Dans certaines régions, le traceur peut (a) ne pas avoir de récepteur, voire même (b) s'accumuler sans ressortir de la région. Le cas (a) correspond à fixer les paramètres a2=b2=0, et le cas (b) à fixer a2=b1=b2=0. L'étape 34 d'estimation est ensuite effectuée de la manière suivante. Soit {0j}1 L'étape 36 de déduction de la fonction d'entrée finale et des paramètres pharmacocinétiques finaux est enfin effectuée. Une fois l'ensemble des E estimations faites, la valeur finale de chaque paramètre O; de la fonction d'entrée et ei de la concentration artérielle du métabolite est choisie comme la médiane des E valeurs estimées pour ce paramètre. Les autres paramètres 4;,k sont ré-estimés comme décrit pour l'étape d'estimation 34, mais en fixant les paramètres O; de la fonction d'entrée et 8; de la concentration artérielle auxdites valeurs finales estimées, permettant une meilleure estimation des paramètres pharmacocinétiques. En ce qui concerne la réduction du biais sur les paramètres estimés, le système est ici contraint par le choix d'une forme analytique paramétrique pour la fonction d'entrée, mais également pour la cinétique du métabolite dans l'artère. Les connaissances a priori sur les fractions vasculaires Vb; des régions i et sur le volume de distribution Vd du traceur introduites dans la fonction de coût d sous la forme Vb 77. ù DI des termes s pi, )2 et y AI+A22 ) contraignent la minimisation. 6vb, 6vd S'il est connu que la concentration du traceur s'équilibre entre la veine et l'artère, L on peut ajouter à la fonction de coût d le terme Sxl((Cp(tl)+Cal(tl))ûCv(tl))Z ou Cv i-i 2x 6ä est la concentration veineuse du traceur mesurée à partir du moment où l'équilibre veine-artère est réalisé et mesurée en L points de temps t, différents, 6v est l'incertitude sur les 15 prélèvements veineux due à la l'imprécision statistique du compteur gamma et ç est le poids mesurant l'importance accordée aux prélèvements veineux. Les prélèvements veineux sont à différencier des prélèvements artériels, ces derniers étant beaucoup plus difficiles et inconfortables à réaliser pour le sujet imagé 18. Si l'on est en mesure d'estimer la cinétique artérielle Cameas du traceur (métabolisé 20 et non-métabolisé) directement à partir de l'image et d'une manière fiable, alors la quantité Cp(t)+Cal(t) peut être remplacée par Cameas dans l'expression du modèle de cinétique C;(t), ce qui permet de diminuer de trois le nombre de paramètres à estimer. Une fois les autres paramètres estimés, ces trois paramètres liés à la concentration artérielle du métabolite pourront être estimés en fixant tous les autres paramètres. 25 L'étape d'estimation 34 est réalisée plusieurs fois pour chaque jeu de cinétiques, par exemple cinq fois, afin de trouver les différentes solutions du modèle de cinétique défini qui correspondent aux différents minimas de la fonction de coût d lorsque celle-ci en présente plusieurs. Les jeux de cinétiques déterminés sont tous différents les uns des autres, certains 30 jeux de cinétiques comportant le même nombre I de cinétiques mesurées et d'autres jeux de cinétiques comportant un nombre 1 différent de cinétiques mesurées. 2 Par exemple ici, le nombre I de cinétiques mesurées par jeu de cinétiques varie 0 entre 3 et 0, ce qui correspond à un nombre total maximal de IClo jeux de cinétiques, l=3 ce qui donne 968 jeux de cinétiques pour 0=10 et 32647 jeux de cinétiques pour 0=15. Un nombre aussi important d'estimations n'est pas nécessaire, car le gain en précision en ajoutant un jeu de cinétiques supplémentaire diminue avec le nombre de jeux de cinétiques déjà utilisés. Par exemple, le nombre de jeux de cinétiques peut être limité à au maximum 100. Concernant l'estimation itérative, les paramètres à extraire sont ici estimés en une unique itération. Lorsque le modèle de cinétique n'est pas connu, des modèles de cinétique simplifiés x sont testés au moyen des étapes 32x et34x. L'étape 60 de choix d'un modèle de cinétique optimal permet ensuite de choisir, parmi tous les modèles x testés, le modèle le plus simple qui rend le mieux compte des processus physiologiques sous-jacents aux cinétiques mesurées. Par exemple, le modèle de cinétique choisi pour chaque région i est celui vérifiant : IxK\ K = arg max d ({Y, L Les paramètres fixés du procédé (on rappelle qu'il est mathématiquement impossible d'estimer la fonction d'entrée et/ou les paramètres pharmacocinétiques sans fixer certains paramètres du modèle de cinétique) sont les fractions vasculaires Vb; des régions i, supposées connues pour le cerveau humain et estimées sur une population pour le rongeur, ainsi que le temps de délai z entre l'injection du traceur et la disponibilité du traceur aux régions i. Ce temps z est très faible chez le rongeur, mais plus long chez l'homme. La Figure 6, qui représente les résultats obtenus pour la fonction d'entrée pour différentes valeurs de bruit sur des données simulées type cerveau humain et compare ces résultats à la fonction d'entrée réelle, indique que la fonction d'entrée est bien estimée quelle que soit l'importance du bruit qui affecte les cinétiques d'entrée, et ce sans recourir à des prélèvement veineux et/ou artériels. Les paramètres pharmacocinétiques K1, k2, k3 et k4 sont tirés aléatoirement dans un intervalle cohérent avec des valeurs connues. Pour les données cérébrales, k4 a été choisi petit pour simuler des cinétiques d'examen au traceur FDG. La Figure 7 représente l'erreur relative moyenne sur chacun des paramètres pharmacocinétiques en fonction du niveau de bruit avec des données simulées type cerveau humain. L'erreur relative importante obtenue sur k4 est due à la très faible valeur de ce paramètre, une faible erreur absolue donnant une forte erreur relative. Le macro-paramètre K est un paramètre permettant de mesurer la concentration régionale de glucose : métabolisme régional du glucose (CMRgIu) = kK" Cgp/LC=KxCgp/LC 2 3 où Cgp est la concentration plasmatique de glucose, aisée à obtenir par contrôle de la glycémie durant l'examen, et LC est la constante localisée (ou « lumped constant » en anglais) qui permet de passer du K du modèle du FDG au K du modèle du glucose. La grandeur K = " est donc la principale grandeur d'intérêt. Une erreur de k2 3 10% sur K est raisonnable compte tenu du fait que l'on ne recourt pas à des prélèvements veineux et/ou artériels. Le procédé a également été testé sur des données expérimentales type cerveau humain et avec le traceur FDG et comparé à des prélèvements artériels. Les fractions vasculaires Vb; des régions i sont ici fixées à des valeurs supposées connues, mais devront dans les faits être estimées plus précisément pour améliorer les résultats. Le délai z a été fixé au milieu de l'entrée du signal TEP dans le cerveau. Une méthode d'estimation plus fine de tau doit également être mise en place afin d'améliorer les résultats. La Figure 8, qui représente les résultats obtenus pour la fonction d'entrée sur des données expérimentales type cerveau humain et compare ces résultats à la fonction d'entrée réelle obtenue par des prélèvements artériels, qui constituent la référence, indique que la fonction d'entrée est bien estimée, avec une légère sous-estimation dans sa partie centrale, en particulier au moment du second pic de la fonction d'entrée correspondant au « lavage de la seringue », pic qui peut être évité si l'on n'a pas besoin de connaître précisément l'activité injectée, le « lavage de la seringue » consistant à ajouter du sérum physiologique dans la seringue et à l'injecter pour pousser le traceur résiduel resté dans les cathéters. La Figure 9 représente les résultats obtenus pour chacun des paramètres pharmacocinétiques ainsi que pour le macro-paramètre K avec des données expérimentales type cerveau humain. On ne dispose pas ici de référence, mais les paramètres peuvent être comparés entre régions. Il paraît raisonnable que le facteur K soit très proche pour le cortex, le thalamus, le caudé et le putamen, et plus bas pour la matière blanche et le globus pallidus, avec une valeur intermédiaire pour le cervelet. D'autre part, le fait que K1 soit plus faible pour la matière blanche que pour le cortex paraît également plausible. L'invention propose donc un procédé à caractère invasif minimal qui permet d'extraire directement à partir d'une image de mesure de la cinétique d'un principe actif et avec le moins de biais possible la fonction d'entrée simultanément avec les paramètres pharmacocinétiques de ce principe actif. En minimisant une distance entre plusieurs cinétiques mesurées et les cinétiques obtenues comme solution d'un modèle à compartiments, en multipliant les estimations à partir d'un nombre important de jeux de cinétiques différents et en introduisant des a priori dans l'estimation autres que des prélèvements artériels, le procédé selon l'invention permet de réduire les erreurs commises lors d'une estimation et d'améliorer ainsi la robustesse de l'estimation. De plus, la prise en compte du comportement pharmacocinétique des éventuels métabolites du traceur permet d'appliquer ce procédé à la quasi-totalité des traceurs existants et à venir, tant dans des études corps entier qu'en intracérébral ou d'autres organes. Par ailleurs, le procédé selon l'invention permet de modifier la nature même de l'information apportée par les examens TEP ou TEMP et qui est d'intérêt pour le praticien (médecin, pharmacologue, biologiste, etc.), qui ne donnera plus uniquement des mesures de radioactivité, mais directement les paramètres du modèle qui sous-tendent ces mesures de radioactivité. Plusieurs applications médicales peuvent être envisagées : la pharmacologie, le développement de médicaments, le diagnostic médical, la réduction du nombre d'examens, etc. La pharmacologie est une branche de la recherche médicale dont le but est la compréhension de l'effet du médicament sur l'organisme (pharmacodynamique) et de l'effet de l'organisme sur le médicament (pharmacocinétique). Elle permet ainsi de caractériser des pathologies, des médicaments et des traceurs. L'obtention de paramètres pharmacocinétiques fins est une étape essentielle de la pharmacocinétique. L'enjeu du présent procédé pour la pharmacocinétique est crucial, car les prélèvements artériels sont de moins en moins pratiqués pour des raisons de sécurité et de confort du sujet imagé. Environ 39% des médicaments n'atteignent pas le marché pour des raisons de mauvaise pharmacocinétique. L'obtention de paramètres pharmacocinétiques dès la phase préclinique et à chaque étape de la phase clinique pourrait permettre, avec une grande richesse d'information et des mesures physiologiquement quantitatives, d'écarter ou de modifier aussitôt que possible les molécules candidates. Cet élément est d'autant plus crucial pour le développement des médicaments que les coûts de développement restent actuellement les mêmes alors que de moins en moins de médicaments arrivent sur le marché, du fait des demandes croissantes de preuves scientifiques et de sûreté de la part des organismes officiels. La possibilité d'obtenir des paramètres pharmacocinétiques sans recourir à des prélèvements artériels ouvre à la pharmacocinétique les portes du diagnostic clinique. Les conséquences, en termes d'indices, de sensibilité et de spécificité diagnostiques, sont susceptibles de révolutionner l'usage de la TEP. La directive 86/609/EEC sur la protection des animaux utilisés pour des buts scientifiques expérimentaux et autres demandes à ce que les expérimentations animales soient aussi réduites que possible. Que ce soit chez l'homme ou chez l'animal, le coût d'un examen TEP ou TEMP est de plus extrêmement élevé, et la réduction du nombre d'examens s'accompagnerait d'une réduction substantielle des coûts pour la recherche médicale. Or, il est connu que, pour la mesure d'un même effet d'un médicament avec le même degré de certitude, l'estimation de paramètres pharmacocinétiques permet de réduire d'environ un ordre de grandeur le nombre d'examens nécessaires comparée à une simple mesure de SUV (Standard Uptake Value). REVENDICATIONS1.- Procédé d'extraction simultanée de la fonction d'entrée (FE) et des paramètres pharmacocinétiques (PP) d'un principe actif dans une structure d'un organisme vivant (18) à partir d'une image de mesure de la cinétique du principe actif dans la structure obtenue par un appareil d'imagerie (12), l'image comprenant une pluralité de régions d'intérêt à l'intérieur de chacune desquelles la cinétique du principe actif est mesurée, caractérisé en ce qu'il comprend les étapes suivantes : - détermination (30) d'une pluralité de jeux de cinétiques à partir des cinétiques mesurées, les jeux de cinétiques déterminés étant différents les uns des autres et chaque jeu de cinétiques déterminé comportant plusieurs cinétiques mesurées ; - définition (32) d'au moins un modèle de cinétique, des paramètres du modèle de cinétique défini étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques ; - estimation (34), pour chaque jeu de cinétiques déterminé, de la fonction d'entrée et des paramètres pharmacocinétiques par résolution du modèle de cinétique défini à partir du jeu de cinétiques déterminé ; et - déduction (36) de la fonction d'entrée finale et des paramètres pharmacocinétiques finaux à partir des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés. 2.- Procédé selon la revendication 1, caractérisé en ce que l'étape de détermination (30) de la pluralité de jeux de cinétiques comprend les étapes suivantes : - constitution (40) de tous les jeux de cinétiques possibles à partir des cinétiques mesurées ; et - sélection (42), parmi les jeux de cinétiques constitués, d'un nombre prédéterminé de jeux de cinétiques, cette étape comprenant les étapes suivantes : • évaluation (44), pour chaque jeu de cinétiques constitué, d'un critère de sélection de jeux de cinétiques permettant d'évaluer leur risque d'introduire des erreurs dans l'estimation de la fonction d'entrée et des paramètres pharmacocinétiques ; • classement (46) des jeux de cinétiques constitués selon une règle du critère de sélection maximal ; et • choix (48) du nombre prédéterminé des premiers jeux de cinétiques classés. 3.- Procédé selon la revendication 2, caractérisé en ce que le critère de sélection de jeux de cinétiques permet de sélectionner les jeux de cinétiques dont les cinétiques mesurées sont différentes les unes des autres en termes de forme et/ou d'amplitude,et/ou dont les cinétiques mesurées sont peu affectées par l'effet de volume partiel et/ou dont les cinétiques mesurées présentent un fort rapport signal sur bruit. 4.- Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que le modèle de cinétique défini prend en compte la métabolisation du principe actif dans la structure de l'organisme vivant imagé (18). 5.- Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que l'étape d'estimation de la fonction d'entrée et des paramètres pharmacocinétiques est réalisée plusieurs fois pour chaque jeu de cinétiques déterminé de manière à trouver les différentes solutions du modèle de cinétique défini. 6.- Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que les jeux de cinétiques déterminés comprennent des jeux de cinétiques comportant le même nombre de cinétiques mesurées. 7.- Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que les jeux de cinétiques déterminés comprennent des jeux de cinétiques comportant un nombre différent de cinétiques mesurées. 8.- Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce que l'étape d'estimation (34) de la fonction d'entrée et des paramètres pharmacocinétiques comprend la minimisation d'une fonction de coût (d) permettant d'évaluer la dissemblance entre le jeu de cinétiques déterminé et le modèle de cinétique défini. 9.- Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que l'étape d'estimation (34) de la fonction d'entrée et des paramètres pharmacocinétiques est réalisée suivant une boucle itérative (54) dans laquelle les fonctions d'entrée estimées et les paramètres pharmacocinétiques estimés lors d'une itération permettent de déterminer l'initialisation des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés à l'itération suivante et/ou d'apporter des contraintes sur les fonctions d'entrée estimées et les paramètres pharmacocinétiques estimés à l'itération suivante. 10.- Procédé selon l'une quelconques des revendications 1 à 9, caractérisé en ce que : - l'étape de définition (32) d'au moins un modèle de cinétique comprend la définition (32x) d'une pluralité de modèles de cinétique, des paramètres des modèles de cinétique définis étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques ;- l'étape d'estimation (34) de la fonction d'entrée et des paramètres pharmacocinétiques comprend, pour chaque jeu de cinétiques déterminé, la résolution (34x) des modèles de cinétique définis à partir du jeu de cinétiques déterminé, et le procédé comprend une étape de choix (60), parmi la pluralité de modèles de cinétique, d'un modèle de cinétique optimal en termes de compromis entre l'adéquation du modèle de cinétique aux données et le faible nombre de paramètres de la fonction d'entrée et des paramètres pharmacocinétiques à extraire. 11.- Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que l'étape de définition (32) d'au moins un modèle de cinétique comprend, outre le choix de la forme du modèle de cinétique, l'introduction dans le modèle de cinétique de contraintes résultant de connaissances a priori non liées à des prélèvements artériels effectués sur l'organisme vivant imagé (18), les contraintes permettant de fixer au moins un paramètre du modèle de cinétique tel que les fractions vasculaires (Vb;) des régions d'intérêt. 12.- Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que l'étape d'estimation (34) de la fonction d'entrée et des paramètres pharmacocinétiques comprend l'introduction dans la résolution du modèle de cinétique de contraintes résultant de connaissances a priori non liées à des prélèvements artériels effectués sur l'organisme vivant imagé (18), les contraintes permettant d'ajouter des termes connus dans la résolution du modèle de cinétique tels que des termes prenant en compte des connaissances a priori sur les fractions vasculaires (Vb;) des régions d'intérêt et/ou sur le volume de distribution (Vd) du principe actif sur une population donnée. 13.- Procédé selon l'une quelconque des revendications 1 à 12, caractérisé en ce que l'image de mesure de la cinétique du principe actif dans la structure est obtenue à partir d'une séquence temporelle d'images tridimensionnelles obtenues par tomographie par émission de positons (TEP) ou par tomographie par émission monophotonique (TEMP). 14.- Support de stockage d'information (26) comprenant un code (28) pour extraire simultanément la fonction d'entrée (FE) et les paramètres pharmacocinétiques (PP) d'un principe actif dans une structure d'un organisme vivant (18) à partir d'une image de mesure de la cinétique du principe actif dans la structure obtenue par un appareil d'imagerie (12), l'image comprenant une pluralité de régions d'intérêt à l'intérieur de chacune desquelles la cinétique du principe actif est mesurée, caractérisé en ce que le code (28) comprend des instructions pour :- déterminer une pluralité de jeux de cinétiques à partir des cinétiques mesurées, les jeux de cinétiques déterminés étant différents les uns des autres et chaque jeu de cinétiques déterminé comportant plusieurs cinétiques mesurées ; - définir au moins un modèle de cinétique, des paramètres du modèle de cinétique défini étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques ; - estimer, pour chaque jeu de cinétiques déterminé, la fonction d'entrée et les paramètres pharmacocinétiques par résolution du modèle de cinétique défini à partir du jeu de cinétiques déterminé ; et - déduire la fonction d'entrée finale et les paramètres pharmacocinétiques finaux à partir des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés. 15.- Dispositif (10) destiné à extraire simultanément la fonction d'entrée (FE) et les paramètres pharmacocinétiques (PP) d'un principe actif dans une structure d'un organisme vivant (18) à partir d'une image de mesure de la cinétique du principe actif dans la structure obtenue par un appareil d'imagerie (12), l'image comprenant une pluralité de régions d'intérêt à l'intérieur de chacune desquelles la cinétique du principe actif est mesurée, caractérisé en ce qu'il comprend : - un appareil d'imagerie (12) ; et - un système de traitement de données (14) relié à l'appareil d'imagerie (12) et comprenant : • des moyens (16) pour déterminer une pluralité de jeux de cinétiques à partir des cinétiques mesurées, les jeux de cinétiques déterminés étant différents les uns des autres et chaque jeu de cinétiques déterminé comportant plusieurs cinétiques mesurées ; • des moyens (20) pour définir au moins un modèle de cinétique, des paramètres du modèle de cinétique défini étant formés par les paramètres de la fonction d'entrée et les paramètres pharmacocinétiques ; • des moyens (22) pour estimer, pour chaque jeu de cinétiques déterminé, la fonction d'entrée et les paramètres pharmacocinétiques par résolution du modèle de cinétique défini à partir du jeu de cinétiques déterminé ; et • des moyens (24) pour déduire la fonction d'entrée finale et les paramètres pharmacocinétiques finaux à partir des fonctions d'entrée estimées et des paramètres pharmacocinétiques estimés.