Magnetic resonance method and system forming an isotropic, high resolution, three-dimensional diagnostic image of a subject from two-dimensional image data scans

• US 6,984,981 B2
• Filed: 04/21/2005
• Issued: 01/10/2006
• Est. Priority Date: 03/31/2000
• Status: Expired due to Term

• Current Assignee: VirtualScopics, LLC (Koninklijke Philips N.V.)
• Original Assignee: VirtualScopics, LLC (Koninklijke Philips N.V.)
• Inventors: Parker, Kevin J., Totterman, Saara Marjatta Sofia, Tamez-Peña, José
• Primary Examiner(s): Shrivastav, Brij
• Assistant Examiner(s): Fetzner, Tiffany A.
• Application Number:
  US11/110,717
• Publication Number:
  US 20050184730A1
• Time in Patent Office: 264 Days
• Field of Search: 324300-322, 600/410
• US Class Current: 324/309
• CPC Class Codes: A61B 5/055 involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• CPC Class Codes: A61B 5/4528 Joints A61B5/4533, A61B5/4538 take precedence
• CPC Class Codes: G01R 33/4835 of multiple slices
• CPC Class Codes: G01R 33/5608 Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels image data processing or generation, in general G06T

Abstract

MRI scans typically have higher resolution within a slice than between slices. To improve the resolution, two MRI scans are taken in different, preferably orthogonal, directions. The scans are registered by maximizing a correlation between their gradients and then fused to form a high-resolution image. Multiple receiving coils can be used. When the images are multispectral, the number of spectral bands is reduced by transformation of the spectral bands in order of image contrast and using the transformed spectral bands with the highest contrast.

Claims

1. A method of forming an isotropic, high-resolution, three-dimensional diagnostic image of a subject from two-dimensional image data, the method comprising:

(a) scanning the subject with an imaging scanner device, in a first direction relative to the subject, in order to take image data of a first plurality of slices, the image data of the first plurality of slices being two-dimensional image data and having a low resolution in the first direction and a high resolution in directions orthogonal to the first direction;

(b) scanning the subject with said imaging scanner device in a second direction relative to the subject, which is different from the first direction in order to take image data of a second plurality of slices, the image data of the second plurality of slices being two-dimensional image data having at least one dimension substantially in common with the image data of the first plurality of slices and having a low resolution in the second direction and a high resolution in directions orthogonal to the second direction;

(c) registering the first plurality of slices with the second plurality of slices to define a matrix of isotropic, high-resolution voxels in image space, wherein the matrix has unknown high-resolution voxel values; and

(d) solving for the unknown high-resolution voxel values in the matrix defined in step (c) in accordance with the image data taken in steps (a) and (b) in order to form the isotropic high-resolution three-dimensional diagnostic image in the image space, wherein steps (a) and (b) are performed with magnetic resonance imaging.

4. The method of claim 1, wherein steps (a) and (b) are performed using MRI with a plurality of receiving coils.

5. The method of claim 4, wherein steps (a) and (b) are performed using at least three receiving coils.

6. The method of claim 5, wherein steps (a) and (b) are performed using at least four receiving coils.

7. The method of claim 1, wherein steps (a) and (b) are performed through MRI;

and the image data of the first plurality of slices and the image data of the second plurality of slices comprise image data having a plurality of spectral bands.

8. The method of claim 7, wherein step (d) comprises selecting a sub-plurality from the plurality of spectral bands to form the image.

9. The method of claim 8, wherein the sub-plurality from the plurality of spectral bands is selected by ranking the plurality of spectral bands in order of image contrast and selecting the spectral bands whose image contrast is highest.

10. The method of claim 9, wherein the plurality of spectral bands is ranked in order of image contrast by:

deriving a covariance matrix from the plurality of spectral bands;

deriving a set of orthogonal eigenvectors and a corresponding set of eigenvalues from the covariance matrix; and

ranking the orthogonal eigenvectors in order of their corresponding eigenvalues.

2. A method of forming an isotropic, high-resolution, three-dimensional diagnostic image of a subject from two-dimensional image data, the method comprising:

(a) scanning the subject with an imaging scanner device, in a first direction relative to the subject, in order to take image data of a first plurality of slices, the image data of the first plurality of slices being two-dimensional image data, having gradients and having a low resolution in the first direction and a high resolution in directions orthogonal to the first direction;

(b) scanning the subject with said imaging scanner device in a second direction relative to the subject, which is different from the first direction in order to take image data of a second plurality of slices, the image data of the second plurality of slices being two-dimensional image data having at least one dimension substantially in common with the image data of the first plurality of slices, having gradients and having a low resolution in the second direction and a high resolution in directions orthogonal to the second direction;

(c) registering the first plurality of slices with the second plurality of slices to define a matrix of isotropic, high-resolution voxels in image space, wherein the matrix has unknown high-resolution voxel values; and

(d) solving for the unknown high-resolution voxel values in the matrix defined in step (c) in accordance with the image data taken in steps (a) and (b) in order to form a correlated isotropic high-resolution three-dimensional diagnostic image in the image space, wherein the correlation is a correlation of the gradients of the image data of the first and second pluralities of slices.

3. The method of claim 2, wherein the correlation is a correlation of a sub-sample of the image data of the first plurality of slices with a sub-sample of the image data of the second plurality of slices, the sub-samples being taken in accordance with the gradients of the image data.

11. A system forming an isotropic, high-resolution, three-dimensional diagnostic image of a subject from two-dimensional image data, the system comprising:

scanning means for

(i) scanning the subject in a first direction relative to the subject in order to take image data of a first plurality of slices, the image data of the first plurality of slices being two-dimensional image data, having gradients and having a low resolution in the first direction and a high resolution in directions orthogonal to the first direction, and

(ii) scanning the subject in a second direction relative to the subject which is different from the first direction in order to take image data of a second plurality of slices, the image data of the second plurality of slices being two-dimensional image data having at least one dimension substantially in common with the image data of the first plurality of slices, having gradients and having a low resolution in the second direction and a high resolution in directions orthogonal to the second direction; and

computing means for

(i) registering the first plurality of slices with the second plurality of slices in order to define a matrix of isotropic, high-resolution voxels in image space, wherein the matrix has unknown high-resolution voxel values and

(ii) solving for the unknown high-resolution voxel values in the matrix defined by the computing means in accordance with the image data taken in the first and second directions by the scanning means and thereby form a correlated isotropic, high-resolution, three-dimensional, diagnostic image in the image space,

wherein the scanning means comprises an MRI scanner.

12. The system of claim 11, wherein the correlation is a correlation of the gradients on the image data of the first and second pluralities of slices.

13. The method of claim 12, wherein the correlation is a correlation of a sub-sample of the image data of the first plurality of slices with a sub-sample of the image data of the second plurality of slices, the sub-samples being taken in accordance with the gradients of the image data.

14. The system of claim 11, wherein:

the image data of the first plurality of slices have a low resolution in the first direction and a high resolution in directions orthogonal to the first direction;

the image data of the second plurality of slices have the low resolution in the second direction and the high resolution in directions orthogonal to the second direction; and

the computing means registers the first plurality of slices with the second plurality of slices such that the image has the high resolution in_—

all directions.

15. The system of claim 11, wherein the scanning means comprises an MRI scanner with a plurality of receiving coils.

16. The system of claim 15, wherein the plurality of receiving coils comprises at least three receiving coils.

17. The system of claim 16, wherein the plurality of receiving coils comprises at least four receiving coils.

18. The system of claim 11, wherein the image data of the first plurality of slices and the image data of the second plurality of slices comprise image data having a plurality of spectral bands.

19. The system of claim 18, wherein the computing means selects a sub-plurality from the plurality of spectral bands in order to form the image.

20. The system of claim 19, wherein the sub-plurality from the plurality of spectral bands is selected by ranking the plurality of spectral bands in order of image contrast and selecting the spectral bands whose image contrast is highest.

21. The system of claim 20, wherein the plurality of spectral bands is ranked in order of image contrast by:

deriving a covariance matrix from the plurality of spectral bands;

deriving a set of orthogonal eigenvectors and a corresponding set of eigenvalues from the covariance matrix; and

ranking the orthogonal eigenvectors in order of their corresponding eigenvalues.