Image processing method based on the finite element method for directly solving inverse problems in structural mechanics

10007,983

14/915690

September 03, 2014

The present disclosure relates to an image processing method including the following steps: receiving a deformation or displacement image; defining a meshing of the deformation image; assigning a respective pair of unknown variables to each node of the meshing; detecting one (or more) node(s) located at the boundary between two different materials of the body studied; enriching the node(s) detected; calculating a basic matrix for each basic cell of the meshing; assembling the basic matrixes in order to obtain an overall structural matrix; calculating a resilience image of the body based on the overall structural matrix.

UNIVERSITE GRENOBLE ALPES (St. Martin d'Hères, FR)

Universite Grenoble Alpes (St. Martin d'Heres, FR)

1. A method for image processing for formation of an elasticity image of a body, the method comprising the steps of: receiving a deformation image illustrating a displacement or deformation field of the points of the body as a function of a pressure difference in the body, wherein the deformation image is a deformation elastogram produced by introducing an ultrasonic wave into the body; defining a meshing of the deformation image, the meshing step applying a finite element method, the meshing including a plurality of elementary cells each delimiting a same material and each including at least three nodes, each node belonging to one or several adjacent cells to the meshing; assigning to each node i of the meshing a pair of unknown nodal variables, each pair being representative of elastic properties to be determined; detecting at least one node common to at least two adjacent elementary cells delimiting different materials; enriching the at least one detected common node the enrichment step including replacing the pair of unknown variables assigned to the detected common node with at least two pairs of enrichment unknown variables, each pair of enrichment unknown variables of the detected common node being assigned to a respective elementary cell from among the at least two adjacent elementary cells delimiting different materials; calculating an elementary matrix for each elementary cell of the meshing so as to obtain a plurality of elementary matrices, the elementary matrices being calculated by taking into account pairs of unknown variables assigned to the elementary cells, and pairs of unknown variables assigned to the elementary cells; assembling the elementary matrices of the plurality of elementary matrices in order to obtain a structure matrix; calculating the elasticity image of the body from the deformation image and from the structure matrix; and displaying the elasticity image.

2. The method according to claim 1 , further comprising, prior to the step of defining a meshing, a step consisting of segmenting the deformation image in order to obtain a segmented deformation image including a plurality of representative regions of areas of the body assumed to be in different materials.

3. The method according to claim 2 , wherein the step of defining a meshing is applied on the segmented deformation image, the dimensions and positions of the cells of the meshing being determined depending on the positions and dimensions of the regions of the segmented deformation image.

4. The method according to claim 1 , wherein the step of calculating an elasticity image comprises the resolution of the following matrix equation: { R }=([ Q′] T [Q ′]) −1 [Q′] T {F′} wherein: [Q′] represents the reduced structure matrix, and [Q′] T its transposed, {F′} represents the reduced matrix of the field of forces applied to the nodes, and {R} represents the matrix of the field of elasticities of the body at the nodes.

5. The method according to claim 1 , wherein the step of enriching a detected common node comprises the replacement of the pair of unknown variables assigned to the common node detected by n pairs of enrichment unknown variables, n corresponding to the number of different materials at the boundary of which is located the common node.

6. A system for image processing for formation of an elasticity image of a body, the system comprising: a receiver for receiving a deformation image illustrating a displacement field of the points of the body as a function of a pressure difference in the body, wherein the deformation image is a deformation elastogram produced by introducing an ultrasonic wave into the body; a processor configured to: define a meshing of the deformation image by applying a finite element method, the meshing consisting of a plurality of elementary cells each delimiting a same material and each including at least three nodes, each node belonging to one or several cells adjacent to the meshing; assign to each node i of the meshing a pair of unknown nodal variables ((λ i , μ i) or (E 1 , ν i)), each pair being representative of elastic properties to be determined; detect at least one node common to at least two adjacent elementary cells delimiting different materials; enrich the and at least one detected common node, the enrichment including replacing the pair of unknown variables assigned to the detected common node with at least two pairs of unknown enrichment variables, each pairs of unknown enrichment variables of the detected common node being assigned to a respective elementary cell from among the at least two adjacent elementary cells delimiting different materials; calculate an elementary matrix for each elementary cell of the meshing so as to obtain a plurality of elementary matrices, the elementary matrices being calculated while taking into account pairs of unknown variables assigned to the elementary cells, and pairs of unknown variables assigned to the elementary cells; assemble the elementary matrices from the plurality of elementary matrices for obtaining a structure matrix; and calculate the elasticity image of the body from the deformation image and from the structure matrix; and a display device, wherein the processor is further configured to display the elasticity image on the display device.

7. The system according to claim 6 , wherein the processor is configured to segment the deformation image in order to obtain a segmented deformation image consisting of a plurality of representative regions of areas of the body assumed to be in different materials.

8. A computer program product for executing image processing for formation of an elasticity image of a body comprising a non-transitory computer readable medium storing instructions that, when executed by a computer, cause the computer to perform operations comprising: receiving a deformation image illustrating a displacement or deformation field of the points of the body as a function of a pressure difference in the body, wherein the deformation image is a deformation elastogram produced by introducing an ultrasonic wave into the body; defining a meshing of the deformation image, the meshing step applying a finite element method, the meshing including a plurality of elementary cells each delimiting a same material and each including at least three nodes, each node belonging to one or several adjacent cells to the meshing; assigning to each node i of the meshing a pair of unknown nodal variables, each pair being representative of elastic properties to be determined; detecting at least one node common to at least two adjacent elementary cells delimiting different materials; enriching the at least one detected common node, the enrichment step including replacing the pair of unknown variables assigned to the detected common node with at least two pairs of enrichment unknown variables, each pair of enrichment unknown variables of the detected common node being assigned to a respective elementary cell from among the at least two adjacent elementary cells delimiting different materials; calculating an elementary matrix for each elementary cell of the meshing so as to obtain a plurality of elementary matrices, the elementary matrices being calculated by taking into account pairs of unknown variables assigned to the elementary cells, and pairs of unknown variables assigned to the elementary cells; assembling the elementary matrices of the plurality of elementary matrices in order to obtain a structure matrix; and calculating the elasticity image of the body from the deformation image and from the structure matrix; and displaying the elasticity image.

6165128 December 2000 Cespedes et al.
6277074 August 2001 Chaturvedi
7318804 January 2008 Weitzel
7338452 March 2008 Shiina et al.
8394026 March 2013 Eskandari
8660326 February 2014 Ohayon et al.
8965487 February 2015 Bouma
9554777 January 2017 Kim
2004/0009459 January 2004 Anderson
2010/0160778 June 2010 Eskandari et al.
2011/0093243 April 2011 Tawhai
2011/0282182 November 2011 Ohayon

Bouvier, Adeline, et al.; “A Direct Vulnerable Atherosclerotic Plaque Elasticity Reconstruction Method Based on an Original Material-Finite Element Formulation: Theoretical Framework,” Physics in Medicine and Biology, vol. 58, No. 23, Dec. 7, 2013, pp. 8457-8476. cited by applicant
Céspedes, E. I., et al., “Intraluminal Ultrasonic Palpation: Assessment of Local and Cross-Sectional Tissue Stiffness,” Ultrasound in Med. & Biol., vol. 26, No. 3, 2000, pp. 385-396. cited by applicant
Deleaval, Flavien, et al., “The Intravascular Ultrasound Elasticity-Palpography Technique Revisited: A Reliable Tool for the in vivo Detection of Vulnerable Coronary Atherosclerotic Plaques,” Ultrasound in Med. & Biol., vol. 39, No. 8, Mar. 1, 2013, pp. 1469-1481. cited by applicant
Doyley, M. M., et al.; “Evaluation of an Iterative Reconstruction Method for Quantitative Elastography,” Phys. Med. Biol., 45, 2000, pp. 1521-1540. cited by applicant
Doyley, M. M.; “Model-Based Elastography: A Survey of Approaches to the Inverse Elasticity Problem,” Phys. Med. Biol., 57, Jan. 6, 2012, pp. R35-R73. cited by applicant
Le Floch'H, Simon, et al.; “On the Potential of a New IVUS Elasticity Modulus Imaging Appraoch for Detecting Vulernable Atherosclerotic Coronary Plaques: in vitro Vessel Phantom Study,” Phys. Med. Biol., Sep. 8, 2010, pp. 5701-5721. cited by applicant
Le Floc'H, Simon, et al.: “Vulerable Atherosclerotic Plaque Elasticity Reconstruction Based on a Segmentation-Driven Optimization Procedure Using Strain Measurements: Theoretical Framework,” IEEE Transactions on Medical Imagining, vol. 28, No. 7, Jul. 2009, pp. 1126-1137. cited by applicant
Moës, Nicolas, et al.; “A Finite Element Method for Crack Growth Without Remeshing,” International Journal for Numerical Methods in Engineering, vol. 46, Sep. 10, 1999, pp. 131-150. cited by applicant
Oberai, Assad A., et al.; “Solution of Inverse Problems in Elasticity Imaging Using the Adjoint Method,” Inverse Problems, 19, 2003, pp. 297-313. cited by applicant
Zhu, Yanning, et al.; “A Finite-Element Approach for Young's Modulus Reconstruction,” IEEE Transactions on Medical Imaging, vol. 22, No. 7, Jul. 2003, pp. 890-901. cited by applicant

edspgr.10007983