# US 10,871,539 B2 · Determination of a joint probability distribution of radius and length of anisotropic pores from double pulsed field gradient MRI data

> Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt

- Number: US 10,871,539 B2
- Document ID: US10871539B2
- Status: Granted
- Kind: B2
- Application number: 15/529,412
- Filed: 2015-11-24
- Published: 2018-11-22
- Granted: 2020-12-22
- Earliest claimed: 2014-11-24
- Assignees: The United States of America, as represented by the Secretary, Department of Health and Human Services
- Applicants: The United States of America, as represented by the Secretary, Department of Health and Human Services
- Inventors: Peter J. Basser; Dan H. Benjamini
- Web page: https://patentgrep.com/patent/US10871539B2

## CPC

- G01R 33/567: gated by physiological signals {, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating}
- G01R 33/56341: {Diffusion imaging}
- G16H 10/40: for data related to laboratory analysis, e.g. patient specimen analysis
- G01R: MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES (indicating correct tuning of resonant circuits H03J3/12)
- G16H: HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA

## Abstract

Described herein are exemplary methods for estimating a nonparametric joint radius-length (R-L) distribution of an ensemble of porous elements represented generally by finite cylinders. Some described methods comprise estimating an eccentricity distribution of a group of anisotropic porous elements. For example, disclosed methods can be applied to estimate a nonparametric joint R-L distribution of injured axons in nervous tissue, muscle tissue, plant elements, or other porous materials. Employing a novel three dimensional (3-D) double pulsed-field gradient (d-PFG) magnetic resonance (MR) acquisition scheme, both the marginal radius and length distributions of a population of generally cylindrical porous elements can be obtained. The marginal radius and length distributions can then b e used to constrain and stabilize the estimate of the joint radius-length distribution. Using the marginal distributions as constraints allows the joint R-L distribution to be reconstructed from an underdetermined system (i.e., more variables than equations).

## Claims

### Claim 1

A method for estimating a joint radius-length distribution of an ensemble of pores, the method comprising:
receiving data obtained from three-dimensional (3-D) diffusion-based magnetic resonance (MR) acquisition;
determining a marginal radius distribution (MRD) and a marginal length distribution (MLD) for an ensemble of pores based on the data obtained from 3-D diffusion-based MR; and
determining a joint radius-length distribution for the ensemble of pores based on data obtained from 3-D diffusion-based MR and based on the MRD and the MLD;
wherein the received 3-D diffusion-based MR acquisition data comprises MR data associated with different orthogonal dimensions of the porous elements based on selective application of diffusion gradient vectors with components along different orthogonal dimensions;
wherein the diffusion gradient vectors comprise:
a first diffusion gradient applied in a first plane substantially perpendicular to a porous length dimension of the pores;
a second diffusion gradient applied in a first direction along an axis substantially parallel to the porous length dimension of the pores; and
a third diffusion gradient applied in a second plane perpendicular to the first plane and parallel to the first direction.

### Claim 2 (depends on 1)

The method of claim 1, wherein the joint radius-length distribution is determined based in part on a non-negative least-squared minimization to estimate relative volumetric fractions of pore sizes using the MRD and MLD as constraints.

### Claim 3 (depends on 1)

The method of claim 1, wherein the joint radius-length distribution, Ψ(R, L), is constrained so that:
∑
j
⁢
⁢
Ψ
⁡
(
R
,
L
j
)
=
Ψ
⁡
(
R
)
⁢
⁢
and
⁢
⁢
∑
i
⁢
⁢
Ψ
⁡
(
R
i
,
L
)
=
Ψ
⁡
(
L
)
.

### Claim 4 (depends on 1)

The method of claim 1, wherein the joint radius-length distribution is a nonparametric joint radius-length distribution.

### Claim 5 (depends on 1)

The method of claim 1, wherein 3-D diffusion-based MR data acquisition comprises:
applying dual pulsed field gradient (d-PFG) MR diffusion gradients in an x-y plane substantially perpendicular to a longitudinal length dimension of the pores to determine the MRD of the pores;
applying d-PFG MR diffusion gradients in a direction along a z axis substantially parallel to the longitudinal length dimension of the pores to determine the MLD of the pores; and
applying d-PFG MR diffusion gradients in an x-z plane of the pores to produce an MR signal that is influenced by both the length and the radius of the pores.

### Claim 6 (depends on 5)

The method of claim 5, wherein:
the pores are represented as infinite cylinders when applying d-PFG MR diffusion gradients in the x-y plane;
the pores are represented as parallel plates when applying d-PFG MR diffusion gradients in a direction along the z axis; and
the pores are represented by both length and radius when applying d-PFG MR diffusion gradients in the x-z plane.

### Claim 7 (depends on 1)

The method of claim 1, wherein estimating the nonparametric joint radius-length distribution comprises using a transfer matrix that is undetermined.

### Claim 8 (depends on 1)

The method of claim 1, wherein the pores comprise anisotropic finite-length porous elements.

### Claim 9 (depends on 1)

The method of claim 1, wherein estimating the joint radius-length distribution comprises estimating a joint radius-length eccentricity distribution of the anisotropic finite porous elements.

### Claim 10 (depends on 1)

The method of claim 1, wherein the pores comprise injured axons.

### Claim 11 (depends on 10)

The method of claim 10, wherein the injured axons exhibit local compartment eccentricity variation due to injury.

### Claim 12 (depends on 1)

The method of claim 1, wherein the pores comprise myocytes or other muscle tissue elements.

### Claim 13 (depends on 1)

The method of claim 1, wherein the pores comprise plant sieve tube elements.

### Claim 14 (depends on 1)

The method of claim 1, wherein the pores are modeled as capped cylinders having a finite length.

### Claim 15

A system comprising:
a magnetic resonance (MR) system operable to apply three-dimensional (3-D) diffusion-based MR to an ensemble of pores and acquire resulting MR signal data from the applied 3-D diffusion-based MR;
an MR data processing system configured to determine a marginal radius distribution (MRD) and a marginal length distribution (MLD) for the ensemble of pores based on the data obtained from 3-D diffusion-based MR, and configured to determine a joint radius-length distribution for the ensemble of pores based on data obtained from the applied 3-D diffusion-based MR and based on the MRD and the MLD;
wherein the MR system is operable to acquire MR signal data from different orthogonal dimensions of the pores by selectively applying diffusion gradient vectors along different orthogonal dimensions of the pores;
wherein the diffusion gradient vectors comprise:
a first diffusion gradient applied in a first plane substantially perpendicular to a porous length dimension of the pores;
a second diffusion gradient applied in a first direction along an axis substantially parallel to the porous length dimension of the pores; and
a third diffusion gradient applied in a second plane perpendicular to the first plane and parallel to the first direction.

### Claim 16 (depends on 1)

A storage device storing computer-executable instructions for causing a computer to perform the method of claim 1.

### Claim 17 (depends on 1)

The method of claim 1, wherein receiving data obtained from 3-D diffusion-based MR acquisition comprises receiving data obtained from 3-D dual pulsed field gradient (d-PFG) MR acquisition.

### Claim 18 (depends on 15)

The system of claim 15, wherein the 3-D diffusion-based MR comprises 3-D dual pulsed field gradient (d-PFG) MR.

### Claim 19 (depends on 1)

The method of claim 1, further comprising utilizing the joint radius-length distribution to quantify axonal injury following traumatic brain injury to reveal microstructural changes resulting from injury.