Method and apparatus for automatic muscle segmentation in digital mammograms
First Claim
1. A method of identifying a pectoral boundary in a digital mammogram, said pectoral boundary being characterized by at least one parameter, comprising the steps of:
- determining a region of interest containing said pectoral boundary;
calculating gradient magnitudes inside said region of interest;
accumulating said gradient magnitudes into a parameter plane, represented by H(ρ
,θ
), according to a Hough transform; and
using information in said parameter plane for identifying said at least one parameter of said pectoral boundary;
wherein said step of using information in said parameter plane comprises the steps of;
selecting a parameter plane region of interest corresponding to the range of possible angles normal to said pectoral boundary;
creating a normalized parameter plane HN(ρ
,θ
) by normalizing the value H(ρ
,θ
) by a normalizing factor, said normalizing factor being inversely related to the number of gradient magnitude pixels capable of being accumulated at (ρ
,θ
);
identifying a highest ranking local peak of said normalized parameter plane HN(ρ
,θ
); and
using the location of said highest ranking local peak for determining said at least one parameter of said pectoral boundary; and
wherein said normalizing factor is determined for each pixel (ρ
,θ
) by the steps of;
backprojecting a line corresponding to (ρ
,θ
) into the image plane;
determining the length L of said line; and
setting said normalizing factor equal to 1/sqrt(L) if L is greater than a first number of pixels.
6 Assignments
0 Petitions
Accused Products
Abstract
A method and apparatus for the automatic segmentation of the pectoral muscle boundary in a digital mammogram is disclosed. Parameters of the pectoral boundary line are determined by calculating gradient magnitudes of the digital mammogram in a region of interest, accumulating the gradient magnitudes into a parameter plane according to a Hough transform, and using information in the parameter plane for identifying the pectoral boundary. Prior to being accumulated in the parameter plane, the gradient magnitudes are weighted by a monotonic, multi-bit weighting function designed to exhibit larger variations near more densely populated levels of gradient magnitudes, and designed to exhibit smaller variations near less densely populated levels of gradient magnitudes, for allowing greater robustness against variations in gradient magnitudes among different digital mammograms.
-
Citations
19 Claims
-
1. A method of identifying a pectoral boundary in a digital mammogram, said pectoral boundary being characterized by at least one parameter, comprising the steps of:
-
determining a region of interest containing said pectoral boundary; calculating gradient magnitudes inside said region of interest; accumulating said gradient magnitudes into a parameter plane, represented by H(ρ
,θ
), according to a Hough transform; andusing information in said parameter plane for identifying said at least one parameter of said pectoral boundary; wherein said step of using information in said parameter plane comprises the steps of; selecting a parameter plane region of interest corresponding to the range of possible angles normal to said pectoral boundary; creating a normalized parameter plane HN(ρ
,θ
) by normalizing the value H(ρ
,θ
) by a normalizing factor, said normalizing factor being inversely related to the number of gradient magnitude pixels capable of being accumulated at (ρ
,θ
);identifying a highest ranking local peak of said normalized parameter plane HN(ρ
,θ
); andusing the location of said highest ranking local peak for determining said at least one parameter of said pectoral boundary; and wherein said normalizing factor is determined for each pixel (ρ
,θ
) by the steps of;backprojecting a line corresponding to (ρ
,θ
) into the image plane;determining the length L of said line; and setting said normalizing factor equal to 1/sqrt(L) if L is greater than a first number of pixels. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12)
-
-
13. A method of identifying a boundary between a first portion and a second portion of a digitized radiograph, said second portion corresponding to a tissue area of generally higher density than said first portion, said boundary generally being a linear boundary, comprising the steps of:
-
determining a region of interest containing said boundary; calculating gradient magnitudes inside said region of interest; accumulating said gradient magnitudes into a parameter plane, represented by H(ρ
,θ
), according to a transform function;using information in said parameter plane for identifying said linear boundary; wherein said step of using information in said parameter plane comprises the steps of; selecting a parameter plane region of interest corresponding to the range of possible angles normal to said pectoral boundary; creating a normalized parameter plane HN(ρ
,θ
) by normalizing the value H(ρ
,θ
) by a normalizing factor, said normalizing factor being inversely related to the number of gradient magnitude pixels capable of being accumulated at (ρ
,θ
);identifying a highest ranking local peak of said normalized parameter plane HN(ρ
,θ
); andusing the location of said highest ranking local peak for determining said at least one parameter of said linear boundary; and wherein said normalizing factor is determined for each pixel (ρ
,θ
) by the steps of;backprojecting a line corresponding to (ρ
,θ
) into the image plane;determining the length L of said line; setting said normalizing factor equal to 1/sqrt(L) if L is greater than a first number of pixels. - View Dependent Claims (14, 15, 16)
-
-
17. A method of identifying a pectoral boundary in a digital mammogram D(x,y), said pectoral boundary being a linear boundary characterized by an offset parameter and an angle parameter, comprising the steps of:
-
determining a region of interest containing said pectoral boundary; calculating gradient magnitudes g(x,y) inside said region of interest; weighting said gradient magnitudes g(x,y) by a monotonic, multi-bit weighting function designed to exhibit larger variations near more densely populated levels of gradient magnitudes, and designed to exhibit smaller variations near less densely populated levels of gradient magnitudes, for producing weighted gradient magnitudes; accumulating said weighted gradient magnitudes into a parameter plane H(ρ
,θ
) according to a Hough transform ρ
=x cos (θ
)+y sin (θ
);creating a normalized parameter plane HN(ρ
,θ
) by normalizing the value H(ρ
,θ
) by a normalizing factor, said normalizing factor being inversely related to the number of gradient magnitude pixels capable of being accumulated at (ρ
,θ
);identifying a highest ranking local peak of said normalized parameter plane HN(ρ
,θ
), said highest ranking local peak being at a location (ρ
P,θ
P), wherein said angle parameter is set equal to θ
P and said offset parameter is set equal to ρ
P ;wherein said digital mammogram comprises a tissue portion including a pectoral muscle portion and a breast tissue portion, said pectoral boundary defining said pectoral muscle portion, said pectoral muscle portion encompassing a first corner of said digital mammogram, said step of determining a region of interest including the steps of; determining a center of mass of said tissue area; and selecting said region of interest above a line centered at said center of mass and having a predetermined slope.
-
-
18. A method of identifying a pectoral boundary in a digital mammogram, said pectoral boundary being characterized by at least one parameter, comprising the steps of:
-
determining a region of interest containing said pectoral boundary; calculating gradient magnitudes inside said region of interest; computing gradient angles inside said region of interest; accumulating only gradient magnitudes associated with a predetermined range of angles into a parameter plane according to a Hough transform; and using information in said parameter plane for identifying said at least one parameter of said pectoral boundary; wherein said digital mammogram comprises a tissue portion including a pectoral muscle portion and a breast tissue portion, said pectoral boundary defines a pectoral muscle portion of the digital mammogram, said pectoral muscle portion encompasses a first corner of said digital mammogram, wherein said predetermined range of gradient angles includes an angle corresponding to said first corner and a first range of angles that is approximately 90 degrees and said step of determining a region of interest includes the steps of; determining a center of mass of said tissue area; and selecting said region of interest above a line centered at said center of mass and having a predetermined slope.
-
-
19. A method of identifying a pectoral boundary in a digital mammogram, said pectoral boundary being characterized by at least one parameter, comprising the steps of:
-
determining a region of interest containing said pectoral boundary; calculating gradient magnitudes inside said region of interest; accumulating said gradient magnitudes into a parameter plane, represented by H(ρ
,θ
), according to a Hough transform; andusing information in said parameter plane for identifying said at least one parameter of said pectoral boundary; wherein said step of using information in said parameter plane comprises the steps of; selecting a parameter plane region of interest corresponding to the range of possible angles normal to said pectoral boundary; said parameter plane being represented by H(ρ
,θ
), wherein said step of using information in said parameter plane comprises the steps of;selecting a parameter plane region of interest corresponding to the range of possible angles normal to said pectoral boundary; creating a normalized parameter plane HN(ρ
, θ
) by normalizing the value H(ρ
,θ
) by a normalizing factor, said normalizing factor being inversely related to the number of gradient magnitude pixels capable of being accumulated at (ρ
,θ
);identifying a highest ranking local peak of said normalized parameter plane HN(ρ
,θ
);using the location of said highest ranking local peak for determining said at least one parameter of said pectoral boundary; and said step of identifying a highest ranking local peak further comprises the steps of; determining whether a first set of peaks of HN(ρ
, θ
) above a lower predetermined threshold TL can be found;if such a first set of peaks cannot be found, selecting a NULL value as the highest ranking local peak; if such a first set can be found, determining if a second set of peaks of HN(ρ
,θ
) above an upper predetermined threshold TH can be found;if such a second set of peaks can be found, determining the corresponding area A of a segmented pectoral muscle in image space corresponding to each of said second set of peaks, and selecting that peak having the greatest corresponding area A as the highest ranking local peak; if such a second set of peaks above TH cannot be found, selecting the largest of said first set of peaks as the highest ranking local peak.
-
Specification