Estimating and correcting for contrast agent extravasation in tissue perfusion imaging

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES),
• wherein the tail part of the apparent residue function is determined by one or more of:
• fitting a singly peaked, monotonously decreasing function, L(t), with zero asymptotes to a global maxima peak of the apparent residue function and defining the tail part of the apparent residue function as the part where the amplitude of the apparent residue function is consistently larger than the amplitude of L(t);
• estimating the approximate mean transit time (MTT) from the perfusion image data, and from this estimate, estimating a time, tfp, larger than MTT and using R′(t) for t>tfp as the tail part of the apparent residue function;
• finding the section of the apparent residue function following its peak amplitude with minimum relative variance (variance divided by mean value) over a predefined section length.

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES);
• and
• wherein Ka is estimated as being proportional to a mean value of the tail part of the apparent residue function.

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cps′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES);
• and
• wherein Ka is estimated as being proportional to an asymptote of the apparent residue function.

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES);
• and
• wherein the estimation further comprises estimating a second leakage rate constant kep quantifying an apparent rate of contrast agent reflux EES back to plasma based on the tail part of the apparent residue function.

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES);
• and
• an automated segmentation of imaged tissue into a tumor region and white and/or grey matter, and wherein a mean value of Ka are estimated for the tumor region and for the white and/or grey matter.

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, R′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES);
• and
• wherein Ka is allowed to take both positive and negative values, and wherein the method further comprises separating between T1- and T2*-dominant leakage by determining the sign of Ka, so that negative Ka values reflect T1-dominant leakage effects and positive Ka values reflect T2*-dominant leakage effects.

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES);
• and
• wherein the apparent residue function, Rps′(t), is estimated as:
• R
• ps
• ′
• ⁡
• (
• t
• )
• =
• R
• ps
• ⁡
• (
• t
• )
• +
• e
• =
• F
• ⁡
• [
• R
• ⁡
• (
• t
• 0
• )
• R
• ⁡
• (
• t
• 1
• )
• ⋮
• R
• ⁡
• (
• t
• N
• -
• 1
• )
• ]
• +
• K
• a
• ⁡
• [
• ⅇ
• -
• k
• ep
• ⁢
• t
• 0
• ⅇ
• -
• k
• ep
• ⁢
• t
• 1
• ⋮
• ⅇ
• -
• k
• ep
• ⁢
• t
• N
• -
• 1
• ]
• where e is an error term due to leakage, R(t) is a tissue specific residue function, F is proportional to tissue flow, N is a number of data-points, and kep is a rate constant from the EES to the intravascular space.

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES);
• and
• wherein the apparent residue function, Rps′(t), is estimated as:
• R
• ps
• ′
• ⁡
• (
• t
• )
• =
• R
• ps
• ⁡
• (
• t
• )
• +
• e
• =
• F
• ⁡
• [
• 1
• R
• ⁡
• (
• t
• 1
• )
• ⋮
• R
• ⁡
• (
• t
• C
• )
• 0
• 0
• ⋮
• 0
• ]
• +
• K
• a
• ⁡
• [
• 0
• 0
• ⋮
• 0
• 1
• exp
• ⁡
• (
• -
• K
• trans
• ⁢
• Δ
• ⁢
• ⁢
• t
• /
• v
• e
• )
• ⋮
• exp
• ⁡
• (
• -
• K
• trans
• ⁢
• N
• ⁢
• ⁢
• Δ
• ⁢
• ⁢
• t
• /
• v
• e
• )
• ]
• where e is an error term due to leakage, R(t) is a tissue specific residue function, F is proportional to tissue flow, Δt is a sampling interval, N is a number of data-points, tc is a time index corresponding to a capillary transit time of the CA, Tc, Ktrans is a transfer constant describing the CA flux from the IVS to the EES, and νe is a leakage space per unit volume.

A method for estimating contrast agent extravasation in contrast agent based perfusion imaging methods to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) by deconvolution of the signal with Cp′(t);
• estimating a first leakage rate constant Ka from a tail part of the apparent residue function, where Ka is quantifying an apparent rate of contrast agent extravasation from plasma to extravascular, extracellular space (EES);
• and
• estimating a factor to be used in correction of partial volume (PV) effects in contrast agent based perfusion MR or CT images to be performed by an electronic processor, the method comprising determining a PV correction factor, fpv, as a ratio between parts of the signal following the first passage of the contrast agent from one or more voxels comprising arterial blood and one or more voxels filled with venous blood, respectively.

A method of distinguishing between T1-dominant and T2*-dominant extravasation effects in perfusion imaging to be performed by an electronic processor, the method comprising:

• accessing perfusion image data comprising a signal related to contrast agent concentration as a function of time in a voxel during the first and consecutive passages of the contrast agent through the vascular system;
• representing the signal by a convolution, Rps′(t)Cp′(t), of an apparent residue function, Rps′(t), and an estimated contrast agent concentration in plasma, Cp′(t), and determining the apparent residue function Rps′(t) for each voxel by deconvolution of the signal with Cp′(t);
• determining a sign of the tail part of the apparent residue function;
• and
• generating an image of a parameter reflecting the determined sign and magnitude.