FR 3002732 A1 20140905 FR 1351855 A 20130301 - 1 - DESCRIPTION La présente invention concerne le domaine des techniques d'imagerie et de représentations exploratoires de sujets, plus particulièrement de sujets humains, notamment en temps réel ou quasi temps réel, et a pour objet un procédé automatique de détermination prédictive de la position de la peau d'un sujet, et des déformations/déplacements de cette dernière du fait de la respiration. Actuellement, dans de nombreux contextes médicaux, les mouvements respiratoires constituent un paramètre qu'il est nécessaire de prendre en compte. Par exemple, en radiothérapie, le volume cible planifié (PTV) est volontairement surestimé de manière à tenir compte des déplacements de la tumeur cible liés aux mouvements inconnus dus à la respiration. Si cela permet d'assurer une irradiation complète de la tumeur, cela augmente en même temps l'irradiation de tissus sains. Dans le contexte de ponctions percutanées en radiologie interventionnelle, pour atteindre sa cible, le radiologue doit bloquer la respiration du patient dans la même position que celle de l'image médicale TomoDensitométrique, IMR ou échographique utilisée pour repérer la cible. De ce fait, la durée de l'opération est nettement plus longue que si la position de la tumeur était parfaitement connue. Ces deux exemples illustratifs expliquent la nécessité de développer des méthodes précises de simulation du mouvement des organes induit par la respiration. Pour certaines de ces méthodes, la connaissance en temps réel de la position de la peau du patient est obligatoire. En fait, la position de la peau peut être utilisée pour piloter un modèle numérique qui prédit la position des structures d'intérêt localisées dans les cavités abdominales et/ou thoraciques (cf par exemple : "A realtime predicitive simulation of abdominal viscera positions during quiet free breathing" - "Simulation prédictive en temps réel des positions des viscères abdominaux durant une respiration libre non forcée" - A. Hostettler et al., Progress in Biophysics and Molecular Biology, 103 (2010), p. 169-184, ELSEVIER). - 2 - S'il existe aujourd'hui des méthodes dites de "motion capture" permettant de fournir une modélisation temps réel de la peau d'un sujet en mouvement, elles nécessitent cependant toutes des ressources matérielles importantes, rendant leur implémentation incompatible avec une mise en place dans une salle d'opération. De plus, ces méthodes connues ne fournissent généralement qu'une estimation ou mesure d'une zone fractionnaire de la peau du patient et non pas de la totalité d'une zone d'intérêt (par exemple sur sensiblement 360° et sur une longueur importante dans la direction crânio-caudale). Concrètement, cette information de surface est difficile à obtenir : durant la respiration, la peau se déplace presque partout et dans différentes directions, seule une partie de la surface du dos qui est en contact avec la table restent immobile. De surcroît, dans de nombreuses applications, il est nécessaire de suivre la position de la peau avec une erreur ne devant pas dépasser l'ordre du millimètre. Plus précisément, on connaît actuellement essentiellement trois méthodes différentes aptes à fournir un résultat exploitable, mais présentant toutes des inconvénients majeurs, à savoir : a) Projection de lumière structurée, suivi par caméra et 20 reconstruction de la surface de la peau : Dans ce cas, une projection de lumière structurée doit être effectuée à partir de plusieurs sources de manière à couvrir toute la surface de la peau du patient. Cela pose un problème de recouvrement à la jonction des multiples projections. Les caméras qui permettront de suivre les 25 mouvements de la peau doivent également être nombreuses afin de garantir une visualisation complète de ces projections. Un algorithme de traitement d'image robuste doit être utilisé pour calculer en temps réel la position de la peau qui est alors décrite comme un nuage de points plus ou moins ordonnés. Enfin, ce nuage de points doit le plus souvent être ré- 30 échantillonné de manière à le rendre exploitable. Quoiqu'il en soit, le patient étant couché sur une table, cette méthode ne permet de localiser que la partie de la peau visible avec les caméras et sur laquelle est projetée la lumière structurée. b) Suivi de marqueurs collés sur la peau du patient et calcul de 35 la position de la peau par des méthodes d'interpolation : Cette méthode représente une simplification de la méthode précédente en réduisant de façon significative le nombre de points suivis. - 3 - Bien que plus rapide, elle conserve les mêmes limites : les caméras doivent être suffisamment nombreuses pour suivre tous les marqueurs visibles. De ce fait, seule la partie de la peau visible aux caméras et sur laquelle sont placés les marqueurs pourra être suivie Enfin, l'interpolation sera moins précise que la méthode précédente. c) Suivi de marqueurs électromagnétiques collés sur la peau du patient (comme par exemple le système Aurora, développé par la société NDI) et calcul de la position de la peau par des méthodes d'interpolation : Dans ce cas, il n'y a pas besoin de caméra, mais les problèmes d'interpolation sont similaires. De surcroît, la précision de la mesure de la position fournie par les marqueurs électromagnétiques est plus faible que les méthodes précédentes (de 1 à 5 mm), l'erreur de mesure étant en plus potentiellement augmentée par la présence d'objets métalliques dans les champs de mesure. Cette dernière limitation ajoutée à la limitation actuelle du nombre de marqueurs simultanément utilisables ainsi qu'à la fréquence limitée à 50Hz, rend de fait difficile la mise en oeuvre de cette méthode en routine clinique. Compte tenu de ce qui précède, le but principal de l'invention 20 est de proposer une solution de modélisation des mouvements de la peau d'un sujet humain qui fournisse une simulation tridimensionnelle fiable et précise d'au moins une zone d'intérêt, qui évolue en temps réel, et ne nécessite qu'un minimum de mesures et une puissance de traitement limitée. A cet effet, l'invention a pour objet un procédé automatique de 25 détermination prédictive, par exemple sous la forme d'une simulation numérique, de la position et des déplacements de la peau d'un sujet humain ou animal dans une zone d'intérêt localisée au niveau du thorax et/ou de l'abdomen, ledit sujet étant allongé sur son dos et respirant librement ou de manière assistée, ledit procédé consistant, de manière préalable, à acquérir 30 ou déterminer plusieurs positions et configurations du profil de la peau dans des plans axiaux consécutifs prédéterminés, sensiblement perpendiculaires à l'axe crânio-caudal (du sujet / de la table de l'imageur) et répartis le long de cet axe à des emplacements déterminés, ceci à des instants successifs donnés, durant des phases d'inspiration et d'expiration du sujet, 35 procédé caractérisé en ce qu'il consiste, durant une phase préalable, pour chaque plan axial considéré, à construire au moins un modèle numérique déformable à partir d'informations relatives aux - 4 - différents profils de peau acquis ou déterminés précédemment, à relever ensuite, de manière répétitive, la position réelle d'un point ou d'une zone ponctuelle spécifique de la peau du sujet au niveau de chacun des plans axiaux précités, préférentiellement repéré(e) à l'aide d'un marqueur visible dans le procédé d'imagerie utilisé et dont la position est modifiée de manière significative durant les phases d'inspiration et d'expiration, et à fournir, sensiblement en temps réel, une simulation du profil de la peau dans chaque plan axial, en fonction de la position réelle relevée et en exploitant le modèle numérique déformable correspondant, ainsi qu'une représentation tridimensionnelle évolutive de la peau au niveau de la zone d'intérêt, par interpolation entre les différents plans axiaux. Le principe de base de l'invention repose sur l'exploitation du constat surprenant fait par les inventeurs que, dans le cas d'un sujet humain couché sur le dos et respirant normalement (de manière non forcée), librement ou de manière assistée, une modélisation tridimensionnelle fiable et précise de la peau du sujet (au niveau de la région du thorax et de l'abdomen), permettant une prédiction en temps réel des mouvements de la peau, dans les trois dimensions, pouvait être réalisée en ne relevant de manière continue que les informations relatives au déplacement et à la position d'un nombre restreint de points choisis de la peau, et en exploitant des données de modélisation acquises précédemment. Plus précisément, le repérage et le suivi en temps réel de la position d'un point de la peau permet, en exploitant un profil de la peau acquis précédemment dans un plan axial donné et comprenant ce point, de prédire les déformations de ce profil de la peau dans ce plan axial au cours du temps en fonction de la respiration du sujet. L'invention sera mieux comprise, grâce à la description ci-après, qui se rapporte à un mode de réalisation préféré, donné à titre d'exemple non limitatif, et expliqué avec référence aux dessins 30 schématiques annexés, dans lesquels : les figures lA à 1C sont des vues, obtenues par IRM, selon le plan sagittal de la zone d'intérêt d'un sujet humain pourvu d'un marqueur filaire selon l'invention, respectivement à l'état expiré (Fig. 1A), à l'état inspiré (Fig. 1B) et par superposition des différentes positions prises par le 35 marqueur (par le biais de sa courbe médiane) dans les phases d'inspiration et d'expiration (Fig. 1C) ; - 5 - les figures 2A à 2D sont des vues, obtenues par IRM, selon un plan axial au niveau d'une zone d'intérêt d'un sujet humain pourvu du marqueur des figures 1, illustrant le profil de la peau déterminé à l'état expiré (Fig. 2A) et inspiré (Fig. 2B), le calcul de la courbe de référence à partir des deux courbes de profils extraites des figures 2A et 2B (Fig. 2C) et la génération d'un champ de déformation à partir des couples de points appartenant aux deux courbes de profil de peau des figures 2A et 2B (Fig. 2D) ; la figure 3 est une vue similaire à celle de la figure 1B avec une 10 mise en évidence de la courbe médiane du marqueur ; les figures 4A à 4C sont des vues de détail en coupe axiale du marqueur illustrant les étapes de détermination des contours et de l'inclinaison du marqueur ; les figures 5A à 5D sont des vues en coupe axiale illustrant les 15 étapes de segmentation et d'extraction du profil de peau ; la figure 6 est une vue de détail en coupe axiale du marqueur et de son environnement immédiat, obtenue par superposition et illustrant différentes positions du centre de gravité du marqueur et l'hystérésis entre les positions prises durant l'inspiration et celles prises durant l'expiration ; 20 les figures 7A et 7B sont des représentations schématiques, respectivement de dessus et de côté, du tronc et de la tête du sujet humain équipé d'un marqueur selon l'invention pour le suivi des mouvements respiratoires, et, les figures 8A à 8C sont des représentations d'une simulation 25 numérique tridimensionnelle de la zone d'intérêt d'un patient humain, asservie à la respiration, obtenue et utilisée dans le cadre de l'invention, les figures 8B et 8C montrant simultanément des vues acquises par l'imagerie dans des plans particuliers. L'invention concerne un procédé automatique de détermination 30 prédictive, par exemple sous la forme d'une simulation numérique, de la position et des déplacements de la peau 2 d'un sujet 1 humain ou animal dans une zone d'intérêt localisée au niveau du thorax et/ou de l'abdomen. Dans le cadre de l'invention, on considère que le sujet 1 est allongé sur son dos et respire librement ou de manière assistée. 35 Le procédé consiste, de manière préalable, à acquérir ou déterminer plusieurs positions et configurations du profil 3 de la peau 2 dans des plans axiaux Pa consécutifs prédéterminés, sensiblement - 6 - perpendiculaires à l'axe crânio-caudal X et répartis le long de cet axe à des emplacements déterminés, ceci à des instants successifs donnés, durant des phases d'inspiration et d'expiration du sujet. Conformément à l'invention, ce procédé consiste également, 5 durant une phase préalable, pour chaque plan axial Pa considéré, à construire au moins un modèle numérique déformable à partir d'informations relatives aux différents profils de peau 3 acquis ou déterminés précédemment, à relever ensuite, de manière répétitive, la position réelle d'un point ou d'une zone ponctuelle 4 spécifique de la peau 2 10 du sujet 1 au niveau de chacun des plans axiaux Pa précités, préférentiellement repéré(e) à l'aide d'un marqueur 5 visible dans le procédé d'imagerie utilisé et dont la position est modifiée de manière significative durant les phases d'inspiration et d'expiration, et à fournir, sensiblement en temps réel, une simulation 3' du profil de la peau dans chaque plan axial Pa, 15 en fonction de la position réelle relevée et en exploitant le modèle numérique déformable correspondant, ainsi qu'une représentation tridimensionnelle 2' évolutive de la peau 2 au niveau de la zone d'intérêt, par interpolation entre les différents plans axiaux Pa. De manière avantageuse, et afin de faciliter le suivi du 20 déplacement de la peau et d'en augmenter la précision, il est prévu que les points ou zones ponctuelles 4 de la peau 2, dont la position est suivie dans le temps, soient disposé(e)s dans un plan parallèle proche du ou confondu avec le plan sagittal médian PSM du sujet 1, en étant préférentiellement espacé(e)s mutuellement d'une distance de l'ordre du centimètre. 25 En accord avec le principe de base de l'invention, le procédé consiste, pour chaque plan axial Pa concerné, à générer sensiblement en temps réel, une courbe 3' représentant le profil 3 de la peau 2 dans le plan axial Pa considéré, en appliquant à une courbe de base Cb de profil de peau, déterminée dans un état extrême d'inspiration ou d'expiration du sujet, une 30 déformation, par exemple par déplacement d'au moins certains points espacés de la courbe de base, fonction de la position mesurée d'un point 4 de la peau 2 du sujet 1 situé dans le plan axial Pa considéré. La courbe de base utilisée correspond préférentiellement au profil axial de la peau relevé à l'état d'expiration maximale du sujet 2. 35 Conformément à une autre caractéristique de l'invention, illustrée par les figures 2A à 2D, chaque modèle numérique déformable de profil de peau, associé à l'un des plans axiaux Pa prédéterminés, est obtenu - 7 - par génération d'une courbe de référence Cr située à mi-distance entre les deux courbes de profils de peau correspondant aux états extrêmes d'inspiration et d'expiration, à déterminer les points d'intersection, avec les deux courbes précitées, d'une droite normale à la courbe de référence Cr et passant par un point de cette courbe, ce pour une pluralité de points espacés le long de cette dernière, et à rassembler de manière ordonnée ces couples de points d'intersection pour constituer un champ de déformation correspondant au déplacement de la peau 2 dans le plan axial Pa considéré entre les deux états extrêmes précités, en fonction de l'état respiratoire du sujet 1. Comme indiqué précédemment, le modèle tridimensionnel de simulation globale de la peau est composé d'une pluralité de modèles axiaux élémentaires bidimensionnels et indépendants de la peau, établis dans des plans axiaux PA consécutifs le long de l'axe crânio-caudal. De manière pratique et une fois la peau segmentée dans les images correspondant respectivement aux positions extrêmes d'inspiration et d'expiration du sujet (Fig. 2A et 2B), chaque modèle axial peut être établi en exécutant les quatre étapes opératoires suivantes (voir Fig. 2C et 2D) : 1) Création d'une courbe de référence Cr située exactement à 20 mi-distance entre les deux courbes correspondant aux deux positions extrêmes précitées. Cette étape peut être réalisée en utilisant une carte de distance correspondant à la différence entre le masque d'image relatif à la position de la peau à l'état inspiré et le masque d'image relative à la position de la peau à l'état expiré. 25 2) Ré-échantillonnage de la courbe de référence axiale Cr de manière à obtenir un ensemble de points consécutifs espacés de manière régulière. 3) Génération d'une droite normale à la courbe de référence Cr au niveau de chaque point de cette courbe ré-échantillonnée. 30 4) Détermination des deux points d'intersection correspondant à l'intersection de cette droite normale avec les courbes des états inspiré et expiré. On obtient ainsi un champ de déformation reflétant le mouvement de la surface axiale de la peau entre les deux images originales des états inspiré et expiré. 35 Ainsi, chaque ensemble d'informations résultant d'une acquisition 2D permet la construction d'un modèle axial local de la peau du patient. Compte tenu de la variation douce de la forme de la peau d'un sujet - 8 - classique, la génération de modèles axiaux de la peau séparés entre eux d'environ 1 cm permet d'obtenir par interpolation un modèle tridimensionnel fidèle et précis de la peau du patient au niveau de la zone d'intérêt. Afin de prendre en compte les faits d'hystérésis influençant de manière différenciée les mouvements de la peau durant les cycles consécutifs de respiration, il est prévu dans le cadre du procédé selon l'invention de construire deux modèles numériques déformables de profil de peau pour chaque plan axial Pa, à savoir un modèle pour les phases d'inspiration et un modèle pour les phases d'expiration, et à détecter la phase respiratoire en cours lors de l'élaboration des courbes simulées 3' de profil de peau dans les différents plans axiaux Pa, ce par l'intermédiaire d'un capteur spécifique ou en évaluant le sens de déplacement des points 4. En vue d'aboutir à un repérage automatique aisé et fiable des points 4 de la peau 2 dont les mouvements sont pris en compte dans le cadre du procédé en tant que paramètres d'entrée, les positions de ces points 4, situés dans le plan sagittal médian PSM du sujet 1 et dans les plans axiaux Pa prédéterminés sont déterminées par l'intermédiaire d'au moins un marqueur 5, préférentiellement un marqueur allongé continu, reposant directement ou indirectement sur la peau 2 du patient 1, épousant la surface locale de la peau 2 et disposé dans le plan sagittal médian PSM. Comme le montrent les figures 1 à 4, 6 et 7, le marqueur 5 mis en oeuvre consiste en un corps souple et allongé, réalisé en un matériau particulièrement bien visible dans le système d'imagerie utilisé, s'étendant au moins sur la longueur de la zone d'intérêt dans la direction de l'axe crânio-caudal X, présentant une section de forme et de dimension minimale déterminées et situé à une distance déterminée de la peau 2 du sujet 1, en reposant sur une bande 6 d'un matériau souple et sensiblement non visible dans le système d'imagerie utilisé. Connaissant les dimensions du marqueur 5 et celles de la bande de matériau 6 et en vue de faciliter le repérage fiable de la position du marqueur 5, le procédé peut avantageusement consister à déterminer et à relever, dans le plan sagittal, la courbe Cm correspondant à l'axe médian du marqueur 5 allongé, par exemple formé par la succession des différents centres de gravité en section dudit marqueur 5, et/ou dans les différents plans axiaux Pa concernés, les centres de gravité locaux dudit marqueur 5, et à déterminer les positions des points spécifiques 4 de la peau 2 du sujet 1 - 9 - dont les mouvements sont pris en compte, en exploitant les informations relatives, d'une part, à cette courbe Cm ou à ces centres de gravité locaux, d'autre part, à la section et aux dimensions en section du marqueur 5 et, enfin, à l'épaisseur de la bande de matériau 6 intercalaire. En accord avec une implémentation pratique préférée de l'invention, ressortant des figures annexées, il peut être prévu que, lorsque le procédé d'imagerie auquel est soumis le sujet 1 est basé sur la résonance magnétique, le marqueur allongé 5 est constitué par une barre formée d'eau gélifiée, avec une section carrée ou rectangulaire et une longueur suffisante pour couvrir au moins la zone du sujet dont on souhaite suivre les mouvements induits par la respiration dudit sujet 1 (zone d'intérêt). Dans ce contexte, la segmentation de la peau dans des images sagittales obtenues par résonance magnétique est avantageusement réalisée en exploitant la très bonne visibilité du marqueur 5 dans les images IRM et son contraste par rapport à la bande textile 6 par l'intermédiaire de laquelle il repose sur la peau du sujet (les épaisseurs dudit marqueur et de ladite bande étant des données connues). Plus précisément, cette phase de segmentation sagittale de la peau peut être décomposée en quatre étapes principales, à savoir : 1) Une image de travail 2D sur-échantillonnée et créée à partir de l'image 2D originelle (prise sur le sujet) en mettant en oeuvre une interpolation cubique (cette étape augmente la précision de la détermination de la position du marqueur 5). De manière pratique, et compte tenu des déformations et distorsions géométriques affectant les images RMN 2D (au niveau des deux côtés opposés des extrémités de la zone explorée - voir par exemple Fig. 3), seule la partie centrale des images initiales est exploitée (environ 60 %). 2) Compte tenu des différences de contraste entre des images 2D + t consécutives, il est avantageux de déterminer une image seuil permettant la détermination de la valeur seuil correspondant à la position du marqueur 5. On recherche ensuite dans chaque droite verticale des différentes images la première suite de pixels consécutifs présentant une longueur au moins supérieure à environ 80 % de l'épaisseur du marqueur 5 et on calcule ensuite le centre du marqueur 5 le long de chaque droite verticale (permettant de déterminer la position du marqueur le long de cette droite). La meilleure valeur seuil correspond à celle qui minimise la - 10 - distance moyenne entre la position simulée du marqueur et la position extraite par segmentation des images originelles. 3) On localise ensuite les parties du marqueur allongé 5 qui ne sont pas continues et on modifie localement le tracé de ce dernier en extrapolant les portions régulières du marqueur situé de part et d'autre de la partie discontinue. 4) Une opération de lissage des positions des centres de gravité consécutifs du marqueur est réalisée de manière à corriger les erreurs de segmentation dues à des manques locaux d'information (l'utilisation du centre de gravité permet de limiter l'influence de l'effet de volume partiel). De même, la segmentation du marqueur 5 peut être réalisée dans les plans axiaux PA ayant donné lieu à l'acquisition d'images 2D + t en exécutant les quatre étapes opératoires suivantes : 1) La recherche de la localisation du marqueur 5 à section 15 carrée par exemple est limitée à une zone de recherche définie manuellement en début de procédure et ensuite centrée de manière automatique en se basant sur la position précédente du centre du marqueur. 2) On recherche la meilleure valeur seuil permettant d'obtenir au moins 80 % de la surface théorique de la section du marqueur (on 20 recherche des pixels dont la valeur de gris est inférieure à la valeur seuil). 3) Dès que le masque représentant la section du marqueur est trouvé, on calcule son centre de gravité en utilisant un algorithme connu de calcul de barycentre. 4) On évalue ensuite l'inclinaison du marqueur à partir du 25 masque précédent en effectuant par exemple les opérations suivantes : détermination des quatre portions de disque du plus grand disque centré sur le centre de gravité du marqueur (Fig. 4B) ; détermination des centres de gravité de ces quatre portions de disque (CG1, CG2, CG3, CG4) ; détermination des inclinaisons horizontale et verticale du marqueur par 30 approximation à partir des segments reliant, d'une part, les points CG1 et CG3 et, d'autre part, les points CG2 et CG4. Enfin, la segmentation de la peau dans les images axiales 2D + t est réalisée en utilisant les informations relatives à la position du marqueur 5 fourni précédemment. 35 Cette segmentation axiale peut être réalisée en exécutant les quatre étapes opératoires suivantes (voir Fig. 5A à 5D) : - 11 - 1) Un masque binaire correspondant au seuil de l'image originelle est calculé. 2) On ne conserve que les plus grands éléments connectés du masque et on remplit ce dernier (Fig. 5B et 5C). 3) La courbe correspondant à la limite périphérique du masque est extraite. Puis cette courbe est lissée et complétée pour pallier les erreurs et manques d'information dans l'image originelle. 4) La courbe correspondant à la position de la peau est corrigée à proximité du marqueur 5, connaissant la distance entre le centre du 10 marqueur et la peau. Comme déjà mentionné précédemment, et afin de tenir compte de l'hystérésis affectant la peau entre les phases d'inspiration et d'expiration, il est recommandé selon l'invention de construire un modèle de simulation différent de la peau en fonction de la phase du cycle respiratoire (inspiration 15 et expiration). A cet effet, les courbes illustrant les positions du marqueur 5 doivent être classées en deux catégories relatives l'une à la phase d'inspiration et l'autre à la phase d'expiration. L'effet d'hystérésis est maximal à la fin de la phase d'inspiration 20 et au début de la phase d'expiration, et la discrimination des courbes est triviale à l'exception de l'état d'inspiration extrême (voir Fig. 6). Afin de pouvoir discriminer les positions du marqueur dans cet état critique, on peut avantageusement exploiter simultanément la position du centre du marqueur et la courbure de la peau autour du marqueur (les 25 inventeurs ont constaté que la position latérale du marqueur changeait du fait de la modification de la courbure de la peau). Dans le cadre d'une phase d'acquisition et de traitement préalable du procédé selon l'invention, les courbes représentant le profil 3 de la peau 2 dans des plans axiaux Pa espacés prédéterminés et à différents 30 instants des phases d'inspiration et d'expiration peuvent être obtenues de différentes façons. Ainsi, ces courbes peuvent être : - extraites, par échantillonnage selon ces plans axiaux, de représentations tridimensionnelles de la surface de la peau 2 au niveau de la 35 zone d'intérêt, acquises préalablement et obtenues par repérage et suivi de la peau au moyen de marqueurs (électromagnétiques, optiques ou autres) fixés sur la peau ou par l'intermédiaire de lumière structurée projetée sur la peau ; - 12 - - extraites, par application d'une technique de segmentation, de représentations tridimensionnelles obtenues par un procédé d'imagerie basé sur la résonnance magnétique ou la tomodensitométrie ; - extraites, par segmentation, d'images du type 2D + t 5 effectuées dans lesdits plans axiaux Pa et obtenues par un procédé d'imagerie basé sur la résonnance magnétique ou la tomodensitométrie. Dans le cadre d'une application médicale avantageuse des résultats fournis par l'invention, la représentation tridimensionnelle évolutive 2' de la peau 2 au niveau de la zone d'intérêt est exploitée, en 10 relation avec au moins une image tridimensionnelle de l'abdomen, pour réaliser une simulation prédictive en temps réel des positions des organes et/ou viscères de la cavité abdominale durant des phases de respiration libre normale du sujet 1, en mettant par exemple en oeuvre l'enseignement de la publication citée dans la partie introductive de la présente description. 15 Ainsi, grâce à l'invention, il est possible de mettre à disposition une simulation tridimensionnelle évoluant en temps réel de la peau d'un patient, en ne relevant de manière courante que la position d'un nombre limité de points de la peau, donc en limitant les ressources de traitement nécessaires, tout en fournissant une précision suffisante pour des 20 applications médicales. Une mise en oeuvre pratique de l'invention a permis d'atteindre une précision de 0,5 cm et un rafraîchissement de la simulation à une fréquence de 80 Hz. De plus, par des évaluations comparatives, les inventeurs ont 25 pu estimer que la simulation produite par l'invention présentait une précision au moins équivalente à celle d'une reconstruction de la peau directement à partir d'images acquises par le système d'imagerie. L'invention porte également sur un système d'imagerie, notamment du type IRM, intégrant les moyens techniques matériels et 30 logiciels autorisant la mise en oeuvre du procédé selon l'invention, ces moyens découlant de manière évidente de la description du procédé ci-dessus. Ce système peut, par exemple, comprendre un dispositif d'imagerie par résonance magnétique du type MAGNETOM Aera 1.5 Tesla de la société SIEMENS. 35 Bien entendu, l'invention n'est pas limitée au mode de réalisation décrit et représenté aux dessins annexés. Des modifications restent possibles, notamment du point de vue de la constitution des divers - 13 - éléments ou par substitution d'équivalents techniques, sans sortir pour autant du domaine de protection de l'invention. REVENDICATIONS1. Procédé automatique de détermination prédictive, par exemple sous la forme d'une simulation numérique, de la position et des déplacements de la peau d'un sujet humain ou animal dans une zone d'intérêt localisée au niveau du thorax et/ou de l'abdomen, ledit sujet étant allongé sur son dos et respirant librement ou de manière assistée, ledit procédé consistant, de manière préalable, à acquérir ou déterminer plusieurs positions et configurations du profil de la peau dans des plans axiaux consécutifs prédéterminés, sensiblement perpendiculaires à l'axe crâniocaudal et répartis le long de cet axe à des emplacements déterminés, ceci à des instants successifs donnés, durant des phases d'inspiration et d'expiration du sujet, procédé caractérisé en ce qu'il consiste, durant une phase préalable, pour chaque plan axial (Pa) considéré, à construire au moins un modèle numérique déformable à partir d'informations relatives aux différents profils de peau (3) acquis ou déterminés précédemment, à relever ensuite, de manière répétitive, la position réelle d'un point ou d'une zone ponctuelle (4) spécifique de la peau (2) du sujet (1) au niveau de chacun des plans axiaux (Pa) précités, préférentiellement repéré(e) à l'aide d'un marqueur (5) visible dans le procédé d'imagerie utilisé et dont la position est modifiée de manière significative durant les phases d'inspiration et d'expiration, et à fournir, sensiblement en temps réel, une simulation (3') du profil de la peau dans chaque plan axial (Pa), en fonction de la position réelle relevée et en exploitant le modèle numérique déformable correspondant, ainsi qu'une représentation tridimensionnelle (2') évolutive de la peau (2) au niveau de la zone d'intérêt, par interpolation entre les différents plans axiaux (Pa). 2. Procédé selon la revendication 1, caractérisé en ce que les points ou zones ponctuelles (4) de la peau (2), dont la position est suivie dans le temps, sont disposé(e)s dans un plan parallèle proche du ou confondu avec le plan sagittal médian (PSM) du sujet (1), en étant préférentiellement espacé(e)s mutuellement d'une distance de l'ordre du centimètre. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce qu'il consiste, pour chaque plan axial (Pa) concerné, à générer sensiblement en- 15 - temps réel, une courbe (3') représentant le profil (3) de la peau (2) dans le plan axial (Pa) considéré, en appliquant à une courbe de base (Cb) de profil de peau, déterminée dans un état extrême d'inspiration ou d'expiration du sujet, une déformation, par exemple par déplacement d'au moins certains points espacés de la courbe de base, fonction de la position mesurée d'un point (4) de la peau (2) du sujet (1) situé dans le plan axial (Pa) considéré. 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que chaque modèle numérique déformable de profil de peau, associé à l'un des plans axiaux (Pa) prédéterminés, est obtenu par génération d'une courbe de référence (Cr) située à mi-distance entre les deux courbes de profils de peau correspondant aux états extrêmes d'inspiration et d'expiration, à déterminer les points d'intersection, avec les deux courbes précitées, d'une droite normale à la courbe de référence (Cr) et passant par un point de cette courbe, ce pour une pluralité de points espacés le long de cette dernière, et à rassembler de manière ordonnée ces couples de points d'intersection pour constituer un champ de déformation correspondant au déplacement de la peau (2) dans le plan axial (Pa) considéré entre les deux états extrêmes précités, en fonction de l'état respiratoire du sujet (1). 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'il consiste à construire deux modèles numériques déformables de profil de peau pour chaque plan axial (Pa), à savoir un modèle pour les phases d'inspiration et un modèle pour les phases d'expiration, et à détecter la phase respiratoire en cours lors de l'élaboration des courbes simulées (3') de profil de peau dans les différents plans axiaux (Pa), ce par l'intermédiaire d'un capteur spécifique ou en évaluant le sens de déplacement des points (4). 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que les positions des points (4) de la peau (2) situés dans le plan sagittal médian (PSM) du sujet (1) et dans les plans axiaux (Pa) prédéterminés sont déterminées par l'intermédiaire d'au moins un marqueur (5), préférentiellement un marqueur allongé continu, reposant directement ou indirectement sur la peau (2) du patient (1), épousant la surface locale de la peau (2) et disposé dans le plan sagittal médian (PSM). 7. Procédé selon la revendication 6, caractérisé en ce que le marqueur (5) mis en oeuvre consiste en un corps souple et allongé, réalisé en un matériau particulièrement bien visible dans le système d'imagerie- 16 - utilisé, s'étendant au moins sur la longueur de la zone d'intérêt dans la direction de l'axe crânio-caudal (X), présentant une section de forme et de dimension minimale déterminées et situé à une distance déterminée de la peau (2) du sujet (1), en reposant sur une bande (6) d'un matériau souple et sensiblement non visible dans le système d'imagerie utilisé. 8. Procédé selon l'une quelconque des revendications 6 et 7, caractérisé en ce qu'il consiste à déterminer et à relever, dans le plan sagittal, la courbe (Cm) correspondant à l'axe médian du marqueur (5) allongé, par exemple formé par la succession des différents centres de gravité en section dudit marqueur (5), et/ou dans les différents plans axiaux (Pa) concernés, les centres de gravité locaux dudit marqueur (5), et à déterminer les positions des points spécifiques (4) de la peau (2) du sujet (1) dont les mouvements sont pris en compte, en exploitant les informations relatives, d'une part, à cette courbe (Cm) ou à ces centres de gravité locaux, d'autre part, à la section et aux dimensions en section du marqueur (5) et, enfin, à l'épaisseur de la bande de matériau (6) intercalaire. 9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que les courbes représentant le profil (3) de la peau (2) dans des plans axiaux (Pa) espacés prédéterminés et à différents instants des phases d'inspiration et d'expiration sont extraites, par échantillonnage selon ces plans axiaux, de représentations tridimensionnelles de la surface de la peau (2) au niveau de la zone d'intérêt, acquises préalablement et obtenues par repérage et suivi de la peau au moyen de marqueurs fixés sur la peau ou par l'intermédiaire de lumière structurée projetée sur la peau. 10. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que les courbes représentant le profil (3) de la peau (2) dans des plans axiaux (Pa) espacés prédéterminés et à différents instants des phases d'inspiration et d'expiration sont extraites, par application d'une technique de segmentation, de représentations tridimensionnelles obtenues par un procédé d'imagerie basé sur la résonnance magnétique ou la tomodensitométrie. 11. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que les courbes représentant le profil (3) de la peau (2) dans des plans axiaux (Pa) espacés prédéterminés et à différents instants des phases d'inspiration et d'expiration sont extraites, par segmentation, d'images du type 2D + t effectuées dans lesdits plans axiaux (Pa) et- 17 - obtenues par un procédé d'imagerie basé sur la résonnance magnétique ou la tomodensitométrie. 12. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que la représentation tridimensionnelle évolutive de la peau (2) au niveau de la zone d'intérêt est exploitée, en relation avec au moins une image tridimensionnelle de l'abdomen, pour réaliser une simulation prédictive en temps réel des positions des organes et/ou viscères de la cavité abdominale durant des phases de respiration libre normale du sujet (1). 13. Procédé selon l'une quelconque des revendications 6 à 8, caractérisé en ce que, lorsque le procédé d'imagerie auquel est soumis le sujet (1) est basé sur la résonance magnétique, le marqueur allongé (5) est constitué par une barre formée d'eau gélifiée, avec une section carrée ou rectangulaire et une longueur suffisante pour couvrir au moins la zone du sujet (1) dont on souhaite suivre les mouvements induits par la respiration dudit sujet.