Classifications

• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/48—Diagnostic techniques
  • A61B6/482—Diagnostic techniques involving multiple energy imaging
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/46—Arrangements for interfacing with the operator or the patient
  • A61B6/461—Displaying means of special interest
  • A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
  • A61B6/502—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
  • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
  • A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
  • A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/54—Control of apparatus or devices for radiation diagnosis
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T11/00—2D [Two Dimensional] image generation
  • G06T11/60—Editing figures and text; Combining figures or text
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T5/00—Image enhancement or restoration
  • G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T7/00—Image analysis
  • G06T7/0002—Inspection of images, e.g. flaw detection
  • G06T7/0012—Biomedical image inspection
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T7/00—Image analysis
  • G06T7/10—Segmentation; Edge detection
  • G06T7/136—Segmentation; Edge detection involving thresholding
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
  • A61B6/46—Arrangements for interfacing with the operator or the patient
  • A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2207/00—Indexing scheme for image analysis or image enhancement
  • G06T2207/10—Image acquisition modality
  • G06T2207/10116—X-ray image
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2207/00—Indexing scheme for image analysis or image enhancement
  • G06T2207/20—Special algorithmic details
  • G06T2207/20092—Interactive image processing based on input by user
  • G06T2207/20104—Interactive definition of region of interest [ROI]
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2207/00—Indexing scheme for image analysis or image enhancement
  • G06T2207/30—Subject of image; Context of image processing
  • G06T2207/30004—Biomedical image processing
  • G06T2207/30068—Mammography; Breast
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2210/00—Indexing scheme for image generation or computer graphics
  • G06T2210/41—Medical
• • G—PHYSICS
  • G06—COMPUTING; CALCULATING OR COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2211/00—Image generation
  • G06T2211/40—Computed tomography
  • G06T2211/408—Dual energy

Definitions

• the invention relates to a method, a computer program, a computer-readable medium and a controller for processing a digital image comprising pixels with intensities relating to different energy levels of record. Additionally, the invention relates to an X-ray imaging system.
• transmission images of a region of interest of a patient may be recorded by detecting X-rays that fall onto a sensor or detector after having passed the region of interest. After that these images may be displayed and a physician or similarly skilled person may decide, whether there are malignant or benign changes in the region of interest.
• An X-ray imaging system may comprise an energy discriminating detector, i.e. a detector that is adapted for differentiating between X-rays of different energy and of recording images at different energy levels.
• an energy discriminating detector i.e. a detector that is adapted for differentiating between X-rays of different energy and of recording images at different energy levels.
• a dual energy detector may produce two energy images of exactly the same geometry and anatomical area.
• US 8 165 379 B2 shows a mammography system that is adapted to record low energy and high energy images.
• Two complementing energy images may show information which may help to discriminate whether visible masses are likely malignant or benign. However, the way to compute malignancy features from the dual energy images is not yet clear.
• An aspect of the invention relates to a method for processing a digital image, which comprises pixels with intensities relating to different energy levels. I.e. each pixel of the image may be associated with at least two intensities, which have been recorded at different energy levels.
• the digital image may be a digital X-ray image and the energies may relate to different X-ray energies during the recording of the digital X-ray image.
• the method comprises the steps of: receiving first image data and second image data of the digital image, the first image data encoding a first energy level and the second image data encoding a second energy level; determining a regression model from the first image data and the second image data, the regression model establishing a correlation between intensities of pixels of the first image data with intensities of pixels of the second image data; and calculating residual mode image data from the first image data and the second image data, such that a pixel of the residual mode image data has an intensity based on the difference of an intensity of the second image data at the pixel and a correlated intensity of the pixel of the first image data, the correlated intensity determinate by applying the regression model to the intensity of pixel of the first image data.
• the basic idea of the invention may be seen that two energy images, i.e. a first energy image and a second energy image (represented by the first (energy) image data and the second (energy) image data) are recorded and evaluated such that new information is extracted from the two energy images that may not be directly seen by simply comparing the two energy images. This new information may be visualized with the residual mode image data.
• a correlation between the two energy images may be analyzed, and a regression model may be established based on the correlation.
• the regression model may be used for predicting an intensity of each pixel in the second energy image from the intensity of the same pixel in the first energy image.
• new images may be computed with the aid of the regression model.
• the residual mode image (represented by residual mode image data) may display information in the second image that could not be predicted from the first image, i.e. new information.
• a dominant mode image (represented by dominant mode image data) may display highly correlated components in the two energy images.
• a computation of malignancy features from dual energy images may be improved.
• the residual mode image it may be established which parts of a second energy image significantly add information to a first energy image.
• a computer-readable medium may be a floppy disk, a hard disk, an USB (Universal Serial Bus) storage device, a RAM (Random Access Memory), a ROM (Read Only Memory), an EPROM (Erasable Programmable Read Only Memory) or a FLASH memory.
• a controller may comprise a processor for executing the computer program, which may be stored in a memory of the controller.
• an X-ray imaging system which comprises a detector arrangement for recording first image data with X-rays of a first energy level and for recording second image data with X-rays of a second energy level, a controller, which is adapted for generating residual mode image data from the first image data and the second image data, and a display device for displaying the residual mode image data.
• the X-ray imaging system may be a mammography station.
• the image data is recorded at a first place and processed at a second place, for example at a workstation remotely connected to an X-ray device.
• Fig. 1 schematically shows an X-ray imaging system according to an embodiment of the invention.
• Fig. 2 shows a flow diagram for a method for processing image data according to an embodiment of the invention.
• Fig. 3 schematically shows a display device for an X-ray imaging system according to an embodiment of the invention.
• Fig. 4 shows an example for a correlation of energy intensities of two images.
• Fig. 5 shows an example of a dominant mode image and a residual mode image generated with a method for processing image data according to an embodiment of the invention.
• Fig. 6 shows an example of a dominant mode image and a residual mode image generated with a method for processing image data according to an embodiment of the invention.
• Fig. 1 shows an X-ray imaging system 10 that comprises a detector arrangement 12, a controller 14 and a display device 16.
• the detector arrangement 12 which may be controlled by the controller 14, is adapted for recording digital images at two different energy levels of an object of interest 18, such as a breast.
• the detector arrangement 12 may comprise an energy discriminating detector 20 or an X-ray source 22 that is adapted for generating X-ray radiation at different energy levels at different times.
• Fig. 2 shows a flow diagram for a method for processing a digital image 40 that may be recorded by the detector arrangement 12.
• the method may be executed by the controller 14, which, for example, may comprise a processor running a corresponding computer program.
• step 30 the digital image 40 is recorded and received in the controller 14.
• the controller 14 may control the detector arrangement 12 in such a way that digital image data 42a, 42b is recorded at different energy levels, for example with an energy discriminating detector 20.
• the image 40 and its parts 42a, 42b may be recorded simultaneously or during a short time period. In such a way, the image 40 may show the same geometry of the region of interest 18 and/or the same anatomical area 18.
• the method comprises the step of: recording the first image data 42a and the second image data 42b with an X-ray detector arrangement 12 adapted for acquiring X-rays at different X-ray energy levels.
• Each of the first and second image data 42a, 42b (or energy image data 42a, 42b) comprises pixels associated with intensities relating to the intensity of the respective energy level at the respective pixel.
• the first and second image data 42a, 42b may have the same size in pixels (as representing the same field of view on the region of interest). It has to be noted that the digital image 40 may comprise more than two sets of image data 42a, 42b, associated with more than two energy levels.
• the method comprises the step of: receiving first image data 42a and second image data 42b of the digital image 40, the first image data 42a encoding a first energy level and the second image data 42b encoding a second energy level.
• step 32 a regression model is determined from the first and second image data 42a, 42b.
• Fig. 4 shows an example of a pixel-wise scatter plot of the intensities of two sets of energy image data 42a, 42b, which show a high but non-linear correlation. This correlation between the two sets of energy image data is analyzed, to generate a regression model 44.
• the regression model 44 is a mapping that maps an intensity of a pixel of the first image data 42a to an intensity of the second image data 42b.
• an intensity of each pixel in the second image data 42b may be predicted from the corresponding intensity in the first image data 42a.
• the method comprises the step of: determining a regression model 44 from the first image data 42a and the second image data 42b, the regression model 44 establishing a correlation between intensities of pixels of the first image data 42a with intensities of pixels of the second image data 42b.
• the regression model 44 may be derived either for the whole image area (i.e. all pixels of the image data 42a, 42b) or only for a selected region 50, which may be selected by a user, as will be explained below.
• the regression model 44 is determined from pixels of only a selected region 50.
• the correlation between the two sets of image data 42a, 42b may be non-linear. Therefore, it is possible that a non-linear regression model 44 is established, for example a piece-wise linear model or a Support- Vector- Regression model.
• the regression model 44 is a non-linear model.
• the regression model may also be linear and may be based, for example, on linear decorrelation techniques, such as Principal Component Analysis.
• residual mode image data 46 and/or dominant mode image data 48 are determined from the first and second image data 42a, 42b with the aid of the regression model 44.
• the first and second image data 42a, 42b may be compared and analyzed.
• intensities of pixels of the one set of image data 42a may be mapped with the regression model 44 to intensities that are comparable with intensities of corresponding pixels of the other set of image data 42b.
• a corresponding pixel may be a pixel at the same position, which may have the same coordinates.
• a residual mode image (represented by residual mode image data 46) may be computed containing the residuals of the correlation, i.e. the new information added by the second image data 42b, which cannot have been predicted from the first image data 42a with the aid of the regression model 44.
• the method comprises the step of: calculating residual mode image data 46 from the first image data 42a and the second image data 42b, such that a pixel of the residual mode image data 46 has an intensity based on the difference of an intensity of the second image data 42b at the pixel and a correlated intensity of the pixel of the first image data 42a, the correlated intensity determinate by applying the regression model to the intensity of pixel of the first image data 42a.
• the residual mode image data 42a may be used to visualize which parts of the second image data 42b are truly adding information to the first image data 42a, i.e. are non- redundant and cannot be simply predicted from the first one.
• the parts of the first and second image data 42a, 42b that are predictable from each other may be computed.
• a dominant mode image (represented by dominant mode image data 48) may be computed, containing the dominant mode, i.e. the highly correlated component between the two energy images.
• the method further comprises the step of: calculating dominant mode image data 48, wherein an intensity of a pixel of the dominant mode image data 48 is based on a correlated component of the pixel of the first image data 42a and second image data 42b with respect to the regression model 44.
• Fig. 5 and 6 show examples of dominant mode images 48 and residual mode images 46 that are depicted side by side and that may be displayed in the shown way on a display device 16 as will be explained below.
• the residual mode image data 46 may further be analyzed to find an intensity threshold above which the intensity difference can be considered meaningful, for example to eliminate noise.
• the residual mode image data 46 may contain mainly noise.
• An intensity threshold may be set, above which the deviation from the dominant mode image data 48 can be considered meaningful.
• the such established threshold intensity of the residual mode image data 46 may be used to display only pixels which are above a specific noise level.
• the method further comprises the step of: applying a threshold intensity to the residual mode image data 46, such that pixels of the residual mode image data 46 with an intensity below the threshold intensity are discarded.
• the noise level of the residual mode image data 46 may be analyzed by means of a curve of Euler characteristics, i.e. for each possible threshold the number of holes and blobs is counted, to build an Euler histogram.
• Techniques are known for establishing the Euler characteristic simultaneously for possible thresholds in a single pass over the residual mode image data 46.
• a threshold intensity may be derived from the curve position at which the Euler characteristics drops below a predefined number, indicating larger spatial structures in the residual mode image data 46.
• the threshold intensity is determined such that an Euler characteristic of the residual mode image data 46 drops below a predefined number at the threshold intensity.
• step 38 the residual mode image data 46 and/or the dominant mode image data 48 are displayed on the display device 16.
• Fig. 3 shows an example of a screen content of the display device 16, which may be a CRT or LCD or similar device.
• Each set of image data 42a, 42b, 46, 48 may be displayed on the display device.
• the residual mode image 46 may be presented separately from the dominant mode image 48 as an additional image.
• the residual mode image 46 may be presented as an overlay or alpha-blended to the dominant mode image 48 (or the other energy images 42a, 42b), either for the whole image area or only for a selected region 50.
• the residual mode image 46 is presented an overlay to the dominant mode image 48, for example as a color overlay, wherein the opacity is determined by the intensities of the residual mode image 46.
• the degree of color overlay may be manually controlled by a user, for example with a slider 54 or with a mouse wheel.
• the method further comprises the step of: displaying the residual mode image data 46 together with further image data 42a, 42b, 48 on a display device 16 by overlaying the image data with the further image data 46.
• the further image data may be one of the first image data 42a, second image data 42b, and the dominant image data 48.
• the residual mode image 46 is toggled in place with the further image data 42a, 42b, for example with the aid of a toggle button 52.
• a user may manually control the toggling of the images.
• the method further comprises the step of: displaying the residual mode image data 46 together with further image data 42a, 42b, 48 on a display device 16 by toggling the image data with the residual mode image data 46 by a user command on the display device 16.
• the residual mode image 36 may be gradually alpha- blended in gray-scale and in place with the other image data 42a, 42b, 48.
• the degree of alpha-blending or the binary toggling may be manually controlled by a user, for example with a slider 54 or mouse wheel.
• the method further comprises the step of: displaying the residual mode image data 46 together with further image data 42a, 42b, 48 on a display device 16 by alpha-blending the image data with the residual mode image data 46 by a user command.
• a region 50 of one of the images 42a, 42b, 46, 48 may be selected by a user. For example, a rectangle may be selected with a mouse.
• the selected region 50 may be used for defining a region in the digital image 40 from which the regression model 44 is determined, i.e. the regression model 44 may be determined from the whole image 40 or only from a part of the image 40.
• the selected region 50 also may be used for defining a region in one of the images 42a, 42b, 48, to which the residual mode image data 46 is overlaid or alpha-blended. According to an embodiment of the invention, the residual mode image data 46 is displayed only for pixels of only the selected region 50.
• two images may be linked such that a mouse click or movement in either one of the two images shows a hair cross at the corresponding position in the other image.
• the user may freely move with the selected region over the image 42a, 42b, 46, 48 (like a magic magnifier).
• the method may be applied to the multiple energy image 40 as a whole, or to a local region of interest, interactively steered by a user like a magic magnifying glass.

Landscapes

• Engineering & Computer Science (AREA)
• Health & Medical Sciences (AREA)
• Life Sciences & Earth Sciences (AREA)
• Medical Informatics (AREA)
• Physics & Mathematics (AREA)
• General Health & Medical Sciences (AREA)
• Radiology & Medical Imaging (AREA)
• Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
• Biomedical Technology (AREA)
• Veterinary Medicine (AREA)
• Biophysics (AREA)
• Optics & Photonics (AREA)
• Pathology (AREA)
• High Energy & Nuclear Physics (AREA)
• Public Health (AREA)
• Heart & Thoracic Surgery (AREA)
• Molecular Biology (AREA)
• Surgery (AREA)
• Animal Behavior & Ethology (AREA)
• General Physics & Mathematics (AREA)
• Theoretical Computer Science (AREA)
• Computer Vision & Pattern Recognition (AREA)
• Human Computer Interaction (AREA)
• Dentistry (AREA)
• Oral & Maxillofacial Surgery (AREA)
• Quality & Reliability (AREA)
• Image Processing (AREA)
• Apparatus For Radiation Diagnosis (AREA)

Priority Applications (4)

Applications Claiming Priority (2)

Publications (2)

Family

ID=50549361

Family Applications (1)

Country Status (5)

Cited By (3)

Families Citing this family (3)

Citations (1)

Family Cites Families (14)

• 2014
  • 2014-03-14 CN CN201480015630.4A patent/CN105122300B/zh not_active Expired - Fee Related
  • 2014-03-14 WO PCT/IB2014/059770 patent/WO2014141163A2/en active Application Filing
  • 2014-03-14 US US14/774,205 patent/US9706968B2/en active Active
  • 2014-03-14 JP JP2015562518A patent/JP6251297B2/ja not_active Expired - Fee Related
  • 2014-03-14 EP EP14719356.9A patent/EP2973391B1/en not_active Not-in-force

Patent Citations (1)

Also Published As

Similar Documents

Legal Events