Dynamic magnetic resonance angiography method and system

1. A dynamic magnetic resonance angiography, MRA, method, comprising:

acquiring, by an MR scanning device, a multi-contrast magnetic resonance, MR, sequence of a portion of a body, through a single multi-contrast MR acquisition, without administering a contrast agent to the body;

wherein the multi-contract MR sequence comprises quantification information of the portion;

wherein the quantification information of the portion comprises a longitudinal relaxation rate R 1 , and a transverse relaxation rate R 2 ;

identifying, by a processing circuit, blood vessels of the portion by identifying blood of the portion by comparing predetermined characteristic of blood with the quantification information of the portion;

wherein the predetermined characteristic of blood comprises: a value of a longitudinal relaxation rate R 1 of blood, and a value of a transverse relaxation rate R 2 of blood;

generating, by the processing circuit, a first MRA image frame and a second MRA image frame, based on the multi-contrast MR sequence, respectively visualising a first part and a second part of the identified blood vessels;

generating, by the processing circuit, a dynamic MRA image for visualising a dynamic blood flow through a part of the portion, based on the first and second MRA image frame;

generating a plurality of multi-contrast images based on the multi-contrast MR sequence; and

identifying the blood vessels of the portion by identifying the blood of the portion based on the predetermined characteristic of blood and the plurality of multi-contrast images,

wherein the step of identifying blood vessels comprises:

calculating quantification values of each voxel of a plurality of voxels representing the portion, based on the plurality of multi-contrast images;

wherein each voxel of the plurality of voxels corresponds to a volume of the portion, and

wherein the quantification values comprise: a value of a longitudinal relaxation rate R 1 and a value of a transverse relaxation rate R 2 ;

selecting a first group of voxels of the plurality of voxels representing the identified blood vessels, by comparing the quantification values of each voxel with the predetermined characteristic of blood.

2. The method of claim 1 , wherein the step of generating a first MRA image frame comprises:

determining a first group of values corresponding to the first part of the identified blood vessel to be visualised in the first MRA image frame,

wherein the first group of values comprise a value of longitudinal relaxation rate R 1 and a value of transverse relaxation rate R 2 ;

selecting a second group of voxels of the plurality of voxels representing the first part of the identified blood vessel by comparing the quantification values of each voxel with the first groups of values;

generating the first MRA image frame visualising the first part of the identified blood vessels, based on the second group of voxels.

3. The method of claim 1 , wherein the quantification information of the portion comprises a Proton Density, PD;

wherein the predetermined characteristic of blood comprises: a value of a PD of blood.

4. The method of claim 3 , wherein the quantification values of each voxel comprise: a value of a PD.

5. The method of claim 4 , wherein the first group of values comprise: a value of a PD.

6. The method of claim 1 , wherein the step of identifying blood vessels comprises:

selecting a part of the multi-contrast MR sequence representing the identified blood vessels, based on the predetermined characteristic of blood;

wherein the predetermined characteristic of blood comprises: a representation of a signal intensity of blood.

7. The method of claim 6 , wherein the step of generating a first MRA image frame comprises:

selecting a first group of data of the selected part of the multi-contrast MR sequence corresponding to the first part of the identified blood vessel;

generating the first MRA image frame visualising the first part of the identified blood vessels, based on the multi-contrast MR sequence and the first group of data.

8. The method of claim 1 , wherein the step of generating a first and a second MRA image frame comprises:

visualising the first MRA image frame and the second MRA image frame, by a Maximum Intensity Projection, MIP.

9. The method of claim 8 , wherein the step of generating a first and a second MRA image frame comprises:

generating a synthetic MR image representing the portion based on the multi-contrast MR sequence, as a static background image of the dynamic MRA image;

generating the dynamic MRA image based on the first, the second MRA image frame, and the synthetic MR image.

10. The method of claim 1 , wherein the first part and the second part of the identified blood vessels are:

at least partially overlapping, or

nonoverlapping.

11. The method of claim 1 , wherein the first part of the identified blood vessels comprises early inflow arteries of the identified blood vessels, and

the second part of the identified blood vessels comprises late inflow arteries of the identified blood vessels.

12. The method of claim 1 , wherein the first part of the identified blood vessels comprises arteries of the identified blood vessels.

13. The method of claim 12 , wherein the second part of the identified blood vessels comprises all the identified blood vessels, or veins of the identified blood vessels.

14. The method of claim 1 , wherein the dynamic MRA image consecutively displays the first and second MRA image frame.

15. The method of claim 14 , wherein at least a first portion of the first MRA image frame and a second portion of the second MRA image frame are displayed in different colours or different brightnesses.

16. The method of claim 1 , further comprising:

generating at least a third MRA image frame visualising at least a third part of the identified blood vessels;

generating the dynamic MRA image based on the first, second and at least the third MRA image frame,

wherein the first part, second part and at least the third part of the identified blood vessels are:

at least partially overlapping, or

nonoverlapping.

17. The method of claim 1 , wherein the value of the longitudinal relaxation rate R 1 of blood is in a range of 0.4-0.6 s −1 ; and

the value of transverse relaxation rate R 2 of blood is in a range of 20-30 s −1 .

18. The method of claim 1 , wherein the portion of the body comprises a head and/or a neck.

19. The method of claim 1 , wherein each voxel of the plurality of voxels has an intensity value,

wherein the method further comprises:

calculating a perfusion value by summing the intensity value of each voxel of the first group of voxels.

20. The method of claim 1 , wherein each voxel of the plurality of voxels has an intensity value,

wherein the method further comprises:

calculating a first sum by summing intensity values of voxels of the plurality of voxels representing the first part of the identified blood vessels;

calculating a second sum by summing intensity values of voxels of the plurality of voxels representing the second part of the identified blood vessels;

calculating a rate of change of the first sum and second sum as a measure of a mean transit time of blood through the portion of the body.

21. The method of claim 1 , wherein the step of generating a first and a second MRA image frame comprises:

visualising the first MRA image frame and the second MRA image frame, by a Minimum Intensity Projection, MinIP.

22. A non-transitory computer readable recording medium having computer readable program code recorded thereon which when executed on a device having processing capability is configured to perform the method of claim 1 .

23. A dynamic magnetic resonance angiography, MRA, system, comprising:

an MR scanning device configured to:

acquire a multi-contrast magnetic resonance, MR, sequence of a portion of a body, through a single multi-contrast MR acquisition;

wherein the multi-contrast MR sequence comprises quantification information of the portion;

wherein the quantification information of the portion comprises a longitudinal relaxation rate R 1 , and a transverse relaxation rate R 2 ;

a processing circuit configured to:

identify blood vessels of the portion by identifying blood of the portion by comparing predetermined characteristic of blood with the quantification information of the portion;

wherein the predetermined characteristic of blood comprises: a value of a longitudinal relaxation rate R 1 of blood, and a value of a transverse relaxation rate R 2 of blood;

generate a first MRA image frame and a second MRA image frame, based on the multi-contrast MR sequence, respectively visualising a first part and a second part of the identified blood vessels;

generate a dynamic MRA image for visualising a dynamic blood flow through a part of the portion, based on the first and second MRA image frame;

generate a plurality of multi-contrast images based on the multi-contrast MR sequence; and

identify the blood vessels of the portion by identifying the blood of the portion based on the predetermined characteristic of blood and the plurality of multi-contrast images,

wherein the processing circuit is further configured to:

calculate quantification values of each voxel of a plurality of voxels representing the portion, based on the plurality of multi-contrast images;

wherein each voxel of the plurality of voxels corresponds to a volume of the portion, and wherein the quantification values comprise: a value of a longitudinal relaxation rate R 1 and a value of a transverse relaxation rate R 2 ; and

select a first group of voxels of the plurality of voxels representing the identified blood vessels, by comparing the quantification values of each voxel with the predetermined characteristic of blood.

24. The system as claimed in claim 23 , further comprising:

a user interface configured to display the dynamic MRA image.