IP Library Granted Patent US 12,223,666
Granted Patent B2
US 12,223,666 · App. 18/448,289 · Granted Feb 11, 2025

Methods for optimizing the planning and placement of probes in the brain via multimodal 3D analyses of cerebral anatomy

Inventors: Nitin Tandon (Houston, TX); Mehmet Kadipasaoglu (Houston, TX); Kevin Pham (Sugar Lang, TX); Cristian Donos (Bucharest, RO); Kiefer Forseth (Houston, TX); Patrick Sarahan Rollo (Houston, TX)
Assignee: THE BOARD OF REGENTS OF THE UNIVERSITY OF TEXAS SYSTEM
G06T7/37A61B34/10G01R33/5607G06T7/11G06T7/38A61B2034/102A61B2034/107G06T2207/10081G06T2207/10088G06T2207/10104G06T2207/30016
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,223,666
App. No.
18/448,289
Granted
Feb 11, 2025
Kind
B2
Abstract

A method includes obtaining a first imaging scan and a second imaging scan of a single subject brain. The first imaging scan is converted to a first dataset, and the second imaging scan is converted to a second dataset. A sequence-adaptive multimodal segmentation algorithm is applied to the first dataset and the second dataset. The sequence-adaptive multimodal segmentation algorithm performs automatic intensity-based tissue classification to generate a first labelled dataset and a second labeled dataset. The first labeled dataset and the second labeled dataset are automatically co-registered to each other to generate a transformation matrix based on the first labeled dataset and the second labeled dataset. The transformation matrix is applied to align the first dataset and the second dataset.

Claims (109)

1. A method, comprising:

obtaining a first imaging scan and a second imaging scan of a single subject brain;

converting the first imaging scan to a first dataset, and the second imaging scan to a second dataset;

applying an adaptive segmentation algorithm to the first dataset and the second dataset, wherein the adaptive segmentation algorithm performs automatic intensity-based tissue classification to generate a first labeled dataset and a second labeled dataset;

automatically co-registering the first labeled dataset and the second labeled dataset to each other to generate a transformation matrix based on the first labeled dataset and the second labeled dataset; and

applying the transformation matrix to align the first dataset and the second dataset.

2. The method of claim 1 , wherein:

the single subject brain comprises a cerebrospinal fluid region;

the first imaging scan is a contrast weighted scan and the first dataset is a contrast weighted dataset; and

the method further comprises:

selecting voxels of the first dataset identified as belonging to the cerebrospinal fluid region;

applying a filtering algorithm to identify voxels of the first dataset representing blood vessels and assign a vesselness weight value to each voxel;

integrating the vesselness weight values into the first dataset; and

after aligning the first dataset and the second dataset, converting the first dataset to a surface anatomical mesh model.

3. The method of claim 1 , wherein the adaptive segmentation algorithm assigns a numeric label value to each voxel of the first dataset or the second dataset.

4. The method of claim 1 wherein the method further comprises:

defining predicted target point coordinates and entry point coordinates for a probe based on target point coordinates and entry point coordinates of previously implanted probes or by user defined target and entry points;

defining a trajectory for the probe based on a mean target coordinates and mean entry point coordinates;

adjusting the trajectory to intersect with a nearest voxel assigned a label of an anatomical region of interest;

checking proximity of the trajectory to critical structures based on user defined constraints and/or user defined modification of the trajectory to satisfy the user defined constraints; and

superimposing the trajectory on the second data set to form a planning dataset.

5. The method of claim 1 , wherein:

the first imaging scan is an anatomical scan, the second imaging scan is a post-implantation imaging scan, the first dataset is an anatomical dataset, and the second dataset is a post-implantation imaging dataset; and

the method further comprises:

obtaining a third imaging scan used to guide electrode implantation during surgery;

converting the third imaging scan to a third dataset;

aligning a third dataset with the first dataset;

obtaining a trajectory implant data file created during the electrode implantation;

generating a planned trajectory dataset, based on trajectory implant data file, that includes dummy objects disposed at locations of electrode geometry;

aligning the planned trajectory dataset to the imaging dataset; and

automatically identifying and labelling electrodes in an electrode dataset based on the dummy objects of the trajectory implant data file.

6. A non-transitory computer-readable medium encoded with instructions that are executable by one or more processors to:

obtain a first imaging scan and a second imaging scan of a single subject brain;

convert the first imaging scan to a first dataset, and the second imaging scan to a second dataset;

apply an adaptive segmentation algorithm to the first dataset and the second dataset, wherein the adaptive segmentation algorithm performs automatic intensity-based tissue classification to generate a first labelled dataset and a second labeled dataset;

automatically co-register the first labeled dataset and the second labeled dataset to each other to generate a transformation matrix based on the first labeled dataset and the second labeled dataset; and

apply the transformation matrix to align the first dataset and the second dataset.

7. The non-transitory computer-readable medium of claim 6 , wherein:

the single subject brain comprises a cerebrospinal fluid region;

the first imaging scan is a contrast weighted scan and the first dataset is a contrast weighted dataset; and

the instructions are executable by the one or more processors to:

select, based on the labeled dataset, voxels of the first dataset identified as belonging to the cerebrospinal fluid region;

apply a filtering algorithm to identify voxels of the first dataset representing blood vessels and assign a vesselness weight value to each voxel;

integrate the vesselness weight values into the first dataset; and

after aligning the first dataset and the second dataset, convert the first dataset to a surface anatomical mesh model.

8. The non-transitory computer-readable medium of claim 6 , wherein the adaptive segmentation algorithm assigns a numeric label value to each voxel of the first dataset or the second dataset.

9. The non-transitory computer-readable medium of claim 6 , wherein:

the first imaging scan is a contrast weighted scan and the second imaging scan is an anatomical scan; and

the instructions are executable by the one or more processors to:

define predicted target point coordinates and entry point coordinates for a probe based on target point coordinates and entry point coordinates of previously implanted probes or by user defined target and entry points;

define a trajectory for the probe based on a mean target coordinates and mean entry point coordinates;

adjust the trajectory to intersect with a nearest voxel assigned a label of an anatomical region of interest;

check proximity of the trajectory to critical structures based on user defined constraints and/or user defined modification of the trajectory to satisfy the user defined constraints; and

superimpose the trajectory on the second data set to form a planning dataset.

10. The non-transitory computer-readable medium of claim 6 , wherein:

the first imaging scan is an anatomical scan, the second imaging scan is a post-implantation imaging scan, the first dataset is an anatomical dataset, and the second dataset is a post-implantation dataset; and

the instructions are executable by the one or more processors to:

obtain a third imaging scan used to guide electrode implantation during surgery;

convert the third imaging scan to a third dataset;

align a third dataset with the first dataset;

obtain a trajectory implant data file created during the electrode implantation;

generate a planned trajectory dataset, based on trajectory implant data file, that includes dummy objects disposed at locations of electrode geometry;

align the planned trajectory dataset to the post-implantation dataset; and

automatically identify and label electrodes in an electrode dataset based on the dummy objects of the trajectory implant data file.

11. A system, comprising:

one or more processors; and

a memory coupled to the one or more processors, wherein the memory stores instructions that configure the one or more processors to:

obtain a first imaging scan and a second imaging scan of a subject brain;

convert the first imaging scan to a first dataset, and the second imaging scan to a second dataset;

apply an adaptive segmentation algorithm to the first dataset and the second dataset, wherein the adaptive segmentation algorithm performs automatic intensity-based tissue classification to generate a first labeled dataset and a second labeled dataset;

automatically co-register the first labeled dataset and the second labeled dataset to each other to generate a transformation matrix based on the first labeled dataset and the second labeled dataset; and

apply the transformation matrix to align the first dataset and the second dataset.

12. The system of claim 11 , wherein:

the single subject brain comprises a cerebrospinal fluid region;

the first imaging scan is a contrast weighted scan and the first dataset is a contrast weighted dataset; and

the instructions configure the one or more processors to:

select, based on the labeled dataset, voxels of the first dataset identified as belonging to the cerebrospinal fluid region;

apply a filtering algorithm to identify voxels of the first dataset representing blood vessels and assign a vesselness weight value to each voxel;

integrate the vesselness weight values into the first dataset; and

after aligning the first dataset and the second dataset, convert the first dataset to a surface anatomical mesh model.

13. The system of claim 11 , wherein the adaptive segmentation algorithm assigns a numeric label value to each voxel of the first dataset or the second dataset.

14. The system of claim 11 , wherein the instructions configure the one or more processors to:

extract voxels from the first dataset having a label corresponding to a subcortical region of interest;

form a third dataset containing the voxels extracted from the first dataset;

convert the third dataset into a first subcortical surface mesh model;

compute curvature and sulcal features of the first subcortical surface mesh model;

align the first subcortical surface mesh model to a subcortical atlas of the region of interest using the curvature and sulcal features; and

overlay the first subcortical surface mesh model aligned to an atlas of the subcortical region of interest on a second subcortical surface mesh model, the second subcortical surface mesh model having a standardized number of nodes that enables a one-to-one correspondence between node identity and atlas location; and

assign coordinates of nodes of the first subcortical surface mesh model to the second subcortical structure surface mesh model such that the second subcortical surface mesh model assumes a topology of the first subcortical structure surface mesh model.

15. The system of claim 11 , wherein:

the first imaging scan is a contrast weighted scan and the second imaging scan is an anatomical scan; and

the instructions configure the one or more processors to:

define predicted target point coordinates and entry point coordinates for a probe based on target point coordinates and entry point coordinates of previously implanted probes or by user defined target and entry points;

define a trajectory for the probe based on a mean target coordinates and mean entry point coordinates;

adjust the trajectory to intersect with a nearest voxel assigned a label of an anatomical region of interest;

check proximity of the trajectory to critical structures based on user defined constraints and/or user defined modification of the trajectory to satisfy the user defined constraints; and

superimpose the trajectory on the second data set to form a planning dataset.

16. The system of claim 11 , wherein:

the first imaging scan is an anatomical scan, the second imaging scan is a post-implantation imaging scan, the first dataset is an anatomical dataset, and the second dataset is a post-implantation imaging dataset; and

the instructions configure the one or more processors to:

obtain a third imaging scan used to guide electrode implantation during surgery;

convert the third imaging scan to a third dataset;

align a third dataset with the first dataset;

obtain a trajectory implant data file created during the electrode implantation;

generate a planned trajectory dataset, wherein:

the planned trajectory dataset is based on trajectory implant data file; and

the planned trajectory dataset includes dummy objects disposed at locations of electrode geometry;

align the planned trajectory dataset to the post-implantation imaging dataset; and

identify electrodes in an electrode dataset based on the dummy objects of the trajectory implant data file.

Continuity (3)
Continuation 17181345 · Feb 22, 2021
Provisional Application 62978868 · Feb 20, 2020
Related Publication 20240221192A1 · Jul 4, 2024
References Cited (33)
US 7813581B1 · Fitzpatrick · 2010 [cited by examiner]
US 10149618B1 · Tandon · 2018 [cited by examiner]
US 20080292194A1 · Schmidt · 2008 [cited by examiner]
US 20100014732A1 · Vija · 2010 [cited by examiner]
US 20110295515A1 · Grady et al. · 2011 [cited by applicant]
US 20120093381A1 · Fan et al. · 2012 [cited by applicant]
US 20150141804A1 · Rooney et al. · 2015 [cited by applicant]
US 20160343127A1 · Miller · 2016 [cited by examiner]
US 20170147908A1 · Chen et al. · 2017 [cited by applicant]
US 20170193161A1 · Sapiro et al. · 2017 [cited by applicant]
US 20200116808A1 · Taher · 2020 [cited by examiner]
JP 201361196 · 2013 [cited by applicant]
JP 2014518516 · 2014 [cited by applicant]
JP 2016517288 · 2016 [cited by applicant]
JP 2019500179 · 2019 [cited by applicant]
WO WO2012135190 · 2012 [cited by applicant]
WO WO2014139024 · 2014 [cited by applicant]
WO WO2015010745 · 2015 [cited by applicant]
WO WO2017072706 · 2017 [cited by applicant]
Puonti, Oula (“Fast and sequence adaptive whole-brain segmentation using parametric Bayesian modeling” NeuroImage (Year: 2016). [cited by examiner]
Awate, Suyash P., et al. “Adaptive Markov modeling for mutual-information-based, unsupervised MRI brain-tissue classification.” [cited by applicant]
Conner, Christopher R., et al. “Category specific spatial dissociations of parallel processes underlying visual naming.” [cited by applicant]
Frangi, Alejandro F., et al. “Multiscale vessel enhancement filtering.” [cited by applicant]
Freesurfer. Retrieved from: https://surfer.nmr.mgh.harvard.edu/fswiki/mri_tessellate; Feb. 21, 2021. [cited by applicant]
International Search Report mailed May 4, 2021, issued in counterpart International Application No. PCT/ US2021/019014. [cited by applicant]
Kadipasaoglu, C. M., et al. “Surface-based mixed effects multilevel analysis of grouped human electrocorticography.” [cited by applicant]
Menze, Bjoern H., et al. “The multimodal brain tumor image segmentation benchmark (BRATS).” [cited by applicant]
Partial Supplementary European Search Report issued in European Application No. 21757383.1, dated Feb. 13, 2024. [cited by applicant]
Puonti, Oula, Juan Eugenio Iglesias, and Koen Van Leemput. “Fast and sequence-adaptive whole-brain segmentation using parametric Bayesian modeling.” [cited by applicant]
Puonti, Oula, Juan Eugenio Iglesias, and Koen Van Leemput. “Fast, sequence adaptive parcellation of brain MR using parametric models.” [cited by applicant]
Saad, Ziad S., and Richard C. Reynolds. “Suma.” [cited by applicant]
Vascular Modeling ToolKit software; Retrieved from VMTK: http://www.vmtk.org/vmtkscripts/vmtkmarchingcubes.html on Feb. 21, 2021. [cited by applicant]
Office Communication dated Oct. 29, 2024, issued in counterpart Japanese Application No. 2022-549631. [cited by applicant]