Method for determining midpoint coordinates of an image of a point-symmetrical structure

• US 8,005,274 B2
• Filed: 10/23/2006
• Issued: 08/23/2011
• Est. Priority Date: 10/24/2005
• Status: Expired due to Fees

• Current Assignee: Siemens Healthcare GMBH (Siemens AG)
• Original Assignee: Siemens AG
• Inventors: Gühring, Jens, Littmann, Arne
• Primary Examiner(s): Arana; Louis M
• Application Number:
  US11/584,991
• Publication Number:
  US 20070092131A1
• Time in Patent Office: 1,765 Days
• Field of Search: 382/123, 382/151, 382/154, 382/173, 324/312, 324/317
• US Class Current: 382/123
• 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
• CPC Class Codes: G01R 33/56518 due to eddy currents, e.g. caused by switching of the gradient magnetic field
• CPC Class Codes: G01R 33/56563 caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0 G01R33/56509, G01R33/56518, G01R33/56536 take precedence
• CPC Class Codes: G01R 33/56572 caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field G01R33/56509, G01R33/56518, G01R33/56536 take precedence
• CPC Class Codes: G01R 33/5659 caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field G01R33/56509, G01R33/56518, G01R33/56536 take precedence
• CPC Class Codes: G01R 33/583 Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency G01R33/246 takes precedence

Abstract

The invention relates to a method for determining midpoint coordinates of an image of a point-symmetrical structure stored in a volume data record, comprising: extracting a partial volume data record from the volume data record which contains the image of the point-symmetrical structure, selecting a set of points whose coordinates lie in the partial volume data record, determining a measure for each point of the set wherein the measure for each of the points of the set characterizes in each case the symmetry of the partial volume data record in relation to this point, selecting a point as a midpoint of the point-symmetrical structure that point of the set in relation to which the partial volume data record has the greatest symmetry on the basis of the previously determined measures, and storing the coordinates of the midpoint in a computer unit.

Claims

1. A method for determining coordinates of a midpoint of a medical image having a point-symmetrical structure stored in a volume data record, comprising:

extracting a portion of the volume data record containing the image of the point-symmetrical structure;

selecting a set of points having coordinates lying in the portion of the volume data record;

measuring a symmetry of the portion of the volume data record with respective to each of the set points;

identifying a point from the set points as the midpoint of the point-symmetrical structure based on the measuring; and

storing the coordinates of the midpoint in a computer unit.

2. The method as claimed in claim 1, wherein the volume data record is filtered before the extracting step.

3. The method as claimed in claim 1, wherein the set of points is a subset of voxels of the portion of the volume data record located around a volume central area of the portion of the volume data record.

4. The method as claimed in claim 3, wherein the midpoint is identified in a resolution which is less than a resolution predetermined by the voxels.

5. The method as claimed in claim 1, wherein the set of points lie on a straight line.

6. The method as claimed in claim 5,

wherein a plurality of two-dimensional layer pictures are extracted from the portion of the volume data record,

wherein a plurality of intensity-value-weighted centers of area are determined for each of the two-dimensional layer pictures respectively, and

wherein the straight line is positioned as a regression line at the intensity-value-weighted centers of area.

7. The method as claimed in claim 1, wherein the symmetry of the portion of the volume data record for a point of the set points is measured by:

mirroring the portion of the volume data record with respective to the point, and

evaluating a similarity between the portion of the volume data record and the mirrored portion of the volume data record.

8. The method as claimed in claim 7, wherein the portion of the volume data record comprises a voxel having an intensity value and the evaluating is performed by comparing the intensity value of the voxel with an intensity value of a mirrored voxel that is a mirroring of the voxel with resective to the point.

9. The method as claimed in claim 8, wherein the intensity value of the mirrored voxel is calculated by interpolating from intensity values of voxels surrounded the mirrored voxel.

10. The method as claimed in claim 8, wherein the similarity between the portion of the volume data record and the mirrored portion of the volume data record is calculated based on a difference from the intensity value of the voxel and the intensity value of the mirrored voxel.

11. The method as claimed in claim 1, wherein a plurality of images comprising a plurality of point-symmetrical structures respectively are stored in the volume data record and coordinates of a midpoint for each of the point-symmetrical structures are determined.

12. The method as claimed in claim 11,

wherein the volume data record is preprocessed prior to the extracting step,

wherein the point-symmetrical structures are differentiable from an interfering structure in the preprocessed volume data record, and

wherein a portion of the volume data record containing one of the images of one of the point-symmetrical structures respectively is automatically selected and extracted based on the preprocessed volume data record.

13. The method as claimed in claim 12, wherein the interfering structure with a geometry lying outside of a tolerance range of the point-symmetrical structures is removed from the volume data record during the preprocessing.

14. The method as claimed in claim 13, wherein the tolerance range is adapted to a geometry of the point-symmetrical structures.

15. The method as claimed in claim 13, wherein the interfering structure is removed from a two-dimensional layer picture of the volume data record.

16. The method as claimed in claim 13, wherein the interfering structure is removed from a three-dimensional volume data record.

17. The method as claimed in claim 11, wherein the point-symmetrical structures are a spherical structure.

18. The method as claimed in claim 11, wherein the point-symmetrical structures are arranged in a phantom.

19. The method as claimed in claim 18, wherein the phantom is a magnetic resonance tomography image.

20. A device for determining coordinates of a midpoint of a medical image having a point-symmetrical structure stored in a volume data record, comprising:

an image recording unit that records the volume data record;

a data processing unit that is configured to;

extract a portion of the volume data record containing the image of the point-symmetrical structure,

select a set of points whose coordinates lie in the portion of the volume data record,

measure a symmetry of the portion of the volume data record with respective to each of the set points, and

identify a point from the set points as the midpoint of the point-symmetrical structure based on the measuring; and

a storage unit that stores the coordinates of the midpoint.