[go: up one dir, main page]

EP3707672A1 - Méthode de segmentation d'un objet tridimensionnel dans une image de rayonnement médical - Google Patents

Méthode de segmentation d'un objet tridimensionnel dans une image de rayonnement médical

Info

Publication number
EP3707672A1
EP3707672A1 EP18796079.4A EP18796079A EP3707672A1 EP 3707672 A1 EP3707672 A1 EP 3707672A1 EP 18796079 A EP18796079 A EP 18796079A EP 3707672 A1 EP3707672 A1 EP 3707672A1
Authority
EP
European Patent Office
Prior art keywords
voxel
distance map
distance
slice images
slice
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP18796079.4A
Other languages
German (de)
English (en)
Inventor
Asma Ouji
Yoni De Witte
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Agfa HealthCare NV
Original Assignee
Agfa HealthCare NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Agfa HealthCare NV filed Critical Agfa HealthCare NV
Publication of EP3707672A1 publication Critical patent/EP3707672A1/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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 OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/174—Segmentation; Edge detection involving the use of two or more images
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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/10072—Tomographic images
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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 OR CALCULATING; 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

Definitions

  • the present invention relates to a method of segmenting a 3D object in a medical radiation image such as a Computed Tomography (CT) or a Magnetic Resonance image .
  • a medical radiation image such as a Computed Tomography (CT) or a Magnetic Resonance image .
  • Typical objects that need to be segmented in an image are multi-tissue organs (e.g. kidney, liver) and low contrasted tumors (e.g. brain tumor , liver tumor , etc.).
  • multi-tissue organs e.g. kidney, liver
  • low contrasted tumors e.g. brain tumor , liver tumor , etc.
  • the invention is valuable for diagnosis, medical assessment and follow up. It is indeed important to segment a tumor or an organ to measure its size and compare it with previous measurements over time . It is also important to segment such objects to better visualize and analyse their shape and morphological aspects .
  • a 3D mask is created via geometrical linear interpolations of contour distance maps defined in the 2D slice images .
  • a plane may be a slice of a tomographic representation of an object or it may be a modified slice, e.g. a rotated slice or any other 2D image representation.
  • contours are drawn by the user around the region to be segmented in a few planes, at least two planes being required (e.g. in Figure l . a) .
  • Contour points of the area to be segmented or a definition of the area defining the region to be segmented can be determined in other ways .
  • the algorithm creates a 3D segmentation mask by
  • This algorithm is designed to interpolate contours defined in parallel planes, but it can also handle non-parallel contours.
  • a 3D segmentation mask is computed as follows: a. A 2D distance map is computed in each plane where a contour is defined, b. A 3D distance map is created via a linear interpolation of the 2D distance maps, c. Each voxel is classified as in/out the segmentation mask depending on its corresponding distance map value.
  • the present invention is generally implemented in the form of a computer program product adapted to carry out the method steps of the present invention when run on a computer combined with user interaction to define the initial contours in the image planes.
  • the computer program product is commonly stored in a computer readable carrier medium such as a DVD.
  • the computer program product takes the form of an electric signal and can be communicated to a user through electronic communication.
  • Fig. 1 illustrates the interpolation of 3 parallel contours defining a kidney
  • Fig. 2 is a slice image of the interpolated mask of figure 1
  • Fig . 3 illustrates the interpolation of contours defining a
  • Fig. 4 is an example of a 2D distance map computed on a plane containing one contour
  • Fig. 5 is a bounding box of two parallel contours
  • Fig. 6 is an illustration of the 3D distance map interpolation .
  • Radiation images of 3D objects are typically generated by applying image recording techniques such as Computer Tomography or Magnetic Resonance Imaging .
  • the medical image representation generated by these techniques consists of a number of 2D slice images (also called 'planes' ) obtained by scanning the obj ect .
  • the slice images are generated by exposing the object and recording images from different angles so as to produce cross- sectional (tomographic) images (virtual "slices" ⁇ of specific areas of a scanned object .
  • CT imaging the results form an exposure are also slice images.
  • the radiation used for imaging can thus be of different types such as x-rays, radio waves etc.
  • the invention is generally applicable to 3D imaging techniques that produce slice images or image planes .
  • the input to the method of the present invention is a set of slice images, at least two, of a 3D image representation of an obj ect .
  • the slice images may be original slice images but likewise it may be modified slice images, e.g. rotated slice images .
  • these slice images not only comprise the object but also include pixels surrounding the object .
  • a segmentation process to separate the object from the remainder of the image is desired for some applications such as diagnosis , medical assessment and follow up.
  • the user performs a first and single action on the display of a number of slice images out of the acquired image representation of the 3D image .
  • This step is the only interactive one (the only step which requires user interaction) .
  • the results of the user action are fed into an image processing device that runs a software implementation of the method of the present invention as explained further on.
  • the slice images that will be taken into account are displayed on a monitor . At least two slice images are required to be able to perform the method of the present invention.
  • contours Figure . 1. a and Figure 3. a
  • the contours are drawn on the displayed 2D images (planes) . These planes can have any space direction but they should be preferably parallel .
  • the segmentation mask is computed inside the region defined by the bounding box of all contours .
  • a bounding box in the context of this invention is a volume that contains all defined contours .
  • this is the smallest bounding box containing all
  • a 2D distance map is computed in the region limited by the projection of the bounding box on this plane .
  • This distance map is defined such as all points inside a plane' s contour are positive and all remaining points are negative .
  • the absolute value of the 2D distance map at a given point is linearly dependent on its distance to the closest contour point ; i.e. for each point in the plane, the distance map absolute value gets bigger as its distance from the closest contour (defined in the same plane) is bigger .
  • Figure 4 illustrates the distance map computation algorithm.
  • the 3D distance map is defined in the region delimited by the bounding box containing all contours (example of bounding bo in Figure 5) .
  • the two closest planes are fetched. These two planes must surround (be located in opposite sides towards the voxel) the voxel if they are parallel .
  • dl be the distance of the voxel to the first closest lane and d2 its distance to the second one .
  • Total is the sum of dl and d2.
  • the voxel is projected on the two planes respectively.
  • costl and cost2 be the 2D distance map value of the projection point on the first and second planes respectively .
  • the interpolated distance map value is: costl x (1 - dl / total) + cost2 x (1 - d2 / total) .
  • the final segmentation mask is obtained by thresholding this distance map: all voxels with a positive distance map value belong the interpolation mask, and all others are outside of it.
  • voxels classified as being in or out the segmented area can be used for further processing, analysis, display etc.

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

Selon la présente invention, sur la base d'une entrée d'un utilisateur, un ensemble de points de contour de l'objet 3D est défini dans un certain nombre d'images de tranches en 2D représentant l'objet en 3D. Une carte de distance en 2D est calculée dans chaque plan où un contour est défini. Ensuite, une carte de distance en 3D est créée par l'intermédiaire d'une interpolation linéaire des cartes de distance en 2D. Chaque voxel est classé comme étant à l'intérieur ou à l'extérieur du masque de segmentation en fonction de sa valeur de carte de distance correspondante.
EP18796079.4A 2017-11-09 2018-11-09 Méthode de segmentation d'un objet tridimensionnel dans une image de rayonnement médical Withdrawn EP3707672A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP17200804 2017-11-09
PCT/EP2018/080725 WO2019092167A1 (fr) 2017-11-09 2018-11-09 Méthode de segmentation d'un objet tridimensionnel dans une image de rayonnement médical

Publications (1)

Publication Number Publication Date
EP3707672A1 true EP3707672A1 (fr) 2020-09-16

Family

ID=60293873

Family Applications (1)

Application Number Title Priority Date Filing Date
EP18796079.4A Withdrawn EP3707672A1 (fr) 2017-11-09 2018-11-09 Méthode de segmentation d'un objet tridimensionnel dans une image de rayonnement médical

Country Status (4)

Country Link
US (1) US20200286240A1 (fr)
EP (1) EP3707672A1 (fr)
CN (1) CN111295690A (fr)
WO (1) WO2019092167A1 (fr)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL2019905B1 (en) * 2017-11-14 2019-05-20 Mat Systems and methods for segmenting images
CN111986138B (zh) * 2019-05-22 2023-08-25 杭州依图医疗技术有限公司 一种获取肋骨定位的方法和装置
EP3832596A1 (fr) * 2019-12-06 2021-06-09 Microsoft Technology Licensing, LLC Segmentation d'images 3d
CN114202719A (zh) * 2021-11-12 2022-03-18 中原动力智能机器人有限公司 视频样本的标注方法、装置、计算机设备及存储介质

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL80364A (en) * 1986-10-20 1990-03-19 Elscint Ltd Three dimensional image construction using binary space interpolation
US6606091B2 (en) * 2000-02-07 2003-08-12 Siemens Corporate Research, Inc. System for interactive 3D object extraction from slice-based medical images
WO2004114063A2 (fr) * 2003-06-13 2004-12-29 Georgia Tech Research Corporation Reconstruction des donnees a l'aide de techniques d'interpolation directionnelle
CN1301494C (zh) * 2004-06-07 2007-02-21 东软飞利浦医疗设备系统有限责任公司 一种医学图像的三维分割方法
US8571277B2 (en) 2007-10-18 2013-10-29 Eigen, Llc Image interpolation for medical imaging
US8253739B2 (en) * 2008-04-03 2012-08-28 Siemens Aktiengesellschaft Method for interpolating an intermediate polygon p from two polygons p1 and p2
US9129391B2 (en) * 2011-09-28 2015-09-08 Siemens Aktiengesellschaft Semi-automated preoperative resection planning
CN105659289B (zh) 2013-10-30 2019-04-12 爱克发医疗保健公司 血管分割方法

Also Published As

Publication number Publication date
WO2019092167A1 (fr) 2019-05-16
US20200286240A1 (en) 2020-09-10
CN111295690A (zh) 2020-06-16

Similar Documents

Publication Publication Date Title
US12094123B2 (en) Image data processing apparatus and method
EP2916738B1 (fr) Systèmes et procédés d'imagerie de poumon, de lobe et de fissure
EP2144199B1 (fr) Procédé et appareil pour la visualisation multimodale d'ensembles de données de volume
EP2877980B1 (fr) Procédé et système de calcul d'erreur de déformation de dose
CN102292744B (zh) 心脏图像处理和分析
EP3707672A1 (fr) Méthode de segmentation d'un objet tridimensionnel dans une image de rayonnement médical
EP3555852B1 (fr) Procédé de segmentation d'un objet 3d dans une image de rayonnement médical
EP2476102A2 (fr) Améliorations apportées à une reconstruction planaire curviligne
CN102132322A (zh) 用于确定对象的尺寸改变的设备
EP3389006B1 (fr) Déploiement de nervure à partir d'images à résonance magnétique
EP4505398B1 (fr) Détection automatique de fracture de côtes à partir d'images d'acquisition dépliées
JP4991697B2 (ja) データセットにおいて構造体を分割する方法、システム及びコンピュータプログラム
CN102498501B (zh) 用于确定模糊图像中模糊的性质的系统和方法
JP4473578B2 (ja) 体の構造の孤立した視覚化されたものを形成する方法及び装置
Uemura et al. Generation of fractal dimension images and its application to automatic edge detection in brain MRI
Dong et al. Multiresolution cube propagation for 3-D ultrasound image reconstruction
Tran et al. An improved method for building a 3D model from 2D DICOM
EP3109824B1 (fr) Système et méthode pour la manipulation des données d'image
Göcke et al. Fast volume rendering methods for voxel-based 2D/3D registration-A comparative study
Liu et al. Study and application of medical image visualization technology
Patra et al. Medical Image Processing in Nuclear Medicine and Bone Arthroplasty
EP4322111A1 (fr) Procédé d'épluchage volumétrique pour traçage de chemin monte-carlo
US11443476B2 (en) Image data processing method and apparatus
EP3229210B1 (fr) Procédé et dispositif de traitement de données d'image médicale permettant de déterminer un modèle géométrique d'un objet de volume anatomique
JPH0236305A (ja) 被検体対象領域体積計測装置

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20200609

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: AGFA HEALTHCARE NV

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20210112