Image enhancement method and system for fiducial-less tracking of treatment targets

• US 7,522,779 B2
• Filed: 06/30/2004
• Issued: 04/21/2009
• Est. Priority Date: 06/30/2004
• Status: Expired due to Fees

• Current Assignee: Accuray Incorporated
• Original Assignee: Accuray Incorporated
• Inventors: Kuduvalli, Gopinath, Fu, Dongshan
• Primary Examiner(s): Wu; Jingge
• Assistant Examiner(s): TSAI, TSUNG YIN
• Application Number:
  US10/881,208
• Publication Number:
  US 20060002615A1
• Time in Patent Office: 1,756 Days
• Field of Search: 382/254, 382/172
• US Class Current: 382/254
• CPC Class Codes: A61B 6/505 for diagnosis of bone
• CPC Class Codes: A61B 6/5235 combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
• CPC Class Codes: A61N 2005/1061 using an x-ray imaging system having a separate imaging source
• CPC Class Codes: A61N 2005/1062 using virtual X-ray images, e.g. digitally reconstructed radiographs [DRR]
• CPC Class Codes: A61N 5/1049 for verifying the position of the patient with respect to the radiation beam
• CPC Class Codes: G06T 2207/10116 X-ray image
• CPC Class Codes: G06T 2207/30004 Biomedical image processing
• CPC Class Codes: G06T 5/30 Erosion or dilatation, e.g. thinning
• CPC Class Codes: G06T 5/94 based on local image properties, e.g. for local contrast enhancement
• CPC Class Codes: G06T 7/33 using feature-based methods
• CPC Class Codes: G06T 7/35 using statistical methods
• CPC Class Codes: G06T 7/38 Registration of image sequences

Abstract

A method and system are presented for enhancing one or more images of an object, so as to increase the visibility within the images of one or more structures within the object. The object may be an anatomical region of a patient, and may include one or more reference structures, for example skeletal structures or vertebral structures, and one or more treatment targets, for example tumors or lesions. An operator, for example a top-hat filter operator, selects at least a first neighborhood and a second neighborhood within the images. The operator selects within each neighborhood one or more pixels having an optimal pixel value, and eliminates the remaining pixels in these neighborhoods. When the operator is applied to the selected neighborhoods, only the pixels having the greatest pixel values remain in the selected neighborhoods, and the remaining pixels are eliminated in the selected neighborhoods. As a result, desired features can be located and enhanced in the images.

Claims

1. A method of enhancing an image of an object so as to increase the visibility in the image of at least one structure within the object, wherein the image is characterized by a plurality of pixels, each pixel having an associated pixel value that represents the image intensity of a corresponding unit volume of the object, the method comprising:

a) selecting at least a first neighborhood and a second neighborhood within the image,

b) constructing a non-linear top-hat filter operator configured to select, within the first and second neighborhoods, one or more pixels having an optimal pixel value, and to eliminate the remaining pixels in the neighborhoods, wherein a mathematical expression for the top-hat filter operator comprises

f_e=f+w×

[f−

γ

_B(f)]−

b×

[φ

_B(f)−

f],

wherein f represents the image, f_e represents an enhanced image resulting from the application of the top-hat filter to selected pixels within the image, w and b represent weighting coefficients, γ

_B(f) represents a first structural element for the opening of the image f, and φ

_B(f) represents a second structural element for the closing of the image f; and

c) applying the non-linear top-hat filter operator to the selected neighborhoods so that only the pixels having the optimal pixel values remain in the selected neighborhoods, and the remaining pixels in the selected neighborhoods are eliminated.

2. A method in accordance with claim 1, wherein the act of constructing the operator comprises:

determining empirical values for the first structural element γ

_B(f) and the second structural element φ

_B(f); and

adaptively determining the values of the weighting coefficients w and b using the amplitudes of the first and second structural elements.

3. A method in accordance with claim 2, wherein the object comprises an anatomical region of a patient, and the structure comprises at least one skeletal structure and at least one treatment target.

4. A method in accordance with claim 3, wherein the anatomical region comprises at least one of:

a cervical region of a spine;

a lumbar region of the spine; and

a thoracic region of the spine; and

wherein the reference structure comprises a skeletal structure.

5. A method in accordance with claim 3, wherein the anatomical region comprises a cervical region of a spine, and wherein the value of the weighting coefficient w is about 1, and the value of the weighting coefficient b is about 1.

6. A method in accordance with claim 3, wherein the anatomical region comprises a lumbar region of a spine, and wherein the value of the weighting coefficient w is greater than 2, and the value of the weighting coefficient b greater than 2.

7. An apparatus for enhancing an image of an object so as to increase the visibility in the image of at least one structure within the object, the image being characterized by a plurality of pixels, each pixel having an associated pixel value that represents the image intensity of a corresponding unit volume of the object, the apparatus comprising:

a neighborhood selector configured to select at least a first neighborhood and a second neighborhood within the image;

an operator generator for constructing an operator configured to select, within the first and second neighborhoods, one or more pixels having an optimal pixel value, and to

eliminate the remaining pixels in the selected neighborhoods, wherein a mathematical operation for the operator comprises;

f_e=f+w×

[f−

γ

_B(f)]−

b×

[φ

_B(f)−

f],

wherein f represents the image, f_e represents an enhanced image resulting from the application of the operator to one or more pixels within the image, w and b represent weighting coefficients, γ

_B(f) represents the first structural element for the opening of the image f, and φ

_B(f) represents the second structural element for the closing of the image f; and

c) a controller for applying the operator to the selected neighborhoods so that only the pixels having the optimal pixel values remain in the selected neighborhoods, and the remaining pixels in the selected neighborhoods are eliminated.

8. An apparatus in accordance with claim 7,

wherein the first structural element γ

_B(f) and the second structural element φ

_B(f) are determined empirically, and

wherein the weighting coefficients w and b are determined adaptively using the amplitudes of the empirically determined first and second structural elements.