Automated method for digital image quantitation
First Claim
1. A method for quantitatively analyzing a diagnostic image of a left ventricle including the steps of:
- a. generating a diagnostic image frame in a digital format containing rows and columns of pixels of a left ventricle having a posterior epicardial border, and an anterior epicardium;
b. determining a first tentative center-point for the left ventricle and the position of the posterior epicardial border of the left ventricle depicted in the diagnostic image by frame repeatedly filtering the image frame with a set of circular arc filters until a maximum value for the set of filters is obtained;
c. determining a second tentative center-point for the left ventricle and the position of the anterior epicardium of the left ventricle depicted in the diagnostic image frame by repeatedly filtering the image frame with a set of matched elliptical arc filters until a maximum value for the set of matched filters is obtained;
d. determining a final center-point for the left ventricle depicted in the diagnostic image frame as midway between the anterior epicardial border and the posterior epicardial border along a line through the first and second tentative center-points; and
e. indicating the position of the final center-point in connection with the diagnostic image frame of the left ventricle.
0 Assignments
0 Petitions
Accused Products
Abstract
A method for quantitatively analyzing digital images of approximately elliptical body organs, and in particular, echocardiographic images. In a first preferred embodiment, the invention automatically determines the center-point of an elliptical organ and a search region for border discrimination by finding the maximum value produced by a series of circular arc filters applied to the image data. In a second preferred embodiment, a pair of elliptical arc filters are used to make a first approximation of the positions of the right and left epicardium so that matched elliptical filters can be used to estimate the anterior epicardial position. The estimated position of the anterior epicardial, together with the posterior epicardial position, determines the center-point of the left ventricle. The positions of the organ borders are then approximated by generating amplitude distributions along each of a number of circular segments based on the center-point. The approximated borders are then compared on a best-fit basis to a set of paired elliptical arcs used to model the borders. The image may then be displayed with the modeled border superimposed on the image. A method of automated video densitometry is also disclosed whereby a sequence of end diastolic images are obtained before and after injection of a contrast agent.
-
Citations
38 Claims
-
1. A method for quantitatively analyzing a diagnostic image of a left ventricle including the steps of:
-
a. generating a diagnostic image frame in a digital format containing rows and columns of pixels of a left ventricle having a posterior epicardial border, and an anterior epicardium; b. determining a first tentative center-point for the left ventricle and the position of the posterior epicardial border of the left ventricle depicted in the diagnostic image by frame repeatedly filtering the image frame with a set of circular arc filters until a maximum value for the set of filters is obtained; c. determining a second tentative center-point for the left ventricle and the position of the anterior epicardium of the left ventricle depicted in the diagnostic image frame by repeatedly filtering the image frame with a set of matched elliptical arc filters until a maximum value for the set of matched filters is obtained; d. determining a final center-point for the left ventricle depicted in the diagnostic image frame as midway between the anterior epicardial border and the posterior epicardial border along a line through the first and second tentative center-points; and e. indicating the position of the final center-point in connection with the diagnostic image frame of the left ventricle. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 25)
-
-
18. A method for quantitatively analyzing a diagnostic digital image of the left ventricle, including the steps of:
-
a. determining a first tentative center-point for and a region of search for the posterior epicardial border of the left ventricle depicted in the diagnostic image by repeatedly filtering the diagnostic image with a set of circular arc filters until a maximum value for the set of filters is obtained; b. determining a second tentative center-point for and a region of search for the anterior endocardial and epicardial borders of the left ventricle depicted in the diagnostic image by repeatedly filtering the diagnostic image with a set of coupled elliptical arc filters until a maximum value for the set of coupled filters is obtained; c. determining a final center-point for the left ventricle depicted in the diagnostic image as midway between the anterior epicardial border and the posterior epicardial border along a line through the first and second tentative center-points; and d. indicating the position of the final center-point in connection with the diagnostic image.
-
-
19. A method for quantitatively analyzing a diagnostic digital image of the heart, said images containing rows and columns of pixels, said method including the steps of:
-
a. identifying a maximum pixel value in a portion of said image where an interface of myocardium and lung is expected; b. selecting two threshold values t1 and t2 which are fractions of said maximum pixel value; c. creating a first image having pixel values by steps comprising; i. assigning a value of zero to locations in the image where the pixel values are less than t1; ii. assigning a value of 1 to locations in the image where the pixel values are between t1 and t2; and iii. assigning a value of 2 to locations in the image where the pixel values are greater than t2; d. determining an expected vertical diameter of said myocardium; e. summing the gray level values of pixel locations along a 180 circular arc having a radius r which varies from 1/4 and 3/4 of said expected diameter; f. calculating a totalsum=sum(r)+2sum(r+step)-sum(r-step)-sum(r-2step)-(r-3step)-sum2(r-step), where sum(r) is the sum of the pixel values along the circular arc, sum2(r) is the sum of the pixel values along a 90 subarc symmetric with respect to the vertical centerline of the left ventricle; g. calculating t=totalsum*min(r, expectedradius); and h. comparing the t values for each row and selecting the row with the highest t value as the anterior epicardial border. - View Dependent Claims (20, 21, 22, 23, 24)
-
-
26. A method of calculating a figure of merit for a left ventricle comprising the steps of:
-
a. defining a small region within a chamber cavity which estimates the position of the epicardium, said region defined having a lower limit equal to the row of the center-point of the left ventricle, an upper limit equal to 1/2 of the anterior radius, and having a width equal to the columns at ±
1/4 the radius to the left and to the right of the column of the center-point;b. calculating the average intensity value of the region; and c. comparing said average to the intensity values of the elliptical arcs defining the epicardium wherein if the epicardium intensity data is less than or equal to the chamber gray level within its octant, then the epicardial data is considered to be poor, and wherein if the gray level of the epicardium within an octant is greater than the chamber mean but less than the mean gray level for the octant plus 1/2 of its standard deviation, then epicardial data is questionable, and wherein if the epicardial gray level for the octant is greater than the mean of its octant plus of its standard deviation, then the estimated position is likely good. - View Dependent Claims (27, 28, 29)
-
-
30. A method of determining a wall motion using a set of time varying echocardiographic image frame comprising the steps of:
-
a. determining the center-point of the left ventricle; b. forming eight one-dimensional amplitude distributions for each of eight 45°
sectors emanating from the center-point; andc. obtaining the cross-correlation for each said amplitude distributions between the frame at end-diastole and the frame at end-systole to detect the wall motion in each of the different regions wherein the amount of radial shift producing the maximum cross-correlation gives an estimate of the amount of wall motion in each sector between the two time frames studied.
-
-
31. A method of detecting dropout across the interventricular septum between the left and right ventricles comprising the steps of:
-
a. comparing the minimum average grey level of a small image region in the cavity of a left ventricle around the center-point with the maximum average grey level of a set of similarly sized regions chosen along a line from the center-point of the ventricle in the direction of the septum; and b. wherein if the maximum ratio of the peak level to the background level for the series of comparisons is less than two, signal dropout is considered to have been detected.
-
-
32. Method of automated video densitometry comprising the steps of:
-
a. obtaining control frames from a baseline period; b. identifying the sides and last row of a sector scan; c. detecting the approximate position of the endocardial and epicardial borders by using the endocardial borders derived during the control period as the region of interest for the left ventricular chamber; d. injecting a contrast agent; e. obtaining post-contrast frames; f. estimating the position of the posterior epicardium on each frame by computing a posterior circular arc filter; g. translating the regions of interest based upon the difference in position of the posterior circular arc filter and the position of the posterior circular arc filter in the first control frame; h. calculating the mean pixel intensity of the left ventricular chamber for the pixels enclosed by the approximated endocardial border, and similarly calculating the mean pixel intensity of the right ventricular outflow tract; i. comparing the mean pixel intensity for the sequence of end diastolic frames with the background pixel intensity means from the time-mean brightness curve for each injection; j. segmenting the regions of the myocardium into regions approximating distributions of coronary arteries; and k. estimating whether perfusion is present in a region by comparing a pre-contrast mean pixel brightness to that calculated after injection of the contrast agents.
-
-
33. A method for quantitatively analyzing a diagnostic image of a left ventricle including the steps of:
-
a. generating a diagnostic image frame of a left ventricle having a posterior epicardial border and an anterior epicardial border; b. determining a first tentative center-point (x0, y0) for said left ventricle and the position of said posterior epicardial border by repeatedly filtering the image frame with a set of circular arc filters until a maximum value for the set of filters is obtained; c. determining a second tentative center-point (x1, y1) for said left ventricle and the position of said anterior epicardium by repeatedly filtering the image frame with a set of matched elliptical arc filters until a maximum value for the set of matched filters is obtained; and d. determining a final center-point for said left ventricle by using said first and second tentative center-points and weighting functions. - View Dependent Claims (34, 35, 36)
-
-
37. A method for determining the optical density of a digital image of a left ventricle said method comprising the steps of:
-
a. obtaining a control frame from a baseline period prior to injecting any contrast agent; b. eliminate any patient data or irrelevant information from said image; c. approximating the position of the endocardial and epicardial borders on said control frame; d. establishing a region of interest for said right ventricular outflow tract; e. injecting contrast agent and obtaining an ED frame; f. calculating the mean pixel intensity of the left ventricular chamber for the pixels enclosed by the approximate endocardial contour, and likewise calculating the mean pixel intensity of the right ventricular chamber to obtain right and left ventricular time-brightness curves; g. determining the maximum right ventricular mean pixel brightness and the area under the right ventricular time-brightness curve above background; h. determining the minimum left ventricular mean pixel brightness and the area under the left ventricular time-brightness curve above background; i. obtaining a time-mean pixel brightness curve by comparing the mean pixel brightness for said ED frame following injection of the contrast agent to the background frame; j. approximating the epicardial and endocardial contours of the peristomal short axis view; and k. calculating the mean pixel brightness of the myocardium of the left ventricle as the mean pixel brightness of the pixel between the epicardial and endocardial borders. - View Dependent Claims (38)
-
Specification