Automatic method of predictive determination of the position of the skin

Abstract

An automatic process of predictive determination of the position and movements of the skin of a subject in a zone of interest, with the subject breathing freely or in an assisted manner, includes preliminarily acquiring multiple configurations of the skin profile in axial planes, at given successive times, in different respiratory positions, and for each axial plane, constructing at least one deformable digital model starting from different skin profiles, then noting, in a repetitive manner, the actual position of a point on the skin at the level of each axial plane, whose position is significantly modified during inhalation and exhalation phases, and providing, essentially in real time, a simulation of the skin profile in each axial plane, as a function of the actual position noted, and an evolving three-dimensional representation of the skin at the level of the zone of interest, by interpolation between the different axial planes.

Claims

1. An automatic process of predictive determination, as a digital simulation, of a position and movements of skin of a human or animal subject in a zone of interest located at a level of a thorax or an abdomen, with said subject being stretched on said subject's back and breathing, said process comprising: in a preliminary phase, acquiring or determining multiple positions and configurations of a profile of the skin in predetermined, consecutive axial planes that are essentially perpendicular to the craniocaudal axis and distributed along the craniocaudal axis at determined locations, at given successive times, in different respiratory positions, during inhalation and exhalation phases of the subject, during a preliminary treatment phase, constructing at least one independent deformable digital model of the profile of the skin in each said axial plane (Pa) from data pertaining to different skin profiles (3) that are acquired or determined previously, noting, in a repetitive manner, an actual position of a point or a specific point zone (4) on the skin (2) of the subject (1) at the level of each of said axial planes (Pa), whose position is modified in a during the inhalation and exhalation phases, and providing, essentially in real time, a simulation (3′) of the profile of the skin in each said axial plane (Pa), as a function of the actual position noted for each point or specific point zone (4) and by exploiting the deformable digital model of the skin profile in the corresponding axial plane (Pa), and, an evolving three-dimensional representation (2′) of the skin (2) at the level of the zone of interest, this on the basis of a three-dimensional model for overall simulation of the skin composed of plural deformable and independent elementary two-dimensional axial models established in the axial planes (Pa) and by interpolation between said axial planes (Pa), wherein the points or specific point zones (4), whose actual position is monitored over time, are distributed along the craniocaudal axis by being aligned or not, with their number being at least equal 10 and their spacing of at most 5 cm in the direction of the craniocaudal axis.

2. Process according to claim 1, wherein the points or specific point zones (4) of the skin (2), whose actual position is monitored over time, are placed in a parallel plane that is close to or merged with the mid-sagittal plane (MSP) of the subject (1).

3. Process according to claim 1, further comprising, for each axial plane (Pa) in question, generating a curve (3′) that represents the profile (3) of the skin (2) in the axial plane (Pa) in question by applying a deformation to a base curve (Cb) of the skin profile, determined in an end inhalation or exhalation state of the subject, by movement of at least points spaced on the base curve, based on the measured position of a point (4) on the skin (2) of the subject (1) located in the axial plane (Pa) in question.

4. Process according to claim 1, wherein each deformable digital model of the skin profile, respectively associated with one of the predetermined axial planes (Pa), is obtained by pairing points in an essentially radial direction, in the axial plane (Pa) in question, between at least two shapes or configurations that are different from the curve of the skin profile, corresponding to at least two different respective states of the respiratory movement of the subject.

5. Process according to claim 1, wherein each deformable digital model of the skin profile, associated with one of the predetermined axial planes (Pa), is obtained by generating a reference curve (Cr), located mid-way between the two skin profile curves corresponding to the end inhalation and exhalation states, to determine the intersection points, with the two skin profile curves, of a normal to the reference curve (Cr) and passing through a point of this curve, for a number of points spaced along the reference curve (Cr), and to assemble in an ordered manner these pairs of intersection points to constitute a deformation field corresponding to the movement of the skin (2) in the axial plane (Pa) in question between the end inhalation and exhalation states, as a function of the respiratory state of the subject (1).

6. Process according to claim 1, further comprising constructing two deformable digital models of the skin profile for each axial plane (Pa), namely a model for the inhalation phases and a model for the exhalation phases, and detecting the respiratory phase in progress during the production of simulated curves (3′) of the skin profile in the different axial planes (Pa) by means of a specific sensor or by evaluating the direction of movement of the points (4).

7. Process according to claim 1, wherein the specific point positions (4) on the skin (2) are noted in real time by an optical sensor device, with or without preliminary installation of optical markers on the patient's skin (2) in the locations of specific points (4), with the latter being selected in such a way as to obtain a modification of their spatial position during the respiratory movement of the subject.

8. Process according to claim 1, wherein the positions of the points (4) of the skin (2) located in the mid-sagittal plane (MSP) of the subject (1) and in the predetermined axial planes (Pa) are determined by at least one marker (5).

9. Process according to claim 8, wherein the specific points (4) on the skin (2) are each highlighted by a marker (5) that is visible in the imagery process that is used, with these specific points (4) being located in such a way that their spatial position is modified during the respiratory movement of the subject.

10. Process according to claim 8, wherein the markers (5) that are used consist of electromagnetic-type markers resting on the skin (2).

11. Process according to claim 8, wherein the one marker (5) that is used comprises a continuously elongated marker, resting directly or indirectly on the skin (2) of the patient (1), conforming in shape to the local surface of the skin (2) and placed in the mid-sagittal plane (MSP).

12. Process according to claim 11, wherein the marker (5) that is used comprises a flexible and elongated body, made of a material that is visible in the imagery system that is used, extending at least over the length of the zone of interest in the direction of the craniocaudal axis (X), exhibiting a cross-section of determined minimal size and shape and located at a determined distance from the skin (2) of the subject (1), by resting on a strip (6) of a material that is flexible in the imagery system that is used.

13. Process according to claim 11, further comprising: determining and noting, in the sagittal plane, the curve (Cm) corresponding to the median axis of the elongated marker (5), formed by the series of different centers of gravity in the cross-section of said marker (5), or in the different axial planes (Pa), the local centers of gravity of said marker (5), and determining the positions of specific points (4) on the skin (2) of the subject (1) whose movements are taken into account, by exploiting the data pertaining, to this curve (Cm) or to these local centers of gravity, to the cross-section and to the cross-sectional dimensions of the marker (5), and to the thickness of the intermediate material strip (6).

14. Process according to claim 11, wherein when the imagery process to which the subject (1) is subjected is based on magnetic resonance, the elongated marker (5) comprises a bar formed by water gel, with a square or rectangular cross-section, and a sufficient length to cover at least the zone of the subject (1) whose movements induced by the breathing of said subject are to be monitored.

15. Process according to claim 1, wherein the curves showing the profile (3) of the skin (2) in spaced axial planes (Pa) that are predetermined and at different times for the inhalation and exhalation phases are extracted, by sampling along these axial planes, from three-dimensional representations of the surface of the skin (2) at the level of the zone of interest, acquired in advance and obtained by highlighting and monitoring the skin by markers attached to the skin or by structured light projected on the skin.

16. Process according to claim 1, wherein the curves showing the profile (3) of the skin (2) in spaced axial planes (Pa) that are predetermined and at different times for the inhalation and exhalation phases are extracted, by application of a segmentation technique, from three-dimensional representations obtained by an imagery process based on magnetic resonance or tomodensitometry.

17. Process according to claim 1, wherein the curves showing the profile (3) of the skin (2) in spaced axial planes (Pa) that are predetermined and at different times for the inhalation and exhalation phases are extracted, by segmentation, from 2D+t-type images made in said axial planes (Pa) and obtained by an imagery process based on magnetic resonance or tomodensitometry.

18. Process according to claim 1, wherein the evolving three-dimensional representation of the skin (2) at the level of the zone of interest is exploited, in connection with at least one three-dimensional image of the abdomen, for carrying out a predictive simulation in real time of the positions of the organs or viscera of the abdominal cavity during normal free breathing phases of the subject (1).

19. An automatic process of predictive determination, as a digital simulation, of a position and movements of skin of a human or animal subject in a zone of interest located at a level of a thorax or an abdomen, with said subject being stretched on said subject's back and breathing, said process comprising: in a preliminary phase, acquiring or determining multiple positions and configurations of a profile of the skin in predetermined, consecutive axial planes that are essentially perpendicular to the craniocaudal axis and distributed along the craniocaudal axis at determined locations, at given successive times, in different respiratory positions, during inhalation and exhalation phases of the subject, during a preliminary treatment phase, constructing at least one independent deformable digital model of the profile of the skin in each said axial plane (Pa) from data pertaining to different skin profiles (3) that are acquired or determined previously, noting, in a repetitive manner, an actual position of a point or a specific point zone (4) on the skin (2) of the subject (1) at the level of each of said axial planes (Pa), whose position is modified in a during the inhalation and exhalation phases, and providing, essentially in real time, a simulation (3′) of the profile of the skin in each said axial plane (Pa), as a function of the actual position noted for each point or specific point zone (4) and by exploiting the deformable digital model of the skin profile in the corresponding axial plane (Pa), and, an evolving three-dimensional representation (2′) of the skin (2) at the level of the zone of interest, this on the basis of a three-dimensional model for overall simulation of the skin composed of plural deformable and independent elementary two-dimensional axial models established in the axial planes (Pa) and by interpolation between said axial planes (Pa), wherein each deformable digital model of the skin profile, respectively associated with one of the predetermined axial planes (Pa), is obtained by pairing points in an essentially radial direction, in the axial plane (Pa) in question, between at least two shapes or configurations that are different from the curve of the skin profile, corresponding to at least two different respective states of the respiratory movement of the subject.

20. An automatic process of predictive determination, as a digital simulation, of a position and movements of skin of a human or animal subject in a zone of interest located at a level of a thorax or an abdomen, with said subject being stretched on said subject's back and breathing, said process comprising: in a preliminary phase, acquiring or determining multiple positions and configurations of a profile of the skin in predetermined, consecutive axial planes that are essentially perpendicular to the craniocaudal axis and distributed along the craniocaudal axis at determined locations, at given successive times, in different respiratory positions, during inhalation and exhalation phases of the subject, during a preliminary treatment phase, constructing at least one independent deformable digital model of the profile of the skin in each said axial plane (Pa) from data pertaining to different skin profiles (3) that are acquired or determined previously, noting, in a repetitive manner, an actual position of a point or a specific point zone (4) on the skin (2) of the subject (1) at the level of each of said axial planes (Pa), whose position is modified in a during the inhalation and exhalation phases, and providing, essentially in real time, a simulation (3′) of the profile of the skin in each said axial plane (Pa), as a function of the actual position noted for each point or specific point zone (4) and by exploiting the deformable digital model of the skin profile in the corresponding axial plane (Pa), and, an evolving three-dimensional representation (2′) of the skin (2) at the level of the zone of interest, this on the basis of a three-dimensional model for overall simulation of the skin composed of plural deformable and independent elementary two-dimensional axial models established in the axial planes (Pa) and by interpolation between said axial planes (Pa), wherein each deformable digital model of the skin profile, associated with one of the predetermined axial planes (Pa), is obtained by generating a reference curve (Cr), located mid-way between the two skin profile curves corresponding to the end inhalation and exhalation states, to determine the intersection points, with the two skin profile curves, of a normal to the reference curve (Cr) and passing through a point of the reference curve (Cr), for a number of points spaced along the reference curve (Cr), and to assemble in an ordered manner these pairs of intersection points to constitute a deformation field corresponding to the movement of the skin (2) in the axial plane (Pa) in question between the end inhalation and exhalation end states, as a function of the respiratory state of the subject (1).

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.

(2) The invention will be better understood using the description below, which relates to two embodiments, provided by way of nonlimiting examples and explained with reference to the accompanying schematic drawings, in which:

(3) FIGS. 1A to 1C are views, obtained by MRI, according to the mid-sagittal plane of a human subject provided with a wire marker according to the invention, respectively in the exhaled state (FIG. 1A), in the inhaled state (FIG. 1B), and by superposition of the different positions taken by the marker (by means of its median curve) in the inhalation and exhalation phases (FIG. 1C);

(4) FIGS. 2A to 2D are views, obtained by MRI, along an axial plane at the level of a zone of interest of a human subject provided with the marker of FIG. 1, illustrating the profile of the skin that is determined in the exhaled state (FIG. 2A) and the inhaled state (FIG. 2B), with the calculation of the reference curve from two curves of profiles extracted from FIGS. 2A and 2B (FIG. 2C), and the generation of a deformation field from pairs of points (pairing of points) belonging to two skin profile curves of FIGS. 2A and 2B (FIG. 2D);

(5) FIG. 3 is a view that is similar to that of FIG. 1B with an illustration of the median curve of the marker;

(6) FIGS. 4A to 4C are detail axial cutaway views of the marker illustrating the stages for determination of the contours and the slope of the marker;

(7) FIGS. 5A to 5D are axial cutaway views illustrating the stages of segmentation and extraction of the skin profile;

(8) FIG. 6 is a detail axial cutaway view of the marker and its immediate environment, obtained by superposition and illustrating different positions of the center of gravity of the marker and the hysteresis between the positions taken during the inhalation and those taken during the exhalation;

(9) FIGS. 7A and 7B are schematic representations, respectively from the top and side, of the torso and the head of the human subject equipped with a marker according to the invention for monitoring the respiratory movements, and

(10) FIGS. 8A to 8C are representations of a three-dimensional digital simulation of the zone of interest of a human patient, controlled with breathing, obtained and used within the framework of the invention, with FIGS. 8B and 8C simultaneously showing views acquired by the imagery in particular planes.

DETAILED DESCRIPTION OF THE INVENTION

(11) The invention relates to an automatic process for predictive determination, for example in the form of a digital simulation, of the position and movements of the skin 2 of a human subject 1 or an animal in a zone of interest located at the level of the thorax and/or the abdomen.

(12) Within the framework of a preferred embodiment of the invention, it is considered that the subject 1 is stretched on his/its back and breathes freely or in an assisted manner.

(13) In agreement with the general principles and the essential characteristics of the invention, the process consists, in a preliminary manner, in acquiring or determining multiple positions and configurations of the skin 2 in at least two positions of the respiratory cycle, in such a way as to be able to observe a respiratory movement whose amplitude is essentially similar to the amplitude corresponding to the breathing. These acquisitions are made using a conventional 3D imager of the CT-scan type. In addition, these acquisitions are to cover the zone of interest of which it is subsequently desired to be able to simulate the movement.

(14) Once these 3D images have been acquired, the next step is to extract from them the position of the skin, which, for example, by applying a simple operator for processing the thresholding-type image, is common practice and can be carried out easily by one skilled in the art. Once the position of the skin has been extracted and calculated in each 3D image, the intersections of this non-planar surface with a certain number of axial planes Pa selected consecutively along the craniocaudal axis of the patient are calculated.

(15) It is possible to use only two 3D acquisitions, i.e., corresponding, for example, respectively to the exhaled position and an inhaled position of the patient, and a series of curve pairs (representing the skin profiles), corresponding respectively to the skin profile of the patient in exhalation and in inhalation in each axial section Pa, is then obtained.

(16) The invention then provides for constructing an independent digital model for each axial section, corresponding, by way of simple example, to a linear interpolation between the position of the exhalation curve and the position of the inhalation curve.

(17) Thus, by monitoring the position in a repetitive manner from a point or a specific point zone 4 of the skin 2 of the subject 1 at the level of each of the above-mentioned axial planes Pa, and of which the position is modified in a significant manner during the respiratory movement, it is possible to provide, essentially in real time, a simulation 3′ of the skin profile in each axial plane Pa by exploiting the previously-described deformable digital model. The invention then makes it possible to create an evolving three-dimensional representation 2′ of the skin 2 at the level of the zone of interest by interpolation between the simulated position of the skin at the level of the different axial planes Pa.

(18) One possible effective way of monitoring a point for each modeled axial section of the patient's skin (skin profile) may be to choose to position all of the axial points to be monitored on the same plane, which can be the mid-sagittal plane. In this case, it is necessary only to monitor the movement of the profile of the skin at the level of the mid-sagittal plane to know the position of all of the points that must be monitored in each axial plane, since these points constitute the profile of the skin at the level of the mid-sagittal plane.

(19) The invention is described below in a more detailed manner, in connection with a particular embodiment, illustrated by the accompanying figures and having variants for certain aspects.

(20) As above, it is considered that the subject 1 is stretched on his/its back and breathes freely or in an assisted manner.

(21) The process consists, in a preliminary manner, in acquiring or determining multiple positions and configurations of the profile 3 of the skin 2 in predetermined consecutive axial planes Pa, essentially perpendicular to the craniocaudal axis X and distributed along this axis at determined locations, at given successive times, in at least two positions of the respiratory cycle, in such a way as to be able to observe a respiratory movement whose amplitude is essentially similar to the amplitude corresponding to normal breathing. The determination of the position of the skin at the level of the axial profiles can be carried out simultaneously for multiple profiles at the same time, if the acquisition device allows it.

(22) In accordance with the invention, this process also consists, during a preliminary phase, for each axial plane Pa in question, in constructing at least one deformable digital model of the skin (profile) starting from data pertaining to the different skin profiles 3 that are acquired or determined previously. Next, in a repetitive manner, the actual current position of a point or a specific point zone 4 of the skin 2 of the subject 1 is noted at the level of each of the above-mentioned axial planes Pa (preferably highlighted using a marker 5 that is visible in the imagery process that is used), whose position is modified in a significant manner by the respiratory movement. Finally, there is provided, essentially in real time, a simulation 3′ of the profile of the skin in each axial plane Pa, based on the actual position noted and by exploiting the corresponding deformable digital model, but also optionally an evolving three-dimensional representation 2′ of the skin 2 at the level of the zone of interest, by interpolation between the different axial planes Pa.

(23) In an advantageous manner, and so as to facilitate the monitoring of the movement of the skin and to increase its precision, it is desirable that the points or specific point zones 4 of the skin 2, whose actual position is monitored over time, be placed in a parallel plane that is close to or merged with the mid-sagittal plane MSP of the subject 1, by preferably being spaced mutually by a distance on the order of one centimeter.

(24) However, it may also be provided that the points or specific point zones 4, whose actual position is monitored over time, are distributed along the craniocaudal axis by being aligned or not, with their number advantageously being at least equal to approximately 10 and their spacing of at most 5 cm in the direction of the above-mentioned axis.

(25) In accordance with the invention, the process preferably consists in, for each axial plane Pa in question, generating essentially in real time a curve 3′, which shows the profile 3 of the skin 2 in the axial plane Pa in question, by applying a deformation to a base curve Cb of the axial profile of the skin, for example by moving at least certain points spaced on the base curve, based on the measured position of a point 4 on the skin 2 of the subject 1 located in the axial plane Pa in question.

(26) The base curve that is used preferably corresponds to the axial profile of the skin noted in the maximum exhalation state of the subject 2.

(27) Advantageously, and as FIG. 2D illustrates by way of example (for inhalation and exhalation positions), each deformable digital model of a skin profile, respectively associated with one of the predetermined axial planes Pa, is obtained by pairing points in an essentially radial direction, in the axial plane Pa in question, between at least two shapes or configurations that are significantly different from the curve of the skin profile, corresponding to at least two different respective states of the respiratory movement of the subject.

(28) In accordance with another characteristic of the invention, specifying the preceding embodiment and illustrated by FIGS. 2A to 2D, each deformable digital model of the skin profile, associated with one of the predetermined axial planes Pa, is obtained by generating a reference curve Cr located mid-way between the two skin profile curves corresponding to the end inhalation and exhalation states to determine the intersection points, with the two above-mentioned curves, of a normal to the reference curve Cr and passing through a point of this curve, for a number of points spaced along the latter, and to assemble in an ordered manner these pairs of intersection points to constitute a deformation field corresponding to the movement of the skin 2 in the axial plane Pa in question between the two above-mentioned end states, as a function of the respiratory state of the subject 1.

(29) The above-mentioned points of the curve Cr are spaced enough so that the normals do not cross before intersecting the two axial profiles of the skin.

(30) As indicated previously, the three-dimensional model of overall simulation of the skin consists of a number of elementary axial models that are two-dimensional and independent of the skin, established in consecutive axial planes PA along the craniocaudal axis.

(31) In practice, once the skin is segmented in the images corresponding respectively to the end inhalation and exhalation positions of the subject (FIGS. 2A and 2B), each axial model can be established by executing the following four operating stages (see FIGS. 2C and 2D):

(32) 1) Creation of a reference curve Cr located exactly mid-way between the two curves corresponding to the two above-mentioned end positions. This stage can be carried out by using a distance chart corresponding to the difference between the image mask relative to the position of the skin in the inhaled state and the image mask relative to the position of the skin in the exhaled state.

(33) 2) Resampling the axial reference curve Cr in such a way as to obtain a set of consecutive points that are spaced in a regular manner.

(34) 3) Generation of a normal to the reference curve Cr at the level of each point of this resampled curve.

(35) 4) Determination of two intersection points corresponding to the intersection of this normal with the curves of the inhaled and exhaled states that frame it. A set of vectors that are normal to the curve Cr and that thus represent a deformation field reflecting the movement of the surface of the skin in a given axial plane Pa between the two original positions of the different respiratory states is then obtained.

(36) Thus, by repeating this operation for all of the consecutive axial sections, an axial model of the respiratory movement is obtained all along the craniocaudal axis; each set of data resulting from a 2D acquisition makes possible the construction of a local axial model of the patient's skin. Taking into account the variation of the shape of the skin of a conventional subject, the generation of axial models of the skin that are separated from one another by approximately 1 cm makes it possible to obtain by interpolation a reliable and precise three-dimensional model of the patient's skin at the level of the zone of interest.

(37) So as to take into account the hysteresis phenomenon of the movement of the skin during breathing, it is provided within the framework of the process according to the invention to construct two deformable digital models of axial profiles of the skin in each axial plane Pa, namely a model for the inhalation phases and a model for the exhalation phases, and to detect the type of respiratory phase in progress during the production of simulated curves 3′ of a skin profile in the different axial planes Pa by means of a specific sensor or by evaluating the direction of movement of the points 4 by considering individually their variations of position relative to their preceding positions.

(38) According to a first variant embodiment, the positions of the specific points 4 on the skin 2 are noted in real time by means of an optical sensor device, in particular with a laser, with or without preliminary installation of optical markers on the patient's skin 2 in the locations of specific points 4, with the latter advantageously being selected in such a way as to obtain a significant modification of their spatial position during the respiratory movement of the subject.

(39) To facilitate the monitoring of the sagittal profile of the skin that corresponds to the concatenation of the different axial points, points 4 on the skin 2 whose movements are taken into account within the framework of the process as initial parameters, the positions of these points 4, located in the mid-sagittal plane MSP of the subject 1 and in the predetermined axial planes Pa, are determined by means of at least one marker 5, preferably a continuously elongated marker, resting directly or indirectly on the skin 2 of the patient 1, conforming in shape to the local surface of the skin 2 and placed in the mid-sagittal plane MSP.

(40) Alternatively, the markers 5 that are used consist of electromagnetic-type markers, such as, for example, miniature coils, resting directly or indirectly on the skin 2, if necessary with a known spacing relative to the latter.

(41) Thus, the specific points 4 on the skin 2 are each highlighted by means of a marker 5 that is particularly readily visible in the imagery process that is used, with these specific points 4 being advantageously located in such a way that their spatial position is modified significantly during the respiratory movement of the subject.

(42) As FIGS. 1 to 4, 6 and 7 show, the marker 5 that is used consists of a flexible and elongated body, made of a material that is particularly readily visible in the imagery system that is used, extending at least over the length of the zone of interest in the direction of the craniocaudal axis X, exhibiting a cross-section of determined minimal size and shape and located at a determined distance from the skin 2 of the subject 1, by resting on a strip 6 of a material that is flexible and essentially invisible in the imagery system that is used.

(43) Knowing the dimensions of the marker 5 and those of the material strip 6 and for the purpose of facilitating the reliable highlighting of the position of the marker 5, the process can advantageously consist in determining and noting, in the sagittal plane, the curve Cm corresponding to the median axis of the elongated marker 5, for example formed by the series of different centers of gravity in the cross-section of said marker 5, and/or in the different axial planes Pa in question, the local centers of gravity of said marker 5, and in determining the positions of specific points 4 on the skin 2 of the subject 1 whose movements are taken into account, by exploiting the data pertaining, on the one hand, to this curve Cm or to these local centers of gravity, on the other hand, to the cross-section and to the cross-sectional dimensions of the marker 5, and finally, to the thickness of the intermediate material strip 6.

(44) In accordance with a preferred practical implementation of the invention, as shown in the accompanying figures, it may be provided that, when the imagery process to which the subject 1 is subjected is based on magnetic resonance, the elongated marker 5 consists of a bar formed by water gel, with a square or rectangular cross-section and a sufficient length to cover at least the zone of the subject whose movements induced by the breathing of said subject 1 are to be monitored (zone of interest).

(45) In this context, the segmentation of the skin in sagittal images obtained by magnetic resonance is advantageously carried out by exploiting the very good visibility of the marker 5 in the MRI images and its contrast relative to the textile strip 6 by means of which it rests on the skin of the subject (the thicknesses of said marker and said strip being known data).

(46) More precisely, this sagittal segmentation phase of the skin can be broken down into four primary stages, namely:

(47) 1) A 2D working image that is oversampled and created from the original 2D image (taken on the subject) by implementing a cubic interpolation (this stage increases the precision of the determination of the position of the marker 5). In practice, and taking into account geometric deformations and distortions that affect the 2D MRI images (at the level of the two opposite sides of the ends of the explored zone—see, for example, FIG. 3), only the central part of the initial images is exploited (approximately 60%).

(48) 2) Taking into account differences in contrast between consecutive 2D+t images, it is advantageous to determine a threshold value corresponding to the position of the marker 5. Then, in each vertical straight line of the different images, the first string of consecutive pixels having a length that is at least greater than approximately 80% of the thickness of the marker 5 is sought, and the center of the marker 5 is then calculated along each vertical straight line (making it possible to determine the position of the marker along this straight line). The best threshold value corresponds to the one that minimizes the mean distance between the simulated position of the marker and the position extracted, by segmentation, from the original images.

(49) 3) The parts of the elongated marker 5 that are not continuous are then located, and the trace of the elongated marker is modified locally by extrapolating the uniform portions of the marker located on both sides of the discontinuous part.

(50) 4) An operation for smoothing positions of the consecutive centers of gravity of the marker is carried out in such a way as to correct the segmentation errors due to local information gaps (the use of the center of gravity makes it possible to limit the influence of the effect of partial volume).

(51) Likewise, the segmentation of the marker 5 can be carried out in the axial planes Pa that are provided for the acquisition of 2D+t images by executing the following four operating stages:

(52) 1) The search for the location of the marker 5 with a square cross-section, for example, is limited to a search zone defined manually at the beginning of the procedure and then centered automatically by being based on the preceding position of the center of the marker.

(53) 2) The best threshold value that makes it possible to obtain at least 80% of the theoretical surface of the cross-section of the marker is sought (pixels whose gray value is less than the threshold value are sought).

(54) 3) As soon as the mask showing the cross-section of the marker is found, its center of gravity is calculated by using a known algorithm for calculating the barycenter.

(55) 4) The slope of the marker is then evaluated starting from the preceding mask by carrying out, for example, the following operations: determination of the four disk portions of the larger disk centered on the center of gravity of the marker (FIG. 4B); determination of the centers of gravity of these four disk portions (CG1, CG2, CG3, CG4); determination of the horizontal and vertical slopes of the marker by approximation starting from segments connecting, on the one hand, the points CG1 and CG3 and, on the other hand, the points CG2 and CG4.

(56) Finally, the segmentation of the skin in the 2D+t axial images is carried out by using the data pertaining to the position of the marker 5 provided previously.

(57) This axial segmentation can be carried out by executing the following four operating stages (see FIGS. 5A to 5D):

(58) 1) A binary mask corresponding to the threshold of the original image is calculated.

(59) 2) Only the largest connected elements of the mask are preserved, and the mask is filled (FIGS. 5B and 5C).

(60) 3) The curve corresponding to the peripheral limit of the mask is extracted. Then, this curve is smoothed and completed to remedy the errors and information gaps in the original image.

(61) 4) The curve corresponding to the position of the skin is corrected close to the marker 5, knowing the distance between the center of the marker and the skin.

(62) As already mentioned previously, and so as to take into account the hysteresis affecting the skin between the inhalation and exhalation phases, it is recommended according to the invention to construct a different simulation model of the skin as a function of the respiratory cycle phase (inhalation and exhalation).

(63) For this purpose, the curves illustrating the positions of the marker 5 should be classified in two categories: one of which relates to the inhalation phase and the other to the exhalation phase.

(64) The hysteresis effect is maximal at the end of the inhalation phase and at the beginning of the exhalation phase, and the discrimination of the curves is trivial except for the end inhalation state (see FIG. 6).

(65) So as to be able to distinguish the positions of the marker in this critical state, it is advantageously possible to exploit simultaneously the position of the center of the marker and the curvature of the skin around the marker (the inventors noted that the lateral position of the marker changed because of the modification of the curvature of the skin).

(66) Within the framework of a preliminary acquisition and processing phase of the process according to the invention, the curves showing the profile 3 of the skin 2 in spaced axial planes Pa that are predetermined and at different times of the inhalation and exhalation phases can be obtained in different ways.

(67) Thus, these curves can be: extracted, by sampling along these axial planes, from three-dimensional representations of the surface of the skin 2 at the level of the zone of interest, acquired in advance and obtained by highlighting and monitoring the skin by means of markers (electromagnetic, optical, or others) attached to the skin or by means of structured light projected on the skin; extracted, by application of a segmentation technique, from three-dimensional representations obtained by an imagery process based on magnetic resonance or tomodensitometry; extracted, by segmentation, from 2D+t-type images made in said axial planes Pa and obtained by an imagery process based on magnetic resonance or tomodensitometry.

(68) Within the framework of an advantageous medical application of the results provided by the invention, the evolving three-dimensional representation 2′ of the skin 2 at the level of the zone of interest is exploited in connection with at least one three-dimensional image of the abdomen, for carrying out a predictive simulation in real time of the positions of the organs and/or viscera of the abdominal cavity during normal free-breathing phases of the subject 1, by implementing, for example, the teaching of the publication cited in the introductory part of this description.

(69) Thus, thanks to the invention, it is possible to make available a three-dimensional simulation evolving in real time of the skin of a patient by continuously noting only the positions of a limited number of points on the skin, i.e., by limiting the necessary treatment resources and limiting the necessary zone that is highlighted while providing sufficient precision for medical applications.

(70) A practical implementation of the invention made it possible to reach a precision of 1 mm and a refresh rate of the simulation at a frequency of 80 Hz.

(71) In addition, by comparative evaluations, the inventors have been able to estimate that the simulation produced by the invention exhibited a precision at least equivalent to that of a reconstruction of the skin starting directly from images acquired by the imagery system.

(72) The invention also relates to an imagery system, in particular of the MRI type, integrating hardware and software assets allowing the implementation of the process according to the invention, with these means originating in an obvious manner from the description of the process above. This system can, for example, comprise a magnetic resonance imaging device of the MAGNETOM Aera 1.5 Tesla type of the SIEMENS Company.

(73) Of course, the invention is not limited to the embodiments described and shown in the accompanying drawings. Modifications remain possible, in particular from the standpoint of the composition of the various elements or by substitution of equivalent techniques, without thereby exceeding the field of protection of the invention.