Adjunctive ultrasound processing and display for breast cancer screening

• US 9,861,342 B2
• Filed: 04/15/2013
• Issued: 01/09/2018
• Est. Priority Date: 11/24/2000
• Status: Expired due to Term

• Current Assignee: U-System Inc. (General Electric Company)
• Original Assignee: U-System Inc. (General Electric Company)
• Inventors: Wang, Shih-Ping, Chin, Donald, Rao, Fangyi, Karssemeijer, Nico
• Primary Examiner(s): Park, Patricia
• Application Number:
  US13/863,259
• Publication Number:
  US 20130289405A1
• Time in Patent Office: 1,730 Days
• Field of Search: None
• CPC Class Codes: A61B 6/463 characterised by displaying multiple images or images and diagnostic data on one display
• CPC Class Codes: A61B 6/502 for diagnosis of breast, i.e. mammography
• CPC Class Codes: A61B 6/5241 combining overlapping images of the same imaging modality, e.g. by stitching
• CPC Class Codes: A61B 6/5247 combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
• CPC Class Codes: A61B 6/563 involving image data transmission via a network
• CPC Class Codes: A61B 8/0825 for diagnosis of the breast, e.g. mammography
• CPC Class Codes: A61B 8/406 using means for diagnosing suspended breasts
• CPC Class Codes: A61B 8/4281 characterised by sound-transmitting media or devices for coupling the transducer to the tissue
• CPC Class Codes: A61B 8/4405 Device being mounted on a trolley
• CPC Class Codes: A61B 8/463 characterised by displaying multiple images or images and diagnostic data on one display composite display in devices for short-range imaging with acoustic waves G01S7/52074
• CPC Class Codes: A61B 8/464 involving a plurality of displays
• CPC Class Codes: A61B 8/465 adapted to display user selection data, e.g. icons or menus
• CPC Class Codes: A61B 8/466 adapted to display 3D data 3D imaging and stereoscopic displays in devices for short-range imaging with acoustic waves G01S15/8993 and G01S7/52068, respectively; 3D image rendering G06T15/00; manipulating 3D models for computer graphics G06T19/00
• CPC Class Codes: A61B 8/467 characterised by special input means
• CPC Class Codes: A61B 8/5215 involving processing of medical diagnostic data
• CPC Class Codes: A61B 8/5238 for combining image data of patient, e.g. merging several images from different acquisition modes into one image for image registration algorithms G06T7/30
• CPC Class Codes: A61B 8/565 involving data transmission via a network

Abstract

An adjunctive ultrasound mammography system and associated methods use an adjunctive ultrasound display configured for quick, intuitive, interactive viewing of volumetric ultrasound scans, displayed near a conventional x-ray mammogram display. Preferred navigations among a thick-slice image array, a selected enlarged thick-slice image, and planar ultrasound views are described, including a mode in which the planar ultrasound views are updated in real time as a cursor is moved across an active thick-slice image. In one example the thick-slice images are inverted prior to display, with non-breast areas of the image preferably segmented out and reset to dark. The inverted thick-slice images are of more familiar significance to radiologists as they are more like conventional x-ray mammograms and allow benign features to be more easily dismissed as compared to non-inverted thick-slice images. Preferred embodiments emphasize larger mass lesions and that compensate for mass lesions that straddle thick-slice region borders.

Claims

1. A method for computing a two-dimensional thick-slice ultrasound image from a volumetric ultrasound representation of a breast, said volumetric ultrasound representation comprising ultrasonic property values for voxels of the breast, said thick-slice ultrasound image comprising values for pixels and corresponding to a first slab-like subvolume of the breast lying between a first border plane and a second border plane thereof, comprising:

identifying for a first pixel location in the thick-slice ultrasound image a first voxel set corresponding to a voxel column in said volumetric ultrasound representation passing through the first pixel location and extending from the first border plane to the second border plane;

identifying a thickness of the first slab-like volume between the first border plane and the second border plane;

identifying a mass size of interest;

identifying a target size to slab thickness ratio based on the identified thickness and mass size of interest;

computing one or more statistical properties of said first voxel set;

identifying a pixel value based on the one or more statistical properties of the first voxel set and the target size to slab thickness ratio; and

setting an output value for the first pixel location to the identified pixel value.

2. The method of claim 1, further comprising computing a respective output value for each additional pixel location using one or more statistical properties of respective additional voxel sets, said one or more statistical properties incurring changes across different pixel locations in mass localities that are more significant for masses greater than the mass size of interest and that are less significant for masses smaller than said mass size of interest, mass lesions greater than said mass size of interest being emphasized and mass lesions smaller than said mass size of interest being de-emphasized in said thick-slice ultrasound image.

3. The method of claim 2, said identifying pixel value comprising:

computing a histogram of values of said first voxel set;

computing at least a portion of a cumulative distribution function from said histogram;

determining a first pixel level for which said cumulative distribution function is equal to the target size to slab thickness ratio times a coefficient value; and

setting said output value equal to said first pixel level, and wherein the method further comprises displaying the thick-slice ultrasound image with the output value for the first pixel location.

4. The method of claim 3, said coefficient value being a fixed value.

5. The method of claim 4, wherein the coefficient value is between 0.20 and 0.45.

6. The method of claim 3, said coefficient value being variable for different pixel locations in said thick-slice ultrasound image.

9. The method of claim 1, further comprising:

identifying for a second pixel location in the thick-slice ultrasound image a second voxel set corresponding to a voxel column in said volumetric ultrasound representation passing through the second pixel location and extending from the first border plane to the second border plane;

computing one or more statistical properties of said second voxel set; and

computing an output value for the second pixel location using said one or more statistical properties of said second voxel set, the target size to slab thickness ratio, and a coefficient value, and where the computed one or more statistical properties of the second voxel set are different than the computed one or more statistical properties of the first voxel set.

7. A method for computing a two-dimensional thick-slice ultrasound image from a volumetric ultrasound representation of a breast, said volumetric ultrasound representation comprising ultrasonic property values for voxels of the breast, said thick-slice ultrasound image comprising values for pixels and corresponding to a first slab-like subvolume of the breast lying between a first border plane and a second border plane thereof, the breast further having a second slab-like subvolume immediately adjacent to the first slab-like subvolume, the method comprising:

identifying for each pixel location in the thick-slice ultrasound image a first voxel set corresponding to a voxel column in said volumetric ultrasound representation passing through that pixel location and extending from the first border plane to the second border plane;

computing one or more statistical properties of said first voxel set;

computing an output value for that pixel location using said one or more statistical properties of said first voxel set;

identifying for each pixel location in the thick-slice ultrasound image a second voxel set corresponding to said voxel column extending from a first intermediate elevation in said first slab-like subvolume to a second intermediate location in said second slab-like subvolume, said first and second voxel sets having the same number of voxels;

computing one or more statistical properties of said second voxel set;

computing an alternative result for that pixel location using said one or more statistical properties of said second voxel set; and

if said alternative result indicates an ultrasound echo intensity less than that indicated by said output value, resetting said output value to said alternative result.

8. The method of claim 7, said volumetric ultrasound representation being formed from scans taken during a sweep of an ultrasound probe across the breast, wherein said second slab-like subvolume is closer to a locus of said ultrasound probe sweep than said first slab-like subvolume.