FR2790123A1 - Procede de traitement d'une sequence d'images fluoroscopiques d'un corps, de facon a ameliorer la qualite des images visualisees - Google Patents

Procede de traitement d'une sequence d'images fluoroscopiques d'un corps, de facon a ameliorer la qualite des images visualisees Download PDF

Info

Publication number
FR2790123A1
FR2790123A1 FR9902032A FR9902032A FR2790123A1 FR 2790123 A1 FR2790123 A1 FR 2790123A1 FR 9902032 A FR9902032 A FR 9902032A FR 9902032 A FR9902032 A FR 9902032A FR 2790123 A1 FR2790123 A1 FR 2790123A1
Authority
FR
France
Prior art keywords
image
images
current
acquired
displacement
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.)
Granted
Application number
FR9902032A
Other languages
English (en)
Other versions
FR2790123B1 (fr
Inventor
Vladislav Boutenko
Thierry Lebihen
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.)
GE Medical Systems SCS
Original Assignee
GE Medical Systems SCS
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 GE Medical Systems SCS filed Critical GE Medical Systems SCS
Priority to FR9902032A priority Critical patent/FR2790123B1/fr
Priority to US09/503,022 priority patent/US7106894B1/en
Priority to EP00301211A priority patent/EP1030268A1/fr
Priority to JP2000037619A priority patent/JP2000244817A/ja
Publication of FR2790123A1 publication Critical patent/FR2790123A1/fr
Application granted granted Critical
Publication of FR2790123B1 publication Critical patent/FR2790123B1/fr
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

Pour chaque image courante acquise (Xn+1 ), on détermine le déplacement (D) de ladite image courante par rapport à l'image précédente acquise (Xn ) dans le plan d'acquisition des images, on élabore une image filtrée précédente décalée, en décalant spatialement l'image filtrée précédente compte tenu dudit déplacement, et on élabore l'image filtrée courante par la moyenne pondérée entre l'image courante acquise (Xn+1 ) et l'image filtrée précédente décalée, de façon à améliorer la qualité des images visualisées.Application à la fluoroscopie.

Description

Procédé de traitement d'une séquence d'images fluoroscopiques
d'un corps, de façon à améliorer la qualité des images visualisées.
L'invention concerne la fluoroscopie, et plus particulièrement le traitement d'une séquence d'images fluoroscopiques d'un corps, en
particulier un corps humain.
L'invention s'applique par ailleurs plus particulièrement à la fluoroscopie cardiaque. Par rapport à la radiographie en mode acquisition o les doses de rayons X sont plus importantes dans le but d'obtenir des images enregistrées de meilleure qualité à des fins de diagnostics, la fluoroscopie s'effectue avec des doses de rayons X plus faibles et est plus particulièrement utilisée dans le domaine interventionnel pour positionner par exemple des endoprothèses coronaires ("stent" en langue
anglaise) à l'aide de cathéters.
En fluoroscopie, le mouvement des objets d'intérêt, comme par exemple des endoprothèses coronaires, ainsi que, notamment en fluoroscopie cardiaque, le mouvement du fond lié par exemple à la respiration du patient, ainsi qu'aux mouvements de la table sur laquelle est placé le patient, provoquent des perturbations dans les images auxquelles s'ajoute le bruit en particulier d'origine électrique. Ce bruit est d'autant plus gênant que l'examen fluoroscopique dure longtemps, typiquement de l'ordre de 5 minutes pour positionner correctement une prothèse intravasculaire, et occasionne par conséquent une gêne visuelle pour le médecin. Si l'on était en présence d'images immobiles, on saurait aisément éliminer le bruit par un simple filtrage temporel. Cependant, en fluoroscopie, on est en présence d'images mobiles qui se traduisent, si l'on effectue simplement un filtrage des images (par exemple un filtrage par une moyenne temporelle) par un mouvement flou ou bien une perte de contraste des objets mobiles (selon la taille des objets). En d'autres termes, on ne fait alors pas la différence entre l'arrivée et le départ d'un objet d'intérêt et un pic de bruit. Actuellement, les algorithmes classiques de traitement d'images en fluoroscopie font appel à un critère de distinction entre une variation due au bruit et une variation due au mouvement. On stoppe ou diminue alors le traitement de filtrage en présence d'un mouvement. Cependant, l'arrêt du filtrage provoque la réapparition du bruit, ce qui se traduit sur
les images par des traînées de bruits derrière les objets mobiles.
L'invention vise à apporter une solution plus satisfaisante à ces problèmes. L'invention propose donc un procédé de traitement d'une séquence d'images fluoroscopiques d'un corps, comprenant l'acquisition d'une séquence d'images, l'élaboration, pour chaque image courante acquise, d'une image filtrée courante à partir de l'image courante acquise et de l'image filtrée précédente, et la visualisation de la séquence d'images filtrées. Selon une caractéristique générale de l'invention, pour chaque image courante acquise, on détermine le déplacement de ladite image courante par rapport à l'image précédente acquise dans le plan d'acquisition des images, on élabore une image filtrée précédente dite "décalée", en décalant spatialement l'image filtrée précédente compte tenu dudit déplacement, et on élabore l'image filtrée courante par la moyenne pondérée entre l'image courante acquise et l'image filtrée
précédente décalée, de façon à améliorer la qualité des images visualisées.
Lorsque le corps est disposé sur une table déplaçable, on détermine avantageusement le déplacement de ladite image courante dans le plan d'acquisition des images à partir de la valeur du déplacement de la table, de l'orientation spatiale du plan d'acquisition par rapport à la table
et de la distance de ce plan d'acquisition par rapport à la table.
On peut également en variante, déterminer le déplacement de ladite image courante dans le plan d'acquisition des images à partir du
contenu de ces images acquises.
D'autres avantages et caractéristiques de l'invention
apparaîtront à l'examen de la description détaillée d'un mode de mise en
oeuvre, nullement limitatif, et des dessins annexés, sur lesquels: - la figure 1 illustre schématiquement un dispositif permettant une mise en oeuvre du procédé selon l'invention; et - la figure 2 illustre schématiquement deux images successives acquises. Sur la figure 1, la référence TB désigne une table déplaçable en translation selon deux directions orthogonales Fx et Fy, à l'aide par
exemple d'une manivelle non représentée ici à des fins de simplification.
Des capteurs de déplacement CDx et CDy, de réalisation connue en soi, permettent de déterminer les valeurs des déplacements selon les directions Fx et Fy, respectivement, et délivrent ces informations à des
moyens de traitement MT comportant un microprocesseur.
Sur la table TB, un patient CP est allongé et subit un examen de fluoroscopie. A cet égard, l'appareillage de fluoroscopie comporte une source S de rayons X émettant des rayons X selon un axe Ax en direction d'un détecteur DC de réalisation classique et connue en soi, relié également aux moyens de traitement MT. Le détecteur DC est un détecteur plan, formant le plan d'acquisition des images. Ce plan est orthogonal à l'axe Ax. Le détecteur DC et la source S font partie d'un appareillage
comportant un bras mobile dans l'espace autour du patient CP.
L'orientation de l'axe Ax est parfaitement connue à chaque instant, de même que la distance entre la table TB et le détecteur DC. Les déplacements selon les directions Fx et Fy de la table TB peuvent donc aisément être convertis de façon classique et connue en soi en des déplacements u et v dans le plan du détecteur DC, c'est-à-dire dans le plan
d'acquisition des images.
On se réfèere maintenant plus particulièrement à la figure 2. On voit que lors de l'examen de fluoroscopie, une séquence d'images Xn sont acquises, typiquement au rythme de 30 images par seconde. Les valeurs des pixels de chaque image sont stockées, au fur et à mesure de leur acquisition dans une mémoire des moyens de traitement MT, afin de
permettre le traitement des images, et en particulier leur filtrage.
Ceci étant, lorsque, au cours de l'examen de fluoroscopie, le médecin déplace la table TB et/ou l'endoprothèse coronaire de façon à maintenir la région considérée du corps dans le champ du rayonnement X, l'objet d'intérêt (par exemple l'endoprothèse coronaire) que l'on a représenté très schématiquement sous la référence A dans l'image Xn, s'est déplacé du vecteur déplacement D dans l'image suivante Xn+1 (POet POD représentent respectivement les centres de la prothèse dans les deux images). Le vecteur déplacement D a des coordonnées u et v dans le plan d'acquisition des images, qui correspondent au déplacement de la table
dans les directions Fx et Fy (on suppose ici que seule la table a bougé).
L'invention prévoit ici d'effectuer un traitement de filtrage en utilisant d'une part l'image acquise courante Xn+l et d'autre part, non pas l'image filtrée précédente correspondant directement à l'image acquise
précédente, mais une image filtrée décalée YD.
Plus précisément, on implémente dans le microprocesseur la loi récursive (1): Yn+l = (1-a) Xn+l1 + a YDn (1) dans laquelle Yn+l représente l'image filtrée courante, Xn+l l'image acquise courante, YDn l'image filtrée précédente décalée par rapport à l'image filtrée précédente Yn, et "a" un coefficient de pondération,
typiquement égal à 0,2 (et éventuellement fonction de Xn+1 et de YDn).
Cette loi (1) se traduit pour chaque pixel de coordonnées i,j dans le plan d'acquisition des images, par la loi (2): Yn+l (ij) = (1-a) xn+ l1 (i,j) + a Yn (i-u, j-v) (2) En d'autres termes, pour chaque pixel de coordonnées i,j de l'image courante acquise, on calcule dans l'image filtrée précédente les coordonnées I,J du pixel décalé compte tenu du mouvement de la table
(I = i-u, J = hj-v) et on applique la loi (2).
Les images filtrées successives Yn+l sont successivement affichées sur les moyens de visualisation MV du dispositif de
fluoroscopie.
Lorsque les coordonnées I et J sont négatives, c'est-à-dire lorsqu'il n'y a pas d'équivalent dans l'image précédente pour une portion de l'image de sortie Yn+l, cette portion d'image est noircie. En d'autres
termes, on fixe alors les valeurs Yn+l (i,j) à zéro.
Ce noircissement n'est pas visible compte tenu de la fréquence d'acquisition des images et des mouvements habituellement rapides de déplacement de la table. En pratique, il a été observé que ce noircissement n'est pas visible lorsqu'il n'excède pas 10% de l'image pour le mouvement
de déplacement le plus rapide.
Dans le mode de mise en oeuvre qui vient d'être décrit, les
coordonnées u et v sont déterminées à partir du mouvement de la table.
Ceci étant, il est également possible en variante de calculer les coordonnées u et v du vecteur déplacement directement entre deux images successives acquises, par des algorithmes classiques de détection de mouvement en utilisant par exemple le critère de maximum de
corrélation entre deux environnements de deux pixels homologues.
Matériellement, le décalage des images filtrées s'effectue très simplement en utilisant des éléments de retard reliés aux lignes et aux
colonnes de la mémoire de stockage des images.

Claims (3)

REVENDICATIONS
1. Procédé de traitement d'une séquence d'images fluoroscopiques d'un corps, comprenant l'acquisition d'une séquence d'images, l'élaboration, pour chaque image courante acquise (Xn+l), d'une image filtrée courante (Yn+l) à partir de l'image courante acquise (Xn+l) et de l'image filtrée précédente, et la visualisation de la séquence d'images filtrées, caractérisé par le fait que pour chaque image courante acquise (Xn+l) on détermine le déplacement (t) de ladite image courante par rapport à l'image précédente acquise dans le plan d'acquisition des images, on élabore une image filtrée précédente décalée (YDn) en décalant spatialement l'image filtrée précédente (Yn) compte tenu dudit déplacement, et on élabore l'image filtrée courante (Yn+l) par la moyenne pondérée entre l'image courante acquise (Xn+l) et l'image filtrée précédente décalée (YDn), de façon à améliorer la qualité des images visualisées.
2. Procédé selon la revendication 1, dans lequel le corps est disposé sur une table déplaçable (TB), caractérisé par le fait qu'on détermine le déplacement de ladite image courante dans le plan d'acquisition des images (DC) à partir de la valeur du déplacement de la table (TB) et de l'orientation spatiale et de la distance du plan
d'acquisition par rapport à la table.
3. Procédé selon la revendication 1, caractérisé par le fait qu'on détermine le déplacement (D) de ladite image courante dans le plan
d'acquisition des images à partir du contenu des images acquises.
FR9902032A 1999-02-18 1999-02-18 Procede de traitement d'une sequence d'images fluoroscopiques d'un corps, de facon a ameliorer la qualite des images visualisees Expired - Fee Related FR2790123B1 (fr)

Priority Applications (4)

Application Number Priority Date Filing Date Title
FR9902032A FR2790123B1 (fr) 1999-02-18 1999-02-18 Procede de traitement d'une sequence d'images fluoroscopiques d'un corps, de facon a ameliorer la qualite des images visualisees
US09/503,022 US7106894B1 (en) 1999-02-18 2000-02-14 Method for reducing X-ray noise
EP00301211A EP1030268A1 (fr) 1999-02-18 2000-02-16 Méthode pour la réduction du bruit de rayons X
JP2000037619A JP2000244817A (ja) 1999-02-18 2000-02-16 X線ノイズを低減するための方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR9902032A FR2790123B1 (fr) 1999-02-18 1999-02-18 Procede de traitement d'une sequence d'images fluoroscopiques d'un corps, de facon a ameliorer la qualite des images visualisees

Publications (2)

Publication Number Publication Date
FR2790123A1 true FR2790123A1 (fr) 2000-08-25
FR2790123B1 FR2790123B1 (fr) 2001-04-27

Family

ID=9542235

Family Applications (1)

Application Number Title Priority Date Filing Date
FR9902032A Expired - Fee Related FR2790123B1 (fr) 1999-02-18 1999-02-18 Procede de traitement d'une sequence d'images fluoroscopiques d'un corps, de facon a ameliorer la qualite des images visualisees

Country Status (4)

Country Link
US (1) US7106894B1 (fr)
EP (1) EP1030268A1 (fr)
JP (1) JP2000244817A (fr)
FR (1) FR2790123B1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8121192B2 (en) 2005-04-11 2012-02-21 General Electric Company Method and system for treating noise in a generated image

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE0003608L (sv) * 2000-10-06 2001-12-17 Xcounter Ab Metod för reducering av högfrekvent brus i bilder med hjälp av medelvärdesbildning av pixlar och parvis addering av pixelpar som uppfyller ett villkor
DE10220295A1 (de) * 2002-05-07 2003-11-20 Philips Intellectual Property Verfahren zur Verbesserung der Bildqualität
FR2884342B1 (fr) * 2005-04-11 2008-01-04 Gen Electric Procede et systeme d'amelioration d'une image digitale generee a partir d'un detecteur de rayons x
JP4924384B2 (ja) * 2007-11-26 2012-04-25 株式会社島津製作所 放射線撮像装置および放射線検出信号処理方法
CN101744628B (zh) * 2008-11-28 2012-05-23 东软飞利浦医疗设备***有限责任公司 一种基于运动信号反馈的降噪装置及方法
JP5053982B2 (ja) * 2008-12-05 2012-10-24 株式会社東芝 X線診断装置および画像処理装置
JP6381198B2 (ja) * 2013-11-08 2018-08-29 キヤノン株式会社 制御装置、制御方法及びプログラム
KR102448069B1 (ko) * 2022-05-27 2022-09-27 주식회사 디알텍 방사선 촬영에 의해 획득된 영상에 포함된 노이즈를 저감하기 위한 영상 처리 방법 및 장치

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4542459A (en) * 1982-11-26 1985-09-17 General Electric Company Matched filter for x-ray hybrid subtraction
FR2560037B1 (fr) 1984-02-28 1988-04-08 Thomson Cgr Procede de controle de positionnement d'un patient par rapport a une installation de radiologie, et installation pour la mise en oeuvre de ce procede
US5309237A (en) 1992-03-31 1994-05-03 Siemens Corporate Research, Inc. Apparatus and method of compensating image-sequences for motion
US6205259B1 (en) * 1992-04-09 2001-03-20 Olympus Optical Co., Ltd. Image processing apparatus
DE69322444T2 (de) * 1992-07-10 1999-06-24 Koninklijke Philips Electronics N.V., Eindhoven Röntgendurchleuchtungsgerät mit Mitteln zur Rauschreduktion
DE4417628C1 (de) 1994-05-19 1995-09-28 Siemens Ag Verfahren zur adaptiven Rauschverminderung für digitale Bildsequenzen
FR2736181A1 (fr) * 1995-06-30 1997-01-03 Philips Electronique Lab Procede de traitement d'images pour la reduction du bruit dans une image d'une sequence d'images numeriques et dispositif mettant en oeuvre ce procede
US5646867A (en) * 1995-07-24 1997-07-08 Motorola Inc. Method and system for improved motion compensation
FR2751825A1 (fr) * 1996-07-24 1998-01-30 Philips Electronics Nv Procede de filtrage temporel du bruit dans une image d'une sequence d'images numerisees et dispositif mettant en oeuvre ce procede
EP0871143B1 (fr) * 1997-04-08 2004-06-02 Koninklijke Philips Electronics N.V. Système de traitement d'une séquence d'images bruitées, et appareil d'examen médical incluant un tel système
FR2772956A1 (fr) * 1997-12-23 1999-06-25 Philips Electronics Nv Procede de traitement d'image incluant des etapes de filtrage spatiotemporel du bruit et appareil d'imagerie medicale mettant en oeuvre ce procede
EP0926626A1 (fr) * 1997-12-23 1999-06-30 Koninklijke Philips Electronics N.V. Procédé de traitement d'image pour réduire le bruit dans une image d'une séquence de trois images temporelles et appareil d'imagerie médicale mettant en oeuvre un tel procédé
EP1046134B1 (fr) * 1998-10-19 2003-12-03 Koninklijke Philips Electronics N.V. Dispositif d'examen aux rayons x, a commande de dose de ceux-ci
EP1064790A1 (fr) * 1999-01-15 2001-01-03 Koninklijke Philips Electronics N.V. Codage et filtrage du bruit d'une sequence d'images
US6442617B1 (en) * 1999-03-31 2002-08-27 3Com Corporation Method and system for filtering multicast packets in a peripheral component environment
WO2001041648A1 (fr) * 1999-12-07 2001-06-14 Koninklijke Philips Electronics N.V. Procede de traitement d'images a ultrasons et systeme d'affichage d'une sequence d'images composites d'un segment arteriel
US7187810B2 (en) * 1999-12-15 2007-03-06 Medispectra, Inc. Methods and systems for correcting image misalignment
US6809758B1 (en) * 1999-12-29 2004-10-26 Eastman Kodak Company Automated stabilization method for digital image sequences
US6436041B1 (en) * 2000-12-22 2002-08-20 Acuson Corporation Medical ultrasonic imaging method with improved ultrasonic contrast agent specificity
US6623428B2 (en) * 2001-10-11 2003-09-23 Eastman Kodak Company Digital image sequence display system and method
JP2003151483A (ja) * 2001-11-19 2003-05-23 Hitachi Ltd 荷電粒子線を用いた回路パターン用基板検査装置および基板検査方法
US20040209237A1 (en) * 2003-04-18 2004-10-21 Medispectra, Inc. Methods and apparatus for characterization of tissue samples
US7499578B2 (en) * 2002-10-18 2009-03-03 Cornell Research Foundation, Inc. System, method and apparatus for small pulmonary nodule computer aided diagnosis from computed tomography scans
WO2004114063A2 (fr) * 2003-06-13 2004-12-29 Georgia Tech Research Corporation Reconstruction des donnees a l'aide de techniques d'interpolation directionnelle

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
AUFRICHTIG R ET AL: "X-RAY FLUOROSCOPY SPATIO-TEMPORAL FILTERING WITH OBJECT DETECTION", IEEE TRANSACTIONS ON MEDICAL IMAGING, vol. 14, no. 4, 1 December 1995 (1995-12-01), pages 733 - 746, XP000548591, ISSN: 0278-0062 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8121192B2 (en) 2005-04-11 2012-02-21 General Electric Company Method and system for treating noise in a generated image

Also Published As

Publication number Publication date
US7106894B1 (en) 2006-09-12
JP2000244817A (ja) 2000-09-08
EP1030268A1 (fr) 2000-08-23
FR2790123B1 (fr) 2001-04-27

Similar Documents

Publication Publication Date Title
EP1869643B1 (fr) Dispositif de traitement d'image et methode d'imagerie du debit sanguin
US11751848B2 (en) Methods and apparatuses for ultrasound data collection
US6088488A (en) Vascular imaging with adaptive averaging
WO2003045263A3 (fr) Systeme de visualisation d'images medicales et procede permettant d'ameliorer des structures dans des images bruyantes
CN111540025B (zh) 预测用于图像处理的图像
US20110235889A1 (en) Method for enhanced visualization of objects in interventional angiographic examinations and device for performing the method
CN101443815A (zh) 重建图像的方法和装置
FR2822273A1 (fr) Procede d'etalonnage pour la reconstruction de modelisations tri-dimensionnelles a partir d'images obtenues par tomographie
JP2007130244A (ja) X線アンギオ撮影装置
FR2790123A1 (fr) Procede de traitement d'une sequence d'images fluoroscopiques d'un corps, de facon a ameliorer la qualite des images visualisees
EP0936575A1 (fr) Procédé de traitement d'images pour l'estimation de mouvement dans une séquence d'images, procédé de filtrage du bruit et appareil d'imagerie médicale utilisant ces procédés
EP1408399A3 (fr) Système et méthode de traitement d'une image numérique pour un visionnement intuitif
WO2013067683A1 (fr) Procédé et système d'acquisition d'images pour le rendu d'images stéréoscopiques à partir d'images monoscopiques
US11436790B2 (en) Passthrough visualization
FR2813973A1 (fr) Procede et dispositif de generation d'images tridimensionnelles et appareil de radiologie associe
JP2013085955A (ja) 連続マルチスケール再構成において詳細画像を補うx線コンピュータ断層撮像装置(x線ct装置)、医用画像処理装置及び医用画像処理方法
US20120057671A1 (en) Data acquisition and visualization mode for low dose intervention guidance in computed tomography
US10504235B2 (en) Method for generating three dimensional images
FR2799029A1 (fr) Procede de reconstruction d'une image tridimensionnelle d'un objet en mouvement, en particulier une image tridimensionnelle de vaisseaux du coeur humain
EP3658031B1 (fr) Reconstruction de valvule cardiaque à compensation de mouvement
WO2017155825A1 (fr) Procédé de reconstruction multivue en 3d au moyen d'un suivi d'éléments et d'un enregistrement de modèle
Krishnan et al. Estimation of distance to texture surface using complex log mapping
CN106470609B (zh) 用于荧光学图像的检测和实时显示的荧光学***
JP2535712B2 (ja) X線ct装置
CN113658332A (zh) 一种基于超声影像的腹直肌智能分割重建方法及装置

Legal Events

Date Code Title Description
ST Notification of lapse

Effective date: 20101029