Echogenicity quantification method and calibration method for ultrasonic device using echogenicity index

• US 10,249,037 B2
• Filed: 12/01/2014
• Issued: 04/02/2019
• Est. Priority Date: 01/25/2010
• Status: Active Grant

• Current Assignee: Amcad Biomed Corporation
• Original Assignee: Amcad Biomed Corporation
• Inventors: Chang, King-Jen, Chen, Chiung-Nein, Wu, Ming-Hsun, Chen, Argon, Wang, Yu-Hsin, Huang, Kuo-Chen
• Primary Examiner(s): Sorey, Robert A
• Application Number:
  US14/557,013
• Publication Number:
  US 20150086094A1
• Time in Patent Office: 1,583 Days
• CPC Class Codes: G06T 2207/10132 Ultrasound image
• CPC Class Codes: G06T 2207/20104 Interactive definition of region of interest [ROI]
• CPC Class Codes: G06T 2207/30096 Tumor; Lesion
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G06T 7/41 based on statistical description of texture

Abstract

An echogenicity quantification method and a calibration method for ultrasonic device using echogenicity index are disclosed. The method includes: receiving an ultrasound image which comprises a plurality of grayscale pixels; choosing a region of interest (ROI); calculating the values of the grayscale pixels in ROI to obtain average value and standard deviation; excluding pixels in ROI with the grayscale value smaller than the sum of the minimum value in the grayscale image and the product of a first scaling factor and the standard deviation, and larger than the sum of the average value and the product of a second scaling factor and the standard deviation; averaging the values of the remaining grayscale pixels in ROI to obtain an average value of interest; choosing a reference region; averaging the values of the remaining grayscale pixels in the reference region to obtain an average value of reference; and calculating the difference between the average value of interest and the average value of reference to obtain an echogenicity index.

Claims

1. An echogenicity quantification method, comprising:

receiving an ultrasound image by an ultrasonic imaging device, and the ultrasound image comprises a plurality of pixel points, each pixel point has a grayscale value;

choosing a region of interest (ROI) In the ultrasound image;

averaging the grayscale values of the pixel points in the ROI to obtain an average value of the ROI and calculating the grayscale scale values of the pixel points in the ROI to obtain a standard deviation of the ROI;

excluding the pixel points with the grayscale values smaller than the sum of the minimum value in the grayscale image and the product of a first scaling factor of a first constant and the standard deviation of the ROI from the ROI;

wherein the first scaling factor is predetermined;

excluding the pixel points with the grayscale values larger than the sum of the average value of the ROI and the product of a second scaling factor of a second constant and the standard deviation of the ROI from the ROI;

wherein the second scaling factor is predetermined;

averaging the grayscale values of the remaining pixel points in the ROI to obtain an average value of interest;

choosing a reference region which is nonoverlapped with the ROI in the ultrasound image;

averaging the grayscale values of the remaining pixel points in the reference region to obtain an average value of reference;

calculating the difference between the average value of interest and the average value of reference to obtain an echogenicity index by a computing device;

setting the grayscale values of an image area outside the ROI to a background value which is related with visual sensitivity by a display device; and

setting the grayscale values of the pixel points within the ROI to a value, wherein the value is the product of the background value and the average value of interest divided by the average value of reference by the display device.

2. The quantification method of claim 1, wherein the first scaling factor is a value within a range of 0 to 0.5, and the second scaling factor is a value within a range of 2 to 5.

3. The quantification method of claim 2, further comprising:

excluding the pixel points with the grayscale values larger than the sum of the average value of the reference region and the product of a fourth scaling factor and the standard deviation of the reference region from the reference region.

4. The quantification method of claim 3 wherein the third scaling factor is a value within a range of 0 to 0.5, and the fourth scaling factor is a value within a range of 0 to 5.

5. The quantification method of claim 1, further comprising:

calculating the grayscale values of the pixel points in the reference region to obtain an average value of the reference region and a standard deviation of the reference region; and

excluding the pixel points with the grayscale values smaller than the sum of the minimum value the grayscale values in the grayscale image and a product of a third scaling factor and the standard deviation of the reference region from the reference region.

6. The quantification method of claim 1, wherein the reference region is situated in an image area which is upside of the centroid of the ROI and not in the ROI.

7. The quantification method of claim 1, further comprising:

calibrating the grayscale values in the ultrasound image by a transform function.

8. The echogenicity quantification system of claim 1, the ultrasonic imaging device further:

calibrating the grayscale values in the ultrasound image by a transform function.

9. A method for deriving the transform function of claim 7 or claim 8, comprising:

imaging a phantom with a plurality of echo intensity areas repeatedly by an ultrasonic imaging device to obtain a plurality of first phantom images having a plurality of pixel points to be corrected, wherein each pixel point to be corrected has a grayscale value;

calculating the grayscale values of the pixel points to be corrected corresponding to one of the echo intensity areas in each first phantom image to obtain an average to be corrected of the echo intensity area;

repeating the aforementioned step to obtain a plurality of averages to be corrected corresponding to different echo intensity areas; and

obtaining a plurality of baseline averages through a baseline ultrasonic imaging device, comparing the average to be corrected and the baseline average of the same echo intensity areas and deriving a transform function of the averages to be corrected and the baseline averages of different echo intensity areas by adding a minimum grayscale range and a maximum grayscale range to derive a trend line.

10. The method of claim 9, wherein the baseline averages are obtained by:

imaging the phantom with a plurality of echo intensity areas repeatedly by the baseline ultrasonic imaging device to obtain a plurality of second phantom images having a plurality of baseline pixel points, wherein each baseline pixel point has a grayscale value;

calculating the grayscale values of the baseline pixel points corresponding to one of the echo intensity areas in each second phantom image to obtain the baseline average of the echo intensity area; and

repeating the aforementioned step corresponding to different echo intensity areas.

11. An echogenicity quantification system, comprising:

an ultrasonic imaging device for receiving an ultrasound image, and the ultrasound image comprises a plurality of pixel points, each pixel point has a grayscale value;

an analyzing device of a region of interest (ROI) connects to the ultrasonic imaging device for choosing a region of interest (ROI) in the ultrasound image;

averaging the grayscale values of the pixel points in the ROI to obtain an average value of the ROI and calculating the grayscale scale values of the pixel points in the ROI to obtain a standard deviation of the ROI;

excluding the pixel points with the grayscale values smaller than the sum of the minimum value in the grayscale image and the product of a first scaling factor of a first constant and the standard deviation of the ROI from the ROI;

excluding the pixel points with the grayscale values larger than the sum of the average value of the ROI and the product of a second scaling factor of a second constant and the standard deviation of the ROI from the ROI; and

averaging the grayscale values of the remaining pixel points in the ROI to obtain an average value of interest;

wherein the first scaling factor and the second scaling factor are predetermined;

an analyzing device of a reference region connects to the ultrasonic imaging device for choosing a reference region which is nonoverlapped with the ROI in the ultrasound image; and

averaging the grayscale values of the remaining pixel points in the reference region to obtain an average value of reference;

a computing device connects to the analyzing device of the ROI and the analyzing device of the reference region for calculating the difference between the average value of interest and the average value of reference to obtain an echogenicity index; and

a display device connects to the computing device for setting the grayscale values of the an image area outside the ROI to a background value which is related with visual sensitivity;

setting the grayscale values of the pixel points within the ROI to a value, wherein the value is the product of the background value and the average value of interest divided by the average value of reference.

12. The echogenicity quantification system of claim 11, the analyzing device of the reference region further:

calculating the grayscale values of the pixel points in the reference region to obtain an average value of the reference region and a standard deviation of the reference region; and

excluding the pixel points with the grayscale values smaller than the sum of the minimum value of the grayscale values (Gmin) in the grayscale image and a product of a third scaling factor and the standard deviation of the reference region from the reference region.

13. The echogenicity quantification system of claim 12, the analyzing device of the reference region further:

excluding the pixel points with the grayscale values larger than the sum of the average value of the reference region and the product of a fourth scaling factor and the standard deviation of the reference region from the reference region.

14. The echogenicity quantification system of claim 13, wherein the third scaling factor is a value within a range of 0 to 0.5, and the fourth scaling factor is a value within a range of 0 to 5.

15. An echogenicity quantification system of claim 11, wherein the first scaling factor is a value within a range of 0 to 0.5, and the second scaling factor is a value within a range of 2 to 5.

16. The echogenicity quantification system of claim 11, wherein the reference region is situated in an image area which is upside of the centroid of the ROI and not in the ROI.