System and method for correcting data for deformations during image-guided surgical procedures

Patent metadata

Field	Value
Publication number	US-9767573-B2
Application number	US-201113574624-A
Country	US
Kind code	B2
Filing date	Jan 21, 2011
Priority date	Jan 22, 2010
Publication date	Sep 19, 2017
Grant date	Sep 19, 2017

Abstract

Systems and methods for collecting and processing physical space data for use while performing an image-guided surgical (IGS) procedure are provided. The system and method includes obtaining a computer model of a non-rigid structure of interest in a patient and performing a rigid alignment of the computer model and surface data in a patient space associated with at least a portion of the non-rigid structure. The system and method also include computing a deformation of the computer model that provides a non-rigid alignment of the computer model and surface data, the deformation computed using a set of boundary conditions defined for each node of the computer model based on the rigid alignment and a kernel function. Additionally, the system and method can include displaying data for facilitating the IGS procedure based on the deformation.

First claim

What is claimed is:

1. A method for collecting and processing physical space data for use while performing an image-guided surgical (IGS) procedure, the method comprising:
obtaining a preoperative image of a surface of a non-rigid structure of interest in a patient, wherein the preoperative image was acquired during a preoperative presentation of the structure;
during an intraoperative presentation of the structure, generating in patient space, via digitization equipment of an IGS system, partial surface data extending over a portion of the surface, wherein the portion, as presented intraoperatively, is less than 50% of the surface and wherein during the intraoperative presentation the surface is deformed relative to the preoperative presentation;
performing, via a processor of the IGS system, a rigid alignment, in the patient space, of the partial surface data to a computer model corresponding to the preoperative image, the computer model comprising a plurality of nodes, the plurality of nodes extending over all of the surface of the non-rigid structure;
computing, via the processor of the IGS system, a deformation of the computer model, wherein the deformation provides a non-rigid alignment, in the patient space, of said computer model and the partial surface data, said deformation computed using boundary conditions calculated for each of the plurality of nodes of said computer model, said boundary conditions calculated using a kernel function and spatial difference values for each of the plurality of nodes derived from said rigid alignment;
applying said deformation of the computer model to the preoperative image and thereby computing a modified image; and
displaying, at a display of the IGS system, the modified image.

2. The method of claim 1 , wherein said computing further comprises:
calculating the spatial difference values between the partial surface data and corresponding portions of the computer model using a correspondence function;
selecting a zero value for the spatial difference values outside the corresponding portions of the computer model;
assigning weights to the spatial difference values for each of the plurality of nodes using the kernel function;
computing the boundary conditions using the spatial difference values and said assigned weights, the boundary conditions for each one of the plurality of nodes computed as a weighted average of the one of the plurality of nodes and surrounding ones of the plurality of nodes;
generating displacement field vector values for each of the plurality of nodes of said computer model using the boundary conditions;
producing an updated computer model using the displacement field vector values;
updating the correspondence function using the updated computer model;
recalculating the spatial difference values between the partial surface data and corresponding portions of the computer model using the updated correspondence function;
repeating the calculating, selecting, assigning, computing, generating, producing, updating, and recalculating steps if the set of spatial difference values fail to meet convergence criteria; and
outputting the updated computer model.

3. The method of claim 2 , wherein said generating said displacement field vector values further comprises accumulating said displacement field vector values from each time the generating step is performed, and wherein said outputting further comprises outputting an accumulated displacement vector field comprising a combination of said accumulated displacement field vector values.

4. The method of claim 2 , further comprising selecting said convergence criteria to be at least one of a minimum threshold value for the set of spatial difference values or a minimum difference value between a current one of the set of spatial difference values and a previous one of the set of spatial difference values, wherein said selecting is performed before said updating, after said updating, or after a pre-defined number of iterations.

5. The method of claim 1 , wherein said computing further comprises:
receiving a location of an object or tool in said patient space;
computing reverse deformation displacement vector field values for each of the plurality of nodes of said computer model corresponding to the deformation;
calculating an envelope of additional reverse deformation vector field values for a portion of said patient space surrounding said deformed computed model;
transforming said location into a computer model space of said computer model using said rigid alignment and said reverse deformation displacement vector field values to yield a transformed location; and
calculating coordinates of the object or tool in an image space of said preoperative image from said transformed location of said object or tool.

6. The method of claim 5 , further comprising:
displaying said preoperative image; and
displaying indicia of said object location in said image space of the preoperative image corresponding to said calculated coordinates in said image space.

Assignees

• Miga Michael I
• Dumpuri Prashanth
• Univ Vanderbilt

Inventors

• Miga Michael I
• Dumpuri Prashanth

Classifications

• A61B5/00

  Measuring for diagnostic purposes (radiation diagnosis A61B6/00; diagnosis by ultrasonic, sonic or infrasonic waves A61B8/00); Identification of persons · CPC title

• G06T7/33 Primary

  using feature-based methods · CPC title

• G06G7/60

  for living beings, e.g. their nervous systems {; for problems in the medical field} · CPC title

• G06T17/00 Primary

  Three-dimensional [3D] modelling for computer graphics · CPC title