FR2978273A1 - Procede de reduction du bruit dans une sequence d'images fluoroscopiques par filtrage temporel et spatial - Google Patents
Procede de reduction du bruit dans une sequence d'images fluoroscopiques par filtrage temporel et spatial Download PDFInfo
- Publication number
- FR2978273A1 FR2978273A1 FR1102306A FR1102306A FR2978273A1 FR 2978273 A1 FR2978273 A1 FR 2978273A1 FR 1102306 A FR1102306 A FR 1102306A FR 1102306 A FR1102306 A FR 1102306A FR 2978273 A1 FR2978273 A1 FR 2978273A1
- Authority
- FR
- France
- Prior art keywords
- image
- acquired
- curvelets
- coefficients
- value
- 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
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000002123 temporal effect Effects 0.000 title claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 27
- 230000001131 transforming effect Effects 0.000 claims abstract description 10
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims 1
- 239000010931 gold Substances 0.000 claims 1
- 229910052737 gold Inorganic materials 0.000 claims 1
- 230000009467 reduction Effects 0.000 abstract description 16
- 238000003384 imaging method Methods 0.000 abstract description 5
- 238000002594 fluoroscopy Methods 0.000 abstract description 4
- 238000011282 treatment Methods 0.000 description 10
- 230000003044 adaptive effect Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 206010002329 Aneurysm Diseases 0.000 description 1
- 206010057469 Vascular stenosis Diseases 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000010109 chemoembolization Effects 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012978 minimally invasive surgical procedure Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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
- G06T2207/10121—Fluoroscopy
-
- 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/20004—Adaptive image processing
- G06T2207/20012—Locally adaptive
-
- 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/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- 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/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
L'invention se situe dans le domaine de l'imagerie médicale par fluoroscopie. Elle concerne un procédé de réduction du bruit dans une séquence d'images fluoroscopiques acquises par un détecteur à rayons X. Le procédé selon l'invention comprend les étapes successives suivantes pour chaque image (x ) de la séquence : ▪ appliquer un filtrage temporel (22) sur l'image (x ) acquise à un instant n, ledit filtrage temporel comportant un traitement de réduction du bruit fluoroscopique (FNR), ▪ appliquer un filtrage spatial (23) sur l'image (y ) acquise à l'instant n, ledit filtrage spatial comportant les sous-étapes suivantes : - transformer (231) l'image (y ) acquise à l'instant n du domaine spatial au domaine des curvelets par une transformation en curvelets, chaque image transformée étant représentée par un ensemble de coefficients, - seuiller (232) les coefficients de l'image (y ) par une fonction de seuillage, la fonction de seuillage annulant les coefficients inférieurs à un troisième seuil prédéterminé, et conservant ou ajustant les coefficients supérieurs au troisième seuil prédéterminé, - transformer (233) l'image (z ) dont les coefficients ont été seuillés du domaine des curvelets au domaine spatial par une transformation en curvelets inverse.
Description
PROCÉDÉ DE REDUCTION DU BRUIT DANS UNE SÉQUENCE D'IMAGES FLUOROSCOPIQUES PAR FILTRAGE TEMPOREL ET SPATIAL L'invention se situe dans le domaine de l'imagerie par rayons X et, plus spécifiquement, celui de l'imagerie médicale par fluoroscopie. Elle concerne un procédé de réduction du bruit dans une séquence d'images fluoroscopiques acquises par un détecteur à rayons X.
Dans le domaine de l'imagerie médicale, l'imagerie par fluoroscopie à rayons X permet de fournir un flux d'images d'un patient à un médecin lors d'interventions chirurgicales dites minimalement invasives, par exemple la chimio-embolisation du foie, la vertébroplastie, la cathérisation des anévrismes ou le traitement des rétrécissements vasculaires. Les images peuvent notamment aider le guidage d'instruments chirurgicaux. Une intervention guidée par fluoroscopie est typiquement réalisée en faisant passer un cathéter à l'intérieur du réseau vasculaire du patient. Un agent de contraste peut ou non avoir été injecté préalablement dans le réseau vasculaire afin de l'opacifier et d'améliorer la visibilité des vaisseaux. Une telle intervention est généralement relativement longue et la dose de rayons X auxquels le patient est soumis doit être limitée afin d'éviter d'occasionner des lésions ou des brûlures des tissus. Du fait de cette limitation de la dose de rayons X utilisée, les images fluoroscopiques obtenues comportent un niveau de bruit relativement élevé, et donc un rapport contraste sur bruit relativement faible, rendant les images difficilement lisibles.
Des traitements de filtrage sont mis en oeuvre afin de diminuer le bruit quantique présent dans ces images et augmenter leur rapport contraste 25 sur bruit. Un traitement connu sous le nom de traitement de réduction du bruit fluoroscopique, ou traitement FNR selon le sigle de l'expression anglo- saxonne "Fluoroscopic Noise Reduction", est ainsi généralement réalisé sur les images acquises par le détecteur à rayons X. L'objectif de ce traitement 30 FNR est de filtrer le bruit présent dans l'image tout en préservant le contraste de l'information présente dans cette image. Le traitement FNR est réalisé en appliquant un filtre temporel aux régions des images où il n'y a pas de déplacement. L'existence ou l'absence de déplacements dans des régions de l'image est détectée à partir de la variation individuelle de l'intensité de chacun des pixels considérés séparément. Un pixel est considéré comme étant en déplacement lorsque sa variation d'intensité entre deux images dépasse un seuil lié à l'écart type du bruit. Ces pixels déterminés comme étant en déplacement ne sont pas, ou très peu, filtrés. En revanche, un traitement par filtrage temporel récursif est appliqué aux pixels déterminés comme étant fixes. Le traitement FNR présente l'avantage de ne pas colorer le bruit spatialement. Autrement dit, il ne fait pas apparaître d'objets fictifs résultant de groupements de pixels d'intensités voisines entre elles. Cependant, le traitement FNR présente une capacité de débruitage relativement limitée et tend à faire disparaître les objets d'intérêts dans l'image. A côté du traitement par filtrage temporel FNR, les traitements par filtrage spatial présente des propriétés intéressantes. Parmi les différentes variantes de traitement par filtrage spatial, le filtrage par transformation de l'image dans le domaine des ondelettes est l'une des techniques les plus utilisées. La théorie des ondelettes, introduite dans les années 1980, a été développée à partir de la transformée de Fourier. La transformée de Fourier est utilisée pour décrire de nombreux phénomènes physiques. Elle constitue en effet un outil mathématique performant pour décrire de nombreux phénomènes physiques, notamment les phénomènes physiques modélisés par des signaux stationnaires. En revanche, les représentations des phénomènes physiques modélisés par des signaux non stationnaires ne sont généralement pas satisfaisantes. En particulier, la transformée de Fourier est incapable de localiser les portions du signal où la fréquence change subitement. Pour y remédier, une transformée de Fourier à fenêtres a été développée. Elle consiste à multiplier le signal par une fenêtre de dimension donnée. L'analyse devient donc locale, la dimension de la fenêtre déterminant la résolution temporelle obtenue, ou la résolution spatiale dans le cas de signaux à deux dimensions. La théorie des ondelettes s'appuie sur la transformée de Fourier à fenêtres. Elle s'en différencie de par le fait qu'elle permet, au cours de la transformation, de changer la dimension de la fenêtre. Pour cette raison, on parle d'analyse multirésolution ou analyse multi-échelle.
Une autre différence est que le signal n'est plus décomposé en une somme de sinus et de cosinus, mais selon des fonctions appelées "ondelettes". Les ondelettes présentent l'avantage de transformer l'image en une représentation dite creuse. Par représentation creuse, on entend une représentation nécessitant un faible nombre de paramètres pour représenter fidèlement l'image. Cependant, les ondelettes sont mal adaptées à la représentation de discontinuités de type ligne, telles que les contours d'objets, car leur forme ne leur permet pas d'exploiter la forme des objets. En conséquence, l'information sur ces discontinuités est contenue dans les coefficients de nombreuses échelles. Pour obtenir une faible erreur de reconstruction, il faut donc calculer un grand nombre de coefficients, entraînant un coût de calcul important. Pour une application de fluoroscopie où les images doivent être traitées "en temps réel", c'est-à-dire dans des durées inférieures à la période de rafraichissement des images, ce coût de calcul implique des ressources matérielles importantes.
D'autres théories ont été développées pour réduire la redondance des informations. Il s'agit notamment des ondelettes isotropes et anisotropes, ou encore des arbres duaux orthogonaux ou non. Au début des années 2000, la théorie des curvelets est apparue. Cette théorie est également fondée sur le principe de l'analyse mufti- échelles. La géométrie des curvelets leur permet de tirer profit de la régularité des contours des objets. Une image peut ainsi être représentée par un nombre relativement faible de coefficients. Cependant, la transformation en curvelets présente l'inconvénient d'introduire des artéfacts dans l'image, notamment une coloration spatiale du bruit. L'introduction d'artéfacts est d'autant plus forte que le rapport contraste sur bruit est faible. Or, dans le domaine de l'imagerie par fluoroscopie, le rapport contraste sur bruit est souvent inférieur à 1. La présence d'artéfacts apparaît donc comme étant inévitable et a limité jusque là son utilisation à des images dont le rapport contraste sur bruit est bien plus élevé.
Un but de l'invention est notamment de proposer un traitement de débruitage adapté aux images fluoroscopiques dont le rapport contraste sur bruit est relativement faible. En particulier, le traitement de débruitage doit présenter une forte capacité de débruitage, engendrer un faible coût de calcul, limiter la perte des objets d'intérêt dans l'image, et ne pas introduire d'artéfacts. A cet effet, l'invention a pour objet un procédé de réduction du bruit dans une séquence d'images fluoroscopiques acquises par un détecteur à rayons X, chaque image étant formée dans le domaine spatial d'une matrice de pixels, chaque pixel ayant une valeur représentative d'un niveau de signal, le procédé comprenant les étapes successives suivantes pour chaque image : ^ appliquer un filtrage temporel sur l'image acquise à un instant n, l'application dudit filtrage temporel comportant les sous-étapes suivantes pour chaque pixel de l'image : - déterminer une différence entre la valeur du pixel considéré dans l'image acquise à l'instant n et la valeur du pixel correspondant dans l'image acquise à un instant n-1 sur laquelle le filtrage temporel a été appliqué, - si la différence est inférieure à un premier seuil prédéterminé, 15 corriger la valeur du pixel correspondant avec un premier ensemble de paramètres, - si la différence est supérieure à un deuxième seuil prédéterminé, corriger la valeur du pixel correspondant avec un deuxième ensemble de paramètres, 20 ^ appliquer un filtrage spatial sur l'image acquise à l'instant n, l'application dudit filtrage spatial comportant les sous-étapes suivantes : - transformer l'image acquise à l'instant n du domaine spatial au domaine des curvelets par une transformation en curvelets, chaque image dans le domaine des curvelets étant représentée par un ensemble de 25 coefficients, - seuiller les coefficients de l'image par une fonction de seuillage, la fonction de seuillage annulant les coefficients inférieurs à un troisième seuil prédéterminé, et conservant ou ajustant les coefficients supérieurs au troisième seuil prédéterminé, 30 - transformer l'image dont les coefficients ont été seuillés du domaine des curvelets au domaine spatial par une transformation en curvelets inverse.
L'étape de transformation de l'image acquise à l'instant n dans le domaine des curvelets utilise avantageusement une transformée en curvelets à 6 ou 9 échelles, et/ou 16 orientations.
Selon une première forme particulière de réalisation, l'étape de transformation de l'image acquise à l'instant n dans le domaine des curvelets utilise une transformée en curvelets discrète par USFFT. Selon une deuxième forme particulière de réalisation, l'étape de 10 transformation de l'image acquise à l'instant n dans le domaine des curvelets utilise une transformée en curvelets discrète par wrapping. La fonction de seuillage du filtrage spatial peut être une fonction de seuillage dur, c'est-à-dire une fonction de seuillage selon laquelle les 15 coefficients sont annulés s'ils sont inférieurs au troisième seuil prédéterminé, et conservés sinon. Le troisième seuil prédéterminé TS du filtrage spatial peut être déterminé par les étapes successives suivantes : 20 ^ créer une image Y de mêmes dimensions que l'image acquise à l'instant n et de valeur uniforme 1 sur tous les pixels, ^ transformer l'image Y par une transformée de Fourier, ^ normaliser l'amplitude F du pic dans la transformée de Fourier de l'image Y par la relation : 6 F IFI _ 5 M x P 25 où M et P désignent les dimensions de l'image Y, ^ transformer la transformée de Fourier de l'image Y du domaine spatial au domaine des curvelets par une transformation en curvelets, chaque image dans le domaine des curvelets étant représentée par un ensemble de coefficients ci,k, 30 ^ déterminer une norme c i des coefficients c l,k par la relation : Y z y _ El cl,l,k M x P où M et P désignent les dimensions de l'image Y et l désigne les orientations des curvelets, ^ déterminer le seuil prédéterminé TS par la relation : TS=r.c j.Q où o- est l'écart type du bruit dans l'image acquise à l'instant n et sur laquelle le filtrage temporel a été appliqué, et où le facteur r prend la valeur 0 pour l'échelle la plus grossière des curvelets, la valeur 4 pour l'échelle la plus fine, et la valeur 1 pour les autres échelles.
Selon une forme particulière de réalisation, l'étape d'application du filtrage temporel sur l'image acquise à l'instant n comprend, à la suite de la sous-étape de détermination de la différence entre la valeur du pixel considéré dans l'image acquise à l'instant n et la valeur du pixel correspondant dans l'image acquise à l'instant n-1 sur laquelle le filtrage temporel a été appliqué, la sous-étape supplémentaire suivante : ^ si la différence est comprise entre le premier seuil prédéterminé et le deuxième seuil prédéterminé, corriger la valeur du pixel correspondant avec un troisième ensemble de paramètres.
En résumé, le procédé de réduction de bruit selon l'invention comporte successivement un filtrage temporel et un filtrage spatial par transformation en curvelets. L'ordre dans lequel sont effectués les filtrages présente une importance en ce qui concerne la coloration du bruit. En effet, le filtrage spatial par transformation en curvelets introduit des artéfacts s'il n'est pas précédé du filtrage temporel. Par ailleurs, l'invention a notamment pour avantage de ne pas nécessiter d'a priori de formes pour le filtrage spatial.
L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description qui va suivre, faite en regard de 30 dessins annexés sur lesquels : - la figure 1 représente un exemple de points pris en compte pour déterminer le contraste d'un objet dans une image fluoroscopique ; - la figure 2 représente des étapes possibles du procédé de réduction de bruit selon l'invention ; - les figures 3A et 3B représentent par des courbes une fonction de seuillage dur et une fonction de seuillage doux, respectivement ;
- la figure 4 illustre, par un graphique, la performance de ce procédé de réduction de bruit pour une séquence d'images fluoroscopiques représentant un objet ayant un mouvement de translation et de rotation. Pour la suite de la description, on considère une séquence d'images fluoroscopiques xn acquises par un détecteur à rayons X, où n désigne les différents instants d'acquisition dans la séquence. Chaque image est composée d'une matrice de pixels ayant chacun une valeur xn(m,p) représentative d'un niveau de signal, avec m et p désignant respectivement la ligne et la colonne de la matrice. La valeur est par exemple une valeur numérique comprise entre 0 et 255. La valeur numérique est ainsi codée sur 8 bits. Elle peut également être codée sur 14 bits. Le rapport contraste sur bruit (CNR) est défini comme suit : 'contraste - bruit' CNR = a où le bruit est de préférence déterminé sur une zone de l'image où il n'y a pas d'objets d'intérêt, où le contraste est de préférence déterminé à partir de points correspondant à un objet d'intérêt dans l'image, et où Q désigne l'écart type du bruit, c'est-à-dire l'écart type des valeurs de pixels correspondant à une zone de l'image ne comportant pas d'objets d'intérêt. La figure 1 représente les points considérés (en blanc) pour déterminer le contraste d'une endoprothèse vasculaire, mieux connue sous la dénomination anglo-saxonne de "stent".
La figure 2 représente des étapes possibles pour le procédé de réduction de bruit dans la séquence d'images fluoroscopiques. Le procédé est décrit pour l'image xn acquise à l'instant n. Il peut cependant être appliqué successivement sur chacune des images de la séquence.
Dans une première étape 21, l'image xn est traitée par une chaîne de prétraitement. Le prétraitement peut notamment comporter des regroupements de pixels adjacents dans l'image. Les pixels peuvent être regroupés par exemple en fonction de la résolution souhaitée ou du niveau de bruit souhaité. Le prétraitement peut également comporter une phase d'étalonnage du détecteur à rayons X. En particulier, le gain et l'offset à appliquer à chaque pixel du détecteur peuvent être déterminés. En outre, l'étalonnage peut également comprendre une identification des pixels prenant une valeur aléatoire ou une valeur constante quel que soit le rayonnement reçu. Ces pixels défectueux peuvent être corrigés en interpolant leur valeur par les valeurs des pixels voisins. Le prétraitement peut aussi comporter une remise à l'échelle des images, c'est-à-dire un ajustement global de la valeur de tous les pixels de chaque image, de manière à ce que la valeur moyenne des pixels soit constante pour toutes les images de la séquence. Enfin, le prétraitement peut comporter une sous-étape d'adaptation des images. Étant donnée la physique des rayons X, le bruit présent dans chaque image xn est un bruit poissonien. Or un tel bruit est inadapté aux traitements de réduction de bruit. De plus, chaque image xn est une image représentative des atténuations des rayons X, alors qu'une image des épaisseurs radiologiques est plus naturelle à analyser. La formule 15 d'Ascombe permet d'obtenir une image des épaisseurs radiologiques avec un bruit pseudo-gaussien : xn (m, p) _ où xo,n(m, p) est la valeur du pixel (m, p) avant la sous-étape d'adaptation de l'image xn. Dans une deuxième étape 22, l'image xn est filtrée par un filtre 20 temporel. Le filtrage temporel est par exemple un filtrage appelé "Fluoroscopic Noise Reduction" (FNR). Ce filtrage est décrit notamment dans le brevet US 6,314,160 B1. Soit yn l'image xn obtenue après le filtrage temporel. On note 0(m, p), la différence entre la valeur xn(m, p) dans l'image dans l'image filtrée xn acquise à l'instant n et la valeur yn_1(m, p) 25 yn_1 correspondant à l'acquisition à l'instant n-1 : 0(m, p) = xn (m, p) - yn-1(m, p) Pour alléger les notations, on notera simplement xn la valeur d'un pixel de l'image xn, yn la valeur d'un pixel de l'image yn et 0 la différence entre les deux valeurs xn et yn_1 du pixel considéré. Le filtrage FNR se décline selon deux modes, appelés mode adaptatif et mode non adaptatif. Le mode 30 adaptatif est appelé ainsi du fait de sa prise en compte du mouvement des objets. Le filtrage FNR consiste à comparer pixel à pixel, la différence 0 à des seuils prédéterminés. Le premier seuil SI correspond à une différence 0 relativement faible. Le deuxième seuil S2 correspond à une différence A relativement forte. Le filtrage FNR est alors paramétré par des coefficients a et [3 en fonction du résultat de cette différence A et du mode considéré. On note a la force du filtre : ~2 a= y a + Qx où a-x est l'écart type dans l'image xn et où Qy est l'écart type dans l'image Yn-1 On définit aussi les constantes AA, BB, CC et DD, comme suit : ^ AA = 1=a 2.ax ^ BB=2a-1 Q(1-a) ^ CC = 2 ^ DD = i Q(1- a)
Pour le mode adaptatif, les seuils SI et S2 sont définis comme suit : ^ Ss=a" ^ S2=2Q et les coefficients a et R sont déterminés comme suit : ^ si IAI < S1, lIa = a 1(3=0 AA.AZ+BB.I oI +CC ^ si SI < 1AI < S2, a = lol /3=0
a=1-II ^ sinon, Pour le mode non adaptatif, les seuils SI et S2 sont définis comme suit : 20 ^ Ss = 26 ^ S2 Ja et les coefficients a et [3 sont déterminés comme suit : ^ si 14I < S1, fa= a = 0 a _ a.°Z ^ si SI < IA1 < S2, 32.a.)2 4 fi' = 3a.a a=1 25 ^sinon, _ a 4.v ~ 3 va Le filtrage FNR proprement dit consiste à déterminer une nouvelle valeur yn pour chaque pixel de l'image à partir de sa valeur xn et des coefficients a et R: yn=a.xn+(1-a)yn-1+fi' Autrement dit, la valeur xn de chaque pixel est corrigée à partir des coefficients a et R. Sur la figure 2, l'étape 22 de filtrage FNR est représentée pour le mode adaptatif. A titre d'exemple, dans une première sous-étape 221, la valeur absolue de la différence 0 est comparée pour chaque pixel xn au seuil S1, à savoir l'écart type a. Si la valeur absolue de la différence à est inférieure au seuil S1, le pixel considéré représente probablement un niveau de bruit et un filtrage dit fort est appliqué dans une sous-étape 222 avec les paramètres correspondants. Sinon, dans une sous-étape 223, la valeur absolue de la différence 0 est comparée pour chaque pixel xn au seuil S2, à savoir 2a. Si la valeur absolue de la différence 0 est supérieure au seuil S2, le pixel considéré correspondant probablement à un point de l'objet ayant bougé et un filtrage dit faible est appliqué dans une sous-étape 224 avec les paramètres correspondants. Enfin, si la valeur absolue de la différence 0 est comprise entre les seuils S1 et S2, une incertitude existe. Afin d'obtenir une fonction continue entre le filtrage fort et le filtrage faible, un filtrage moyen est appliqué dans une sous-étape 225 avec les paramètres correspondants. Dans une troisième étape 23, l'image yn obtenue après filtrage temporel est filtrée par un filtre spatial basé sur la transformée en curvelets. La théorie des curvelets utilise le principe des espaces multirésolution ou multi-échelle, une échelle correspondant à un niveau de partition de l'espace. L'espace est partitionné en curvelets dont les enveloppes sont de dimensions 2j x 2j/2, où j est un entier positif désignant l'échelle. Chaque curvelet dépend du nombre d'échelles considérées, de sa position dans l'image et de son orientation. Pour chaque échelle, la meilleure approximation possible de l'image est calculée. D'une échelle à l'autre, seules les améliorations obtenues sont représentées. Le nombre d'échelles détermine la qualité de la reconstruction de l'image. Avec les curvelets, un nombre d'échelles compris entre 6 et 10 permet de fournir des résultats très satisfaisants en termes de rapport contraste sur bruit.
Dans une première sous-étape 231 de l'étape de filtrage spatial 23, l'image yn est transformée du domaine spatial au domaine des curvelets par une transformation en curvelets. De préférence, une transformation en curvelets discrète est utilisée. Une telle transformation est notamment décrite dans E. Candès, Demanet L., Donoho D. and Ying L., "Fast Discrete Curvelet Transforms", Multiscale Model. Simul., vol. 5, no. 3, 2006, pp. 861-899. Sous leur forme discrète, les curvelets ne peuvent pas être échantillonnées selon les grilles de l'espace de Fourier. Le document de Donoho et al propose ainsi deux solutions pour translater les curvelets à échelle et orientation données sur la grille spatiale. Une première solution pour transformer l'image dans le domaine des curvelets discrètes est appelée transformée en curvelets discrète via USFFT, où USFFT est le sigle anglo-saxon pour "unequally-spaced fast Fourier transform" (USFFT), que l'on peut traduire en français par "transformée de Fourier rapide pour données à répartition irrégulière". Dans cette solution, la grille de translation est inclinée de manière à être alignée dans la direction des curvelets à échelle et angle donnés. Une deuxième solution est appelée transformée en curvelets discrète via wrapping, où le terme "wrapping" peut être traduit en français par "enroulement". Dans cette solution, la grille de translation est identique pour tous les angles d'un même quadrant à une échelle donnée. Cela implique que le nombre d'angles utilisés pour la transformée soit un multiple de 4. La transformée en curvelets via USFFT est plus fidèle à la théorie continue des curvelets que la transformée en curvelets via wrapping. Cependant, son coût de calcul est plus important. Dans les deux solutions, la transformée en curvelets de l'image yn peut être représentée par un ensemble de coefficients c°(j, 1, k), où l'exposant D fait référence à la transformée discrète, et où j, l et k désignent respectivement l'échelle, l'orientation et la position de la curvelet considérée. Selon une forme préférentielle de réalisation du filtre spatial, le nombre d'échelles utilisées est égal à 9. Selon une autre forme de réalisation, le nombre d'échelles utilisées est égal à 6. Cette dernière forme de réalisation implique un nombre plus faible de coefficients à calculer et donc une durée de traitement inférieure. Par ailleurs, le nombre d'orientations utilisées est avantageusement égal à 16. Ce nombre présente un bon compromis entre la résolution en orientation et le coût de calcul des curvelets, notamment pour des images réelles dont les courbes nécessitent un grand nombre d'orientations pour les approcher. De plus, le nombre 16 étant un multiple de 4, il facilite la transformée en curvelets par wrapping. Dans une deuxième sous-étape 232 de l'étape de filtrage spatial 35 23, les coefficients c°(j, 1, k) de l'image yn sont seuillés par une fonction de seuillage. L'étape de seuillage 232 est basée sur l'hypothèse selon laquelle certains coefficients cD(j, 1, k) représentent le bruit alors que d'autres contiennent l'information utile. L'étape de seuillage 232 consiste donc à identifier et à supprimer les coefficients ne contenant pas d'informations utiles. La fonction de seuillage se caractérise d'une part par sa forme et d'autre part par son seuil T. Il existe principalement deux fonctions de seuillage, la fonction de seuillage dit doux et la fonction de seuillage dit dur. En simplifiant la notation des coefficients par en et en appelant en les coefficients seuillés, la fonction de seuillage dur peut être modélisée comme suit : 1cn sinon e _ n sinon La fonction de seuillage doux peut être modélisée comme suit : _ si Icni T en - IO cn - T. signe(cn) sinon La figure 3A représente par une courbe 31 la fonction de seuillage dur et la figure 3B représente par une courbe 32 la fonction de seuillage doux. La fonction de seuillage doux a une capacité de débruitage plus importante que la fonction de seuillage dur. En revanche, elle modifie la valeur des pixels en fonction du seuil T choisi. Il en résulte une perte de précision sur certains détails de l'image. La détermination du seuil T est une étape cruciale car elle détermine la performance du filtrage spatial, tant en ce qui concerne la capacité de débruitage que l'introduction d'artéfacts.
Selon une première forme de réalisation, le seuil T est déterminé comme étant le seuil universel Tu décrit notamment dans D.L. Donoho and Johnstone, I.M., "Ideal spatial adaptation via wavelet shrinkage", Biometrica, vol. 81, 1994, pp. 425-455. Ce seuil Tu est défini par la relation suivante : logN Tu _ 2 N où N est le nombre d'échantillons, c'est-à-dire ici le nombre de pixels dans l'image yn. Le seuil universel TU n'est pas optimal et peut introduire des artéfacts. Néanmoins, des seuils proportionnels au seuil universel TU peuvent présenter de bons résultats si le facteur de proportionnalité est bien choisi. Selon une deuxième forme de réalisation, le seuil T est déterminé par la méthode décrite dans J.L. Starck, Candès I.J. and Donoho L., "The Curvelet Transform for Image Denoising", IEEE Transactions on Image Processing, vol. 11, no. 6, 2002, pp. 670-684. Cette méthode consiste à créer une image, notée Y, de mêmes dimensions que l'image yn à débruiter et de valeur uniforme 1 sur tous les pixels. L'image Y est transformée par la transformée de Fourier. La fréquence de coordonnées (0,0) dans l'image Y est centrée et son amplitude est normalisée par les dimensions de l'image. Soit F l'amplitude de la fréquence centrale. Son amplitude normalisée, notée IFI, est par exemple réalisée selon la relation suivante : F IFI = 5 VM x p où M et P désignent les dimensions de l'image Y, c'est-à-dire également de l'image xn et de l'image yn.
On obtient ainsi une image avec un pic de Dirac centré sur l'image, dont l'amplitude est fonction des dimensions de l'image. Le seuil de cette méthode, noté TS dépend de l'échelle j et de l'orientation 1. Il est déterminé en fonction de la norme des coefficients ci,k des curvelets de l'image Y pour une échelle j et une orientation l données. Cette norme c est déterminée par la relation suivante : 2 y El cjY,I,k - M x P
La fonction de seuillage est une fonction de seuillage dur, définie comme suit : Cn = 1cn sinon Le facteur multiplicatif r est déterminé en fonction de l'échelle j. Pour l'échelle la plus fine (par exemple j = 9), le facteur r est égal à 4. Pour les autres échelles, le facteur r est égal à 3. Ces facteurs présentent une bonne capacité de débruitage pour des images présentant un fort rapport contraste sur bruit. Cependant, pour des images fluoroscopiques, des artéfacts sont introduits. Ces artéfacts peuvent être évités en adaptant le facteur r. En particulier, pour l'échelle la plus grossière (j = 1), le facteur r peut être égal à 0. Il n'y a donc pas de seuillage. Pour l'échelle la plus fine, le facteur r peut encore être égal à 4. Et pour les autres échelles, le facteur r peut être égal à 1. Dans une troisième sous-étape 233 de l'étape de filtrage spatial 23, l'image yn est transformée du domaine des curvelets au domaine spatial (matriciel) par une transformation en curvelets inverse. A cet effet, les 0si Icnl <r.ci .6 coefficients seuillés sont utilisés. L'image obtenue après filtrage temporel et filtrage spatial est notée zn.
La figure 4 illustre par un graphique la performance du procédé de réduction de bruit selon l'invention pour une séquence d'images représentant un stent ayant un mouvement de translation et de rotation. Sur le graphique, l'axe des abscisses représente le CNR moyen de la séquence d'images xn et l'axe des ordonnées représente le CNR moyen de la séquence d'images obtenue après un traitement de réduction de bruit. Une première courbe 41 représente le CNR moyen de la séquence d'images Zn obtenue par le procédé de réduction de bruit selon l'invention, c'est-à-dire par un filtrage temporel suivi d'un filtrage spatial par transformation en curvelets. Une deuxième courbe 42 représente le CNR moyen de la même séquence d'images après traitement par le même filtrage temporel suivi d'un filtrage spatial par transformation en ondelettes. Le graphique montre que pour une séquence d'images à traiter dont le CNR est relativement faible (0,5 - 0,7 - 1), le procédé de réduction de bruit selon l'invention présente une amélioration notable en termes de CNR. La réduction de bruit est encore bien meilleure pour une séquence d'images dont le CNR est plus élevé (1,5 - 2 - 2,5). De plus, le procédé de réduction de bruit selon l'invention permet de réduire significativement l'introduction d'artéfacts.
Claims (9)
- REVENDICATIONS1. Procédé de réduction du bruit dans une séquence d'images fluoroscopiques acquises par un détecteur à rayons X, chaque image (xn) étant formée dans le domaine spatial d'une matrice de pixels, chaque pixel ayant une valeur représentative d'un niveau de signal, le procédé comprenant les étapes successives suivantes pour chaque image (xn) : ^ appliquer un filtrage temporel (22) sur l'image (xn) acquise à un instant n, l'application dudit filtrage temporel comportant les sous-étapes suivantes pour chaque pixel de l'image : - déterminer (221, 223) une différence entre la valeur du pixel considéré dans l'image (xn) acquise à l'instant n et la valeur du pixel correspondant dans l'image (yn_1) acquise à un instant n-1 sur laquelle le filtrage temporel (22) a été appliqué, - si la différence est inférieure à un premier seuil prédéterminé, corriger (222) la valeur du pixel correspondant avec un premier ensemble de 15 paramètres, - si la différence est supérieure à un deuxième seuil prédéterminé, corriger (224) la valeur du pixel correspondant avec un deuxième ensemble de paramètres, ^ appliquer un filtrage spatial (23) sur l'image (yn) acquise à l'instant n, 20 l'application dudit filtrage spatial comportant les sous-étapes suivantes : - transformer (231) l'image (yn) acquise à l'instant n du domaine spatial au domaine des curvelets par une transformation en curvelets, chaque image dans le domaine des curvelets étant représentée par un ensemble de coefficients, 25 - seuiller (232) les coefficients de l'image (yn) par une fonction de seuillage, la fonction de seuillage annulant les coefficients inférieurs à un troisième seuil prédéterminé, et conservant ou ajustant les coefficients supérieurs au troisième seuil prédéterminé, - transformer (233) l'image (yn) dont les coefficients ont été 30 seuillés du domaine des curvelets au domaine spatial par une transformation en curvelets inverse.
- 2. Procédé selon la revendication 1, dans lequel l'étape (231) de transformation de l'image (yn) acquise à l'instant n dans le domaine des curvelets utilise une transformée en curvelets à 6 échelles.
- 3. Procédé selon la revendication 1, dans lequel l'étape (231) de transformation de l'image (yn) acquise à l'instant n dans le domaine des curvelets utilise une transformée en curvelets à 9 échelles.
- 4. Procédé selon l'une des revendications précédentes, dans 10 lequel l'étape (231) de transformation de l'image (yn) acquise à l'instant n dans le domaine des curvelets utilise une transformée en curvelets comprenant 16 orientations.
- 5. Procédé selon l'une des revendications précédentes, dans 15 lequel l'étape (231) de transformation de l'image (A) acquise à l'instant n dans le domaine des curvelets utilise une transformée en curvelets discrète par USFFT.
- 6. Procédé selon l'une des revendications 1 à 4, dans lequel 20 l'étape (231) de transformation de l'image (yn) acquise à l'instant n dans le domaine des curvelets utilise une transformée en curvelets discrète par wrapping.
- 7. Procédé selon l'une des revendications précédentes, dans 25 lequel la fonction de seuillage est une fonction de seuillage dur, c'est-à-dire que les coefficients sont annulés s'ils sont inférieurs au troisième seuil prédéterminé, et conservés sinon.
- 8. Procédé selon la revendication 7, dans lequel le troisième 30 seuil prédéterminé TS est déterminé par les étapes successives suivantes : ^ créer une image Y de mêmes dimensions que l'image (yn) acquise à l'instant n et de valeur uniforme 1 sur tous les pixels, ^ transformer l'image Y par une transformée de Fourier, ^ normaliser l'amplitude F du pic dans la transformée de Fourier de 35 l'image Y par la relation :6 F 1FI _ 5 M x P où M et P désignent les dimensions de l'image Y, ^ transformer la transformée de Fourier de l'image Y du domaine spatial au domaine des curvelets par une transformation en curvelets, chaque image dans le domaine des curvelets étant représentée par un ensemble de coefficients Cjl,k, ^ déterminer une norme ê I des coefficients c jl,k par la relation : Y y zI Cl I,k a sil- M x P où M et P désignent les dimensions de l'image Y et l désigne les orientations des curvelets, ^ déterminer le seuil prédéterminé TS par la relation : TS = r. ci. Q où or est l'écart type du bruit dans l'image (A) acquise à l'instant n et sur laquelle le filtrage temporel (22) a été appliqué, et où le facteur r prend la valeur 0 pour l'échelle la plus grossière des curvelets, la valeur 4 pour l'échelle la plus fine, et la valeur 1 pour les autres échelles.
- 9. Procédé selon l'une des revendications précédentes, dans lequel l'étape (22) d'application du filtrage temporel sur l'image acquise (xn) à l'instant n comprend, à la suite de la sous-étape (221, 223) de détermination de la différence entre la valeur du pixel considéré dans l'image acquise (xn) à l'instant n et la valeur du pixel correspondant dans l'image (yn_1) acquise à l'instant n-1 sur laquelle le filtrage temporel (22) a été appliqué, la sous-étape supplémentaire suivante : ^ si la différence est comprise entre le premier seuil prédéterminé et le deuxième seuil prédéterminé, corriger (225) la valeur du pixel correspondant avec un troisième ensemble de paramètres.25
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1102306A FR2978273B1 (fr) | 2011-07-22 | 2011-07-22 | Procede de reduction du bruit dans une sequence d'images fluoroscopiques par filtrage temporel et spatial |
PCT/EP2012/064347 WO2013014110A1 (fr) | 2011-07-22 | 2012-07-20 | Procédé de réduction du bruit dans une séquence d'images fluoroscopiques par filtrage temporel et spatial |
EP12737811.5A EP2734975B1 (fr) | 2011-07-22 | 2012-07-20 | Procédé de réduction du bruit dans une séquence d'images fluoroscopiques par filtrage temporel et spatial |
US14/234,292 US9058654B2 (en) | 2011-07-22 | 2012-07-20 | Method for reducing noise in a sequence of fluoroscopic images by temporal and spatial filtering |
JP2014520683A JP6120843B2 (ja) | 2011-07-22 | 2012-07-20 | 時間フィルタリングおよび空間フィルタリングによる一連の蛍光透視画像における雑音を低減させるための方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1102306A FR2978273B1 (fr) | 2011-07-22 | 2011-07-22 | Procede de reduction du bruit dans une sequence d'images fluoroscopiques par filtrage temporel et spatial |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2978273A1 true FR2978273A1 (fr) | 2013-01-25 |
FR2978273B1 FR2978273B1 (fr) | 2013-08-09 |
Family
ID=46548470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1102306A Expired - Fee Related FR2978273B1 (fr) | 2011-07-22 | 2011-07-22 | Procede de reduction du bruit dans une sequence d'images fluoroscopiques par filtrage temporel et spatial |
Country Status (5)
Country | Link |
---|---|
US (1) | US9058654B2 (fr) |
EP (1) | EP2734975B1 (fr) |
JP (1) | JP6120843B2 (fr) |
FR (1) | FR2978273B1 (fr) |
WO (1) | WO2013014110A1 (fr) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9672597B2 (en) | 2013-06-28 | 2017-06-06 | Thales | Method for reducing noise in sequences of fluoroscopic images |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2654015A1 (fr) * | 2012-04-21 | 2013-10-23 | General Electric Company | Procédé, système et support lisible sur ordinateur pour le traitement dýune image vidéo médicale |
JP6381198B2 (ja) * | 2013-11-08 | 2018-08-29 | キヤノン株式会社 | 制御装置、制御方法及びプログラム |
FR3018147B1 (fr) * | 2014-03-03 | 2016-03-04 | Sagem Defense Securite | Debruitage video optimise pour systeme multicapteur heterogene |
CN104463325A (zh) * | 2014-12-04 | 2015-03-25 | 中国科学院电子学研究所 | 一种极地探冰雷达原始数据噪声抑制方法 |
CN105259539B (zh) * | 2015-10-29 | 2017-09-29 | 中国科学院大学 | 一种针对极地探冰雷达数据的分部式噪声抑制方法 |
US9934557B2 (en) | 2016-03-22 | 2018-04-03 | Samsung Electronics Co., Ltd | Method and apparatus of image representation and processing for dynamic vision sensor |
US11113791B2 (en) * | 2017-01-03 | 2021-09-07 | Flir Systems, Inc. | Image noise reduction using spectral transforms |
US10983236B2 (en) | 2017-06-20 | 2021-04-20 | Saudi Arabian Oil Company | Super-resolution radon transform based on thresholding |
JP7183746B2 (ja) * | 2018-12-03 | 2022-12-06 | コニカミノルタ株式会社 | 制御装置及びプログラム |
JP7257829B2 (ja) * | 2019-03-15 | 2023-04-14 | ソニー・オリンパスメディカルソリューションズ株式会社 | 画像処理装置、画像処理方法およびプログラム |
CN112051612A (zh) * | 2019-06-05 | 2020-12-08 | 河北地质大学 | 一种消除地震数据随机噪声的方法和系统 |
DE102019122667A1 (de) * | 2019-08-22 | 2021-02-25 | Schölly Fiberoptic GmbH | Verfahren zur Unterdrückung von Bildrauschen in einem Videobildstrom, sowie zugehöriges medizinisches Bildaufnahmesystem und Computerprogrammprodukt |
CN111766306B (zh) * | 2020-05-21 | 2023-06-27 | 广州多浦乐电子科技股份有限公司 | 一种全聚焦成像方法及装置 |
CN112633202B (zh) * | 2020-12-29 | 2022-09-16 | 河南大学 | 一种基于双重去噪联合多尺度超像素降维的高光谱图像分类算法 |
CN112800727B (zh) * | 2021-04-14 | 2021-07-20 | 北京三维天地科技股份有限公司 | 给pdf文件加批注的方法及应用系统 |
CN116883292B (zh) * | 2023-09-07 | 2023-11-28 | 上海海栎创科技股份有限公司 | 基于伪高斯的图像方向场获取方法、系统及计算机设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7542622B1 (en) * | 2003-06-02 | 2009-06-02 | The Trustees Of Columbia University In The City Of New York | Spatio-temporal treatment of noisy images using brushlets |
WO2009100032A1 (fr) * | 2008-02-05 | 2009-08-13 | Ntt Docomo, Inc. | Diminution du bruit et/ou du scintillement dans des séquences vidéo au moyen du traitement spatial et temporel |
EP2151193A1 (fr) * | 2007-05-09 | 2010-02-10 | Hitachi Medical Corporation | Dispositif à ultrasons et procédé d'amélioration de la qualité d'image du dispositif à ultrasons |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06169906A (ja) * | 1992-12-10 | 1994-06-21 | Toshiba Corp | X線画像診断装置 |
JP3326885B2 (ja) * | 1993-06-30 | 2002-09-24 | 株式会社島津製作所 | X線テレビジョン装置 |
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 |
US6335990B1 (en) * | 1997-07-03 | 2002-01-01 | Cisco Technology, Inc. | System and method for spatial temporal-filtering for improving compressed digital video |
US6281942B1 (en) * | 1997-08-11 | 2001-08-28 | Microsoft Corporation | Spatial and temporal filtering mechanism for digital motion video signals |
US6018600A (en) * | 1998-04-15 | 2000-01-25 | Arch Development Corp. | Sampling and reconstruction of signals and images including MR images of multiple regions |
US6314160B1 (en) | 1999-12-17 | 2001-11-06 | General Electric Company | Method and apparatus for performing fluoroscopic noise reduction |
US6907143B2 (en) * | 2001-05-16 | 2005-06-14 | Tektronix, Inc. | Adaptive spatio-temporal filter for human vision system models |
US7840625B2 (en) * | 2005-04-07 | 2010-11-23 | California Institute Of Technology | Methods for performing fast discrete curvelet transforms of data |
CN101536033A (zh) * | 2006-11-08 | 2009-09-16 | 皇家飞利浦电子股份有限公司 | 图像信号的降噪 |
US8553783B2 (en) * | 2008-12-29 | 2013-10-08 | Stmicroelectronics Asia Pacific Pte Ltd. | Apparatus and method of motion detection for temporal mosquito noise reduction in video sequences |
US8942454B2 (en) * | 2009-12-03 | 2015-01-27 | The United States Of America, As Represented By The Secretary, Department Of Health And Human Services | Signal to-noise enhancement in imaging applications using a time-series of images |
-
2011
- 2011-07-22 FR FR1102306A patent/FR2978273B1/fr not_active Expired - Fee Related
-
2012
- 2012-07-20 JP JP2014520683A patent/JP6120843B2/ja active Active
- 2012-07-20 WO PCT/EP2012/064347 patent/WO2013014110A1/fr active Application Filing
- 2012-07-20 EP EP12737811.5A patent/EP2734975B1/fr active Active
- 2012-07-20 US US14/234,292 patent/US9058654B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7542622B1 (en) * | 2003-06-02 | 2009-06-02 | The Trustees Of Columbia University In The City Of New York | Spatio-temporal treatment of noisy images using brushlets |
EP2151193A1 (fr) * | 2007-05-09 | 2010-02-10 | Hitachi Medical Corporation | Dispositif à ultrasons et procédé d'amélioration de la qualité d'image du dispositif à ultrasons |
WO2009100032A1 (fr) * | 2008-02-05 | 2009-08-13 | Ntt Docomo, Inc. | Diminution du bruit et/ou du scintillement dans des séquences vidéo au moyen du traitement spatial et temporel |
Non-Patent Citations (1)
Title |
---|
M. TOMIC ET AL: "Adaptive spatio-temporal denoising of fluoroscopic X-ray sequences", BIOMEDICAL SIGNAL PROCESSING AND CONTROL, vol. 7, no. 2, 5 March 2011 (2011-03-05), pages 173 - 179, XP055021074, ISSN: 1746-8094, DOI: 10.1016/j.bspc.2011.02.003 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9672597B2 (en) | 2013-06-28 | 2017-06-06 | Thales | Method for reducing noise in sequences of fluoroscopic images |
Also Published As
Publication number | Publication date |
---|---|
EP2734975A1 (fr) | 2014-05-28 |
US9058654B2 (en) | 2015-06-16 |
WO2013014110A1 (fr) | 2013-01-31 |
JP6120843B2 (ja) | 2017-04-26 |
FR2978273B1 (fr) | 2013-08-09 |
US20140153842A1 (en) | 2014-06-05 |
JP2014520641A (ja) | 2014-08-25 |
EP2734975B1 (fr) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2978273A1 (fr) | Procede de reduction du bruit dans une sequence d'images fluoroscopiques par filtrage temporel et spatial | |
CN111081354B (zh) | 用于通过深度学习网络对医疗图像进行去噪的系统和方法 | |
US7362845B2 (en) | Method and apparatus of global de-noising for cone beam and fan beam CT imaging | |
US7457451B2 (en) | Systems, methods and apparatus for specialized filtered back-projection reconstruction for digital tomosynthesis | |
FR2752975A1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimentionnelle angiographique | |
FR2779853A1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique | |
FR2924254A1 (fr) | Procede de traitement d'images en radioscopie interventionnelle | |
CN111110260B (zh) | 一种图像重建方法、装置及终端设备 | |
FR2777374A1 (fr) | Procede de recalage de deux images differentes d'un meme objet | |
FR2848007A1 (fr) | Procede et systeme pour accentuer une image de tomosynthese en utilisant un filtrage transversal | |
Maini et al. | Peformance evaluation of various speckle noise reduction filters on medical images | |
Khare et al. | Daubechies complex wavelet transform based multilevel shrinkage for deblurring of medical images in presence of noise | |
Prakash et al. | Medical image denoising based on soft thresholding using biorthogonal multiscale wavelet transform | |
CN115795272A (zh) | 一种基于分数阶迭代离散小波变换的谱线去噪方法 | |
Yang et al. | Image deblurring using empirical Wiener filter in the curvelet domain and joint non-local means filter in the spatial domain | |
Benammar et al. | Sinogram interpolation method for limited-angle tomography | |
FR2835949A1 (fr) | Procede de tomosynthese a reconstruction multiplan | |
Lew et al. | Adaptive Gaussian Wiener Filter for CT-Scan Images with Gaussian Noise Variance | |
Bahendwar et al. | A comparative performance analysis of discrete wavelet transforms for denoising of medical images | |
Saoji et al. | Speckle and rician noise removal from medical images and Ultrasound images | |
Ma | Deconvolution using singular integral regularization and curvelet shrinkage | |
Shakeel et al. | Efficient blind image deblurring method for palm print images | |
Kadri et al. | Power shrinkage—curvelet domain image denoising using a new scale-dependent shrinkage function | |
Sugantharathnam et al. | Removal of poisson noise in medical images using fdct integrated with rudin–osher–fatemi model | |
Chen et al. | Research on laser one-dimensional range profile denoising and deconvolution |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 5 |
|
ST | Notification of lapse |
Effective date: 20170331 |