FR2596221A1 - Systeme de filtre median adaptif - Google Patents

Systeme de filtre median adaptif Download PDF

Info

Publication number
FR2596221A1
FR2596221A1 FR8703916A FR8703916A FR2596221A1 FR 2596221 A1 FR2596221 A1 FR 2596221A1 FR 8703916 A FR8703916 A FR 8703916A FR 8703916 A FR8703916 A FR 8703916A FR 2596221 A1 FR2596221 A1 FR 2596221A1
Authority
FR
France
Prior art keywords
samples
signal
sample
noise
control signal
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
FR8703916A
Other languages
English (en)
Other versions
FR2596221B1 (fr
Inventor
Stuart Stanley Perlman
Sanford Eisenhandler
Paul Wallace Lyons
Michael John Shumila
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.)
RCA Corp
Original Assignee
RCA Corp
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 RCA Corp filed Critical RCA Corp
Publication of FR2596221A1 publication Critical patent/FR2596221A1/fr
Application granted granted Critical
Publication of FR2596221B1 publication Critical patent/FR2596221B1/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/14Picture signal circuitry for video frequency region
    • H04N5/21Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20004Adaptive image processing
    • G06T2207/20012Locally adaptive
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • G06T2207/20032Median filtering

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Filters That Use Time-Delay Elements (AREA)
  • Picture Signal Circuits (AREA)
  • Noise Elimination (AREA)

Abstract

L'INVENTION CONCERNE UN SYSTEME DE FILTRE ADAPTIF MEDIAN. SELON L'INVENTION, IL COMPREND UNE SOURCE D'UN SIGNAL D'ENTREE 5 CONTENANT EVENTUELLEMENT DU BRUIT, UN MOYEN 10 REPONDANT A CE SIGNAL D'ENTREE POUR PRODUIRE DES GROUPES SUCCESSIFS D'ECHANTILLONS LE REPRESENTANT; UN MOYEN 20 POUR LE FILTRAGE MEDIAN ADAPTIF DES ECHANTILLONS EN REPONSE A UN SIGNAL DE COMMANDE; UN MOYEN 30 REPONDANT AU SIGNAL D'ENTREE POUR ESTIMER LA DENSITE RELATIVE DU BRUIT ET UN MOYEN 40 REPONDANT A LA DENSITE ESTIMEE DE BRUIT POUR PRODUIRE LE SIGNAL DE COMMANDE. L'INVENTION S'APPLIQUE NOTAMMENT AU TRAITEMENT DE L'IMAGE.

Description

La présente invention se rapporte à des systèmes
de filtre médian adaptif.
Les développements du traitement de l'image ont rendu possible le traitement d'images sous une grande variété de formes. Par exemple, une photographie peut être traitée sous la forme d'une série bidimensionnelle d'échantillons. Un signal vidéo, par exemple provenant d'une caméra de télévision, peut être traité sous la forme d'une série tridimensionnelle d'échantillons, la troisième 10 dimension étant la dimension temporelle et représentant le mouvement dans l'image. Un algorithme de filtrage que l'on peut appliquer au traitement de telles images est-le
filtre médian.
Un filtre médian réduit le bruit du type impulsion15 nel dans tout signal mais peut être particulièrement efficace pour réduire un bruit dans une image tel qu'il peut être provoqué par une rayure sur une photographie ou une chute du signal dans le canal de communication par lequel le signal vidéo est transmis. Cependant, un filtre 20 médian n'affecte pas de manière néfaste les bords ou coins qui existent dans une image. Dans des filtres médians de l'art antérieur, l'échantillon courant (l'échantillon qui est filtré) était remplacé par un échantillon ayant la valeur médiane des valeurs de cet échantillon et d'un
certain nombre d'échantillons entourant cet échantillon.
Une amélioration de l'algorithme de filtrage médian décrit ci-dessus substitue sélectivement l'échantillon ayant la valeur médiane à l'échantillon courant, uniquement si la valeur de l'échantillon courant se trouve en dehors 30 de seuils prédéterminés. Par exemple, dans l'article "Image Enhancement Using the Median and Interquartile Distance" de Scollar et autres, dans Computervision, Graphics and Image Processing, Volume 25, pages 236-251 (1984), la différence entre la valeur de l'échantillon 35 courant et la valeur médiane est comparée àuneconstante multipliée par la différence entre les quartiles supérieur et inférieur (la distance entre quartiles). Si la différence échantillon courant/médiane est plus importante que la fonction de la distance entre quartiles, alors l'échantillon ayant la valeur médiane est substitué à l'échantillon courant. Scollar et autres ont trouvé cet algorithme avantageux parce qu'il "s'adapte au contraste local dans les données, et que les bords nets seront à
peine brouillés,...".
Selon les principes de la présente invention, une autre amélioration peut être obtenue si les critères 10 se rapportant à la non image sont utilisés pour adapter l'algorithme du filtre médian, en plus des critères se rapportant au contraste décrits dans l'article de Scollar ci-dessus. Comme l'un des avantages principaux d'un filtre médian est sa capacité à atténuer le bruit impulsionnel, un meilleur fonctionnement peut être obtenu si le fonctionnement du filtre médian adaptif est modifié enréponse à la densité relative du bruit impulsionnel dans le signal
qui est filtré.
Un appareil selon les principes de la présente
invention comprend une source d'un signal d'entrée contenant éventuellement du bruit, et un moyen pour produire.
un certain nombre d'échantillonsreprésentant le signal.
Un filtre médian adaptif filtre ces échantillons en réponse à un signal de commande. Un moyen d'estimation de bruit 25 estime la densité relative du bruit dans le signal qui est échantillonné. Un générateur de signaux de commande produit le signal de commande appliqué au filtre adaptif
médian en réponse à la densité estimée du.bruit.
L'invention sera mieux comprise, et d'autres buts, caractéristiques, détails et avantages de celle-ci
apparaîtront plus clairement au cours de la description
explicative qui va suivre faite en référence aux dessins schématiques annexés donnés uniquement à titre d'exemple illustrant plusieurs modes de réalisation de l'invention 35 et dans lesquels: - la figure l donne un schéma-bloc d'un système de filtre médian adaptif selon les principes de la présente invention; - la figure 2 donne un schéma-bloc d'un filtre médian adaptif que l'on peut utiliser dans le système illustré à la figure 1; - la figure 3 donne un schéma-bloc d'un producteur de M- tilesque l'on peut utiliser dans le filtre médian illustré à la figure 2; et - la figure 4 donne un schéma-bloc d'un autre
mode de réalisation d'un système de filtre médian adaptif 10 selon les principes de la présente invention.
Sur la figure 1, le signal d'entrée, passant par le filtre médian,.par exemple un signal vidéo de l'extrémité avant d'un récepteur de télévision (non représenté) est appliqué à une borne d'entrée 5. La borne d'entrée 5 15 est reliée à un moyen 10 de production d'échantillons et un moyen 30 d'estimation du bruit. Le moyen 10 de production d'échantillons produit un certain nombre d'échantillons, comprenant un échantillon courant (l'échantillon qui est filtré à ce moment) représentant un point donné dans l'image qui est traitée et d'autres échantillons représentant des points entourant ce point. Ces échantillons sont appliqués aux bornes respectives d'entrée de données d'un filtre médian adaptif 20. Le filtre 20 produit, à une borne de sortie, un signal filtré de sortie. Une borne de sortie du moyen d'estimation 30 du bruit est coupléeà une borne d'entrée d'un générateur de signaux de commande qui produit un signal de commande à une borne de sortie, lequel est appliqué à une borne d'entrée de commande du
filtre médian adaptif 20.
En fonctionnement, le moyen 10 de production d'échantillons produit des versions du signal d'entrée à la borne 5, qui représentent les points courant et environnants. Dans un système dans lequel le signal d'entrée est un signal vidéo, le moyen 10 de production d'échantil35 lons peut comprendre une ligne à retard à-prises, qui produit un échantillon courant et des échantillons entourant l'échantillon courant soit verticalement, horizontalement, temporellement ou toute combinaison de cela. Le filtre médian adaptif 20 applique sélectivement soit l'échantillon courant ou un échantillon ayant la valeur médiane des échantillons courant et environnants à la borne de sortie du filtre. La décision de substitution d'un échantillon de valeur médiane à l'échantillon courant est basée sur les valeurs des échantillons d'entrée et du signal de
commande appliqué au filtre médian adaptif 20.
Le moyen 30 d'estimation du bruit évalue le signal 10 d'entrée appliqué à la borne d'entrée 5 pour déterminer la densité relative du bruit dans le signal d'entrée. Comme le filtre médian adaptif est particulièrement utile pour atténuer le bruit impulsionnel, le moyen 30 d'estimation du bruit peut, par exemple, détecter la densité relative du bruit impulsionnel dans le signal d'entrée. Le moyen 30 peut comprendre un détecteur de seuils qui produit une impulsion à chaque fois que le signal d'entrée dépasse les seuils, indiquant une impulsion de bruit; et un compteur
pouvant être remis à zéro, qui compte le nombre d'impulsions 20 sur un intervalle prédéterminé de temps.
Le signal de commande est produit par le générateur de signaux de commande, en réponse à la-densité du bruit impulsionnel pour atténuer plus complètement le bruit impulsionnel sans affecter de manière néfaste les renseigne25 ments transmis dans le signal. Le générateur 40 de signaux de commande peut par exemple être un décodeur pour produire
des signaux de commande ayant des valeurs différentes qui correspondent à des plages différentes de densité estimée du bruit. Dans une forme numérique, ce décodeur peut être 30 une mémoire numérique préprogrammée.
La figure 2 illustre un filtre médian adaptif 20 que l'on peut utiliser dans l'appareil illustré à la figure 1. Sur la figure 2, les échantillons du moyen 10 de production d'échantillons de la figure 1 sont appliqués aux bornes d'entrée de données d'un moyen 21 de production de M-tiles. Le signal de commande du générateur 40 est appliqué à une borne d'entrée de commande du moyen de production 21 de M-tiles. Le moyen de production 21 produit un échantillon de M-tile supérieur, un échantillon médian et un échantillon de M-tile inférieur à trois bornes de sortie de données. Les échantillons de M-tiles supérieur et inférieur sont appliqués aux bornes respectives d'entrée d'un soustracteur 23. La borne de sortie du soustracteur 23 est couplée à une première borne d'entrée d'un comparateur 27 par un générateur de fonction 26. L'échantillon courant (à l'entrée du moyen de production 21 de M-tiles) et l'échantillon médian de sortie sont appliqués aux bornes respectives d'entrée d'un soustracteur 25. L.a borne de sortie du soustracteur 25 est couplée à une seconde borne d'entrée du comparateur 27 par un générateur 28 de valeurs absolues. Une borne de sortie du comparateur 27 est
couplée à une borne d'entrée de commande du multiplexeur 29.
Une première borne d'entrée de données du multiplexeur 29 est couplée pour recevoir l'échantillon médian du moyen 21 de production de M-tiles. Une seconde borne d'entrée de données du multiplexeur 29 est couplée pour recevoir l'échantillon courant de l'entrée du moyen de production 21 de Mtiles. La borne de sortie de données du multiplexeur
29 produit l'échantillon sorti médian filtré.
M-tile est un terme général se rapportant à la position relative d'un échantillon dans une liste d'échan25 tillons triés selon leur valeur. La médiane et les quartiles supérieur et inférieur sont des cas spéciaux indiquant des valeurs à la moitié, les trois quarts et
un quart du chemin à travers la liste ordonnée respectivement.
En fonctionnement, le moyen 21 de production de M-tiles trie les échantillons d'entrée pour produire une liste ordonnée d'échantillons. De la liste ordonnée, l'échantillon du milieu a la valeur médiane et est produit à la borne de sortie médian. Les échantillons de deux 35 autres positions dans la liste ordonnée sont produits aux bornes de sortie des Mtiles supérieur et inférieur. Le signal de commande contrôle les positions particulières de la liste ordonnée qui sont choisies pour les M- tiles
supérieur et inférieur.
Dans un article intitulé "The LSI Architecture for a One-Chipe. Video Median Filter" de Demassieux et autres, dans Proceedings IEEE International Conference on Acoustics, Speech and Signal Processing, Volume 3, pages 1001-1004 (1985), est décrit un exemple d'une réalisation numérique d'un filtre médian. Le filtre comprend un élément de triage comprenant une série de dix comparateurs à 8 bits 10 qui trient cinq échantillons de 8 bits en une liste ordonnée d'o une valeur médiane et une valeur de M-tile supérieur et une valeur de M-tile inférieur peuvent être
prélevées à chaque cycle d'horloge.
Le soustracteur 23 soustrait la valeur du M-tile 15 inférieur de la valeur du M-tile supérieur. La différence est la distance entre M-tiles. Une fonction de la distance entre M-tiles est calculée dans le générateur de fonctions 26. Dans l'article ci-dessus mentionné de Scollar et autres,
par exemple, la distance entre quartiles est multipliée 20 par une constante.
D'autres fonctions de la distance entre M-tiles peuvent également être avantageuses. Par exemple, une valeur fixe peut être produite, ou bien une fonction linéaire de la distance entre M-tiles. Une telle fonction 25 est égale à la moitié de la distance entre M-tiles plus une constante. Dans une forme numérique, la fonction peut être préprogrammée dans une mémoire morte programmable (PROM), ayant des bornes d'entrée d'adresse qui sont
couplées à la sortie du soustracteur 23 et des bornes de 30 sortie de données qui sont couplées au comparateur 27.
La fonction de la distance entre M-tiles peut également être rendue adaptive en réponse au signal de commande comme cela est illustré en tracé fantôme sur la figure 2. La fonction particulière de la distance entre
M-tiles choisie dépend de la valeur du signal de commande.
Dans une forme numérique comprenant une PROM, le signal de commande peut être appliqué à d'autres bornes d'entrée d'adresse de la PROM dans laquelle toutes les fonctions
possibles ont été programmées.
La sortie du soustracteur 25 est la différence entre la médiane et l'échantillon courant. Le générateur 28 de valeurs absolues produit la grandeur de la différence qui est appelée distance échantillon courant/médian Cette distance est comparée à la valeur de la fonction de la distance entre M-tiles par le comparateur 27. Si la distance échantillon courant/médian dépasse la distance 10 entre M-tiles traitée, alors le comparateur 27 produit un signal, Co, ayant un premier état et autrement le signal Co a un second état. Si le signal Co est dans son premier état, le multiplexeur 29 applique l'échantillon de valeur médiane du moyen 21 de production de M-tiles à la borne 15 de sortie du filtre; si le signal Co est à son second état, le multiplexeur 29 applique l'échantillon courant à
la borne de sortie du filtre. -
Ce système peut être configuré comme un filtre
médian standard dans lequel chaque échantillon est remplacé 20 par un échantillon ayant la valeur médiane correspondante.
Si la sortie du générateur de fonctions 26 a une valeur si faible que chaque différence échantillon courant/médian dépasse toujours cette valeur (comme 0), alors l'échantillon de valeur médiane remplacera toujours l'échantillon 2' courant. Inversement, le filtrage médian peut être arrêté si le générateur de fonctions 26 produit une valeur si importante qu'il n'y a jamais de différence échantillon courant/médian dépassant cette valeur. Dans ce cas,
l'échantillon courant passera toujours par le multiplexeur 30 29.
La distance entre M-tiles est une mesure du contraste de l'image dans la localité de l'échantillon courant. Si le contraste est fort, la distance entre M-tiles sera importante et inversement. S'il y a un fort contraste 35 dans la localité de l'échantillon courant, alors l'échantillon courant doit s'écarter relativement fortement de la valeur médiane avant d'être remplacé par la valeur médiane. Par ailleurs, si le contraste est relativement faible dans la localité de l'échantillon courant, alors l'échantillon courant ne doit varier de la médiane que d'une relativement faible quantité avant d'être remplacé par la valeur médiane. Par conséquent, un tel filtre médian est adaptif et sensible au contraste local de
l'image dans la zone de l'échantillon courant.
Dans la présente demande, le signal à la sortie du générateur 40 de signaux de commande de la figure 1 10 forcera le moyen 21 de production de M-tiles à produire des valeurs différentes pour les M-tiles supérieur et inférieur. Par exemple, si la densité du bruit est relativement forte, alors les M-tiles supérieur et inférieur peuvent être choisis pour être relativement proches de 15 la médiane. La distance entre M-tiles est ainqi assez faible et l'échantillon courant doit seulement s'écarter de la valeur médiane d'une relativement faible quantité
avant d'être remplacé par l'échantillon de valeur médiane.
Cela aura tendance à remplacer plus d'échantillons par la 20 médiane en réponse à la densité supérieure de bruit. Par ailleurs, si la densité du bruit est relativement faible', alors les M-tiles supérieur et inférieur peuvent être
choisis pour être relativement plus loin de la médiane.
Dans ce cas, la distance entre M-tiles est plus importante, 25 et l'échantillon courant doit s'écarter de la valeur médiane d'une relativement grande quantité avant d'être remplacé par l'échantillon de valeur médiane. Cela aura tendance à remplacer moins d'échantillons par la médiane en réponse à la plus faible densité du bruit. De plus, si 30 la densité du bruit est très faible, il peut être avantageux d'arrêter complètement le filtre médian par un choix approprié de la fonction utilisée par le générateur de
fonctions 26, comme on l'a décrit ci-dessus.
La figure 3 illustre un mode de réalisation possible du moyen 21 de production de M-tiles illustré à la figure 2. Sur la figure 3, les échantillons d'entrée sont appliqués aux bornes respectives d'entrée d'un élément trieur 22. L'élément trieur 22 peut être tel que décrit dans l'article ci-dessus mentionné de Demassieux et autres. L'élément trieur 22 produit à ses bornes respectives de sortie, la liste ordonnée des valeurs des 5 échantillons d'entrée. Par exemple, un échantillon ayant la plus grande valeur des échantillons d'entrée, est produit à la borne de sortie la plus haute et un
Échantillon ayant la plus faible valeur des échantillons d'entrée est produit à la borne de sortie la plus basse.
Un circuit de multiplexage (MUX) 24 a ses bornes respectives d'entrée couplées aux bornes de sortie de l'élément trieur 22. Dans le circuit de multiplexage 24, l'échantillon
du milieu à la sortie de l'élément trieur 22 est appliqué à la borne de sortie médiane du circuit de démultiplexage 15 24 (comme cela est illustré par la ligne en pointillé).
La borne de sortie du M-tile supérieur du circuit de multiplexage 24 est couplée à l'une des bornes d'entrée du circuit 24 en réponse au signal de commande. La borne de sortie du M-tile inférieur du circuit de multiplexage 24 20 est couplée à une borne différente d'entrée du circuit de
multiplexage 24, également selon le signal de commande.
Le circuit de multiplexage 24 peut,par exemple, comprendre une connexion directe de sa borne d'entrée médiane à sa borne de sortie médiane et deux multiplexeurs, 25 un pour le couplage de chaque borne de sortie de M-tile supérieur et de M-tile inférieur respectivement à des bornes d'entrée choisies du circuit de multiplexage 24, chacun recevant le signal de commande de la borne d'entrée
de commande.
En se référant de nouveau à la figure 2, le comparateur 27 produit un signal Co à un "1" logique à chaque fois que l'échantillon courant peut être remplacé par l'échantillon médian. Co à la sortie du comparateur 27
est ainsi une indication que l'échantillon courant peut 35 représenter une impulsion de bruit qu'il faut remplacer.
Ce signal peut être utilisé pour estimer la densité relative du bruit impulsionnel dans le signal d'entrée en déterminant le nombre de présences du bruit impulsionnel
dans un intervalle de temps. Le signal Co est montré en pointillé appliqué à une borne de sortie de la figure 2.
La figure 4 montre un autre mode de réalisation d'un filtre médian adaptif selon les principes de la présente invention. Des éléments similaires à ceux de la figure 1 sont désignés par des chiffres identiques et fonctionnent d'une manière semblable. Sur la figure 4, un signal d'entrée de la borne 5 est appliqué à un moyen 10 10 de production d'échantillons.:Le moyen 10 applique un certain nombre d'échantillons à un bloc de filtre médian adaptif 20'. Le filtre médian adaptif 20' produit, à ses bornes respectives de sortie, un signal médian filtré de sortie et un signal Co qui, comme on l'a décrit ci-dessus, 15 indique qu'une impulsion de bruit a été détectée. Le signal Co est appliqué à une borne d'entrée d'un moyen 50 d'estimation du bruit. Une seconde borne d'entrée du moyen 50 est couplée à une borne d'entrée de signaux de synchronisation. Le moyen 50 d'estimation du bruit produit, à une borne de sortie, une estimation de la densité relative du bruit impulsionnel détectée dans le signal d'entrée. Ce signal d'estimation est appliqué à la borne d'entrée du générateur 40 de signaux de commande
qui produit un signal de commande appliqué à la borne 25 d'entrée de commande du filtre médian adaptif 20'.
Si le signal à la borne d'entrée 5 représente, par exemple, un signal vidéo composite, alors il y a des intervalles de temps dans le signal vidéo,pendant l'intervalle de synchronisation verticale o le signal est à un niveau prédéterminé. Le signal de synchronisation appliqué au moyen 50 d'estimation du bruit peut par exemple être la composante de synchronisation verticale du signal vidéo composite. Pendant l'intervalle de temps o l'on sait que le signal est au niveau prédéterminé, toutes les impulsions 35 indiquées par le signal Co sont supposées être provoquées par le bruit. Le nombre d'impulsions indiquées par le signal Co à la sortie du filtre adaptif médian 20' pendant l'intervalle de synchronisation verticale est compté dans le moyen 50 d'estimation du bruit. Plus le compte
est élevé, plus la densité relative du bruit est importante.
Le moyen 50 d'estimation du bruit peut par exemple comprendre un compteur numérique déclenché. Au début de l'intervalle de synchronisation verticale, le compteur est remis à zéro et le comptage est validé. A la fin de l'intervalle de synchronisation verticale, le comptage
est inhibé et le compte résultant est verrouillé à la 10 borne de sortie du moyen 50 d'estimation du bruit.
Les résultats du comptage des impulsions pendant l'intervalle de synchronisation verticale sont appliqués en tant que signal d'entrée au générateur 40 de signaux de commande qui produit alors le signal de commande pour 15 contrôler les valeurs des M-tiles supérieur et inférieur que l'on utilise pour déterminer les seuils utilisés dans le filtre médian adaptif. Comme on l'a décrit ci-dessus, pour une forte densité du bruit, le signal de commande est produit de manière que les M-tiles soient choisis 20 relativement prochesde la médiane et pour une faible densité du bruit, ils sont choisis relativement plus loin
de la médiane.
Il faut noter que l'agencement illustré sur la figure 4 peut être utilisé pour estimer la densité du bruit 25 impulsionnel sans accomplir aucun filtrage. Un tel appareil comprendra, en tant que sortie, la densité estimée du bruit impulsionnel à la sortie du moyen 50 d'estimation du bruit, comme cela est montré en tracé fantôme à la figure 4. Il peut être souhaitable d'estimer la densité du bruit 30 impulsionnel en tout moment. Dans un tel cas, aucune entrée de signaux de synchronisation ne sera requise pour le moyen 50 d'estimation du bruit. Il peut de plus être souhaitable de produire une indication qu'une impulsion est présente. Dans ce cas, le signal Co du filtre médian 35 adaptif 20' est produit à une borne de sortie, ce qui
est également montré en tracé fantôme.
Les modes de réalisation ci-dessus ont été décrits en référence à un signai vidéo composite. Il faut cependant noter que tout signal qui forme l'objet d'un bruit impulsionnel peut 6tre avantageusement filtré en utilisant le système du filtre médian adaptif décrit ci-dessus. L'invention peut également être mise en oeuvre dans un système de données continues ou échantillonnées et le système de données échantillonnées peut être mis
en oeuvre sous forme analogique ou numérique.

Claims (8)

R E V E N D I C A T I 0 N S
1.- Système de filtre médian adaptif, caractérisé par: une source d'un signal d'entrée (5) contenant éventuellement du bruit; un moyen (10) répondant audit signal d'entrée pour produire des groupes successifs d'échantillons représentant ledit signal d'entrée; un moyen (20) pour le filtrage médian adaptif 10 desdits échantillons en réponse à un signal de commande; un moyen (30) répondant audit signal d'entrée pour estimer la densité relative du bruit; et un moyen (40) répondant à la densité estimée de
bruit pour produire le signal de commande.
2.- Système selon la revendication 1, caractérisé en ce que le moyen formant filtre médian adaptif précité comprend: un moyen pour déterminer un échantillon de valeur médiane et des échantillons de valeurs de Mtiles supérieur 20 et inférieur dudit groupe d'échantillons (21); un moyen répondant à un autre signal de commande pour sélectivement substituer ledit échantillon de valeur médiane à un échantillon courant dudit groupe d'échantillons qui est filtré; et un moyen (23, 25, 27) répondant audit échantillon de valeur médiane, auxdits échantillons de valeurs de M-tiles supérieur et inférieur et audit échantillon courant
pour produire ledit autre signal de commande.
3.- Système selon la revendication 2, caractérisé 30 en ce que ledit autre générateur de signaux de commande comprend: un premier moyen de soustraction (23) répondant aux échantillons des valeurs de M-tiles supérieur et inférieur pour produire la différence entre eux; un second moyen de soustraction (25) répondant audit échantillon courant et audit échantillon à la valeur médiane pour produire la différence entre eux; et un moyen de comparaison (27) relié auxdits premier et second moyens de soustraction pour comparer les différences desdits premier et second moyens de soustraction et produire ledit autre signal de commande si la différence dudit second moyen de soustraction
dépasse la différence dudit premier moyen de soustraction.
4.- Système selon la revendication 3, caractérisé 10 en ce que l'autre générateur de signaux de commande précité comprend: un générateur de fonctions (26) couplé entre le
premier moyen de soustraction et le moyen de comparaison.
5.- Système selon la revendication 4, caractérisé 15 en ce que le générateur de fonctions (26) comprend un moyen pour multiplier un signal à l'entrée du'générateur de fonctions pour une première constante et ajouter une
seconde constante non nulle.
6.- Système selon la revendication 2, caractérisé 20 en ce que le moyen de production de valeurs médianes et d'échantillons de valeurs de M-tiles supérieur et inférieur comprend: un moyen (22) ayant un certain nombre de bornes d'entrée répondant au groupe d'échantillons et un certain 25 nombre de bornes de sortie pour produire le groupe d'échantillons ordonnés selon leurs valeurs; et un moyen de multiplexage (24) ayant des bornes respectives d'entrée de données répondant au groupe ordonné d'échantillons, des bornes de sortie de données pour produire au moins des échantillons des valeurs des M-tiles supérieur et inférieur et une borne d'entrée de commande répondant audit signal de commande, pour coupler certains des échantillons du groupe ordonné aux bornes de sortie de
données d'échantillons de valeurs de M-tiles supérieur et 35 inférieur en réponse à la valeur du signal de commande.
7.- Système de filtre médian adaptif, caractérisé par: une source (5) d'un signal d'entrée contenant éventuellement du bruit; un moyen (10) répondant au signal d'entrée pour produire des groupes successifs d'échantillons représentant le signal d'entrée; un moyen (20') pour le filtrage médian adaptif des échantillons en réponse à un signal de commande et pour produire un signal indiquant qu'un échantillon courant du groupe d'échantillons représente un signal 10 relativement bruyant; un moyen (50) répondant au signal d'indication pour estimer la densité relative du bruit; et un moyen (40) répondant à la densité estimée
du bruit pour produire le signal de commande.
8.- Système de filtre médian adaptif, caractérisé par: une source (5) d'un signal d'entrée contenant des intervalles de temps alternés de synchronisation et d'information et contenant éventuellement du bruit; un moyen (10) répondant au signal d'entrée pour produire des groupes successifs d'échantillons représentant le signal d'entrée; un moyen (20') pour le filtrage médian adaptif des échantillons en réponse à un signal de commande pendant l'intervalle de temps de l'information et pour produire un signal indiquant qu'un échantillon courant du groupe d'échantillons représente un signal relativement bruyant pendant l'intervalle de temps de synchronisation; un moyen (50) répondant au signal d'indication 30 pour estimer la densité relative du bruit; et un moyen (40) répondant à la densité estimée
de bruit pour produire le signal de commande.
FR8703916A 1986-03-21 1987-03-20 Systeme de filtre median adaptif Expired - Lifetime FR2596221B1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US06/842,644 US4682230A (en) 1986-03-21 1986-03-21 Adaptive median filter system

Publications (2)

Publication Number Publication Date
FR2596221A1 true FR2596221A1 (fr) 1987-09-25
FR2596221B1 FR2596221B1 (fr) 1990-12-21

Family

ID=25287894

Family Applications (1)

Application Number Title Priority Date Filing Date
FR8703916A Expired - Lifetime FR2596221B1 (fr) 1986-03-21 1987-03-20 Systeme de filtre median adaptif

Country Status (7)

Country Link
US (1) US4682230A (fr)
JP (1) JPH0783236B2 (fr)
DE (1) DE3709068C2 (fr)
FR (1) FR2596221B1 (fr)
GB (1) GB2188511B (fr)
HK (1) HK146894A (fr)
SG (1) SG25092G (fr)

Families Citing this family (60)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4761829A (en) * 1985-11-27 1988-08-02 Motorola Inc. Adaptive signal strength and/or ambient noise driven audio shaping system
US4933978A (en) * 1986-02-28 1990-06-12 Rca Licensing Corporation Method and apparatus for determining the value of a sample in the mth position of an ordered list of a plurality of samples
US4682230A (en) 1986-03-21 1987-07-21 Rca Corporation Adaptive median filter system
US4782389A (en) * 1987-04-30 1988-11-01 Rca Licensing Corporation Adaptive M-tile sample producer
US4774565A (en) * 1987-08-24 1988-09-27 Polaroid Corporation Method and apparatus for reconstructing missing color samples
DE3884212T2 (de) * 1987-10-21 1994-01-13 James Andrew Bangham Filterungstechniken.
JPH01136418A (ja) * 1987-11-21 1989-05-29 Nec Corp 搬送波位相制御ループ
JPH01260913A (ja) * 1988-04-12 1989-10-18 Toshiba Corp デジタルフィルタ
US4933766A (en) * 1988-06-02 1990-06-12 U.S. Philips Corporation Interpolation filter and receiver provided with such an interpolation filter
JP2870762B2 (ja) * 1988-07-19 1999-03-17 ソニー株式会社 画像信号の高能率符号化装置
US5166794A (en) * 1989-03-10 1992-11-24 Sony Corporation Drop out compensating circuit with interpolation of proximate pixel signals
US5043814A (en) * 1989-10-06 1991-08-27 North American Philips Corporation Adaptive ghost cancellation circuit
DE4001552A1 (de) * 1990-01-20 1991-07-25 Broadcast Television Syst Schaltungsanordnung zur filterung von videosignalen
US5218649A (en) * 1990-05-04 1993-06-08 U S West Advanced Technologies, Inc. Image enhancement system
DE4028322A1 (de) * 1990-09-06 1992-03-19 Bosch Gmbh Robert Verfahren zur verbesserung einer kanalparameterschaetzung und anordnung zu dessen durchfuehrung
DE4031785C2 (de) * 1990-10-08 1996-07-04 Broadcast Television Syst Verfahren zur Reduktion von Rauschen in Videosignalen
US5161014A (en) * 1990-11-26 1992-11-03 Rca Thomson Licensing Corporation Neural networks as for video signal processing
US5517531A (en) * 1991-05-29 1996-05-14 The United States Of America As Represented By The Secretary Of The Navy Kernel adaptive interference suppression system
GB2256991A (en) * 1991-06-07 1992-12-23 Clarion Co Ltd Impulse noise eliminating circuit for a television receiver
FI90176C (fi) * 1991-09-02 1993-12-27 Salora Oy Metod foer att filtrera drop-out -brus
US5327240A (en) * 1991-12-24 1994-07-05 Texas Instruments Incorporated Methods, systems and apparatus for providing improved definition video
FI90477C (fi) * 1992-03-23 1994-02-10 Nokia Mobile Phones Ltd Puhesignaalin laadun parannusmenetelmä lineaarista ennustusta käyttävään koodausjärjestelmään
SE469411B (sv) * 1992-04-13 1993-06-28 Dv Sweden Ab Saett att detektera och avlaegsna fel oeverstigande en specifik kontrast i digitala videosignaler
US5384865A (en) * 1992-06-01 1995-01-24 Eastman Kodak Company Adaptive, hybrid median filter for temporal noise suppression
EP0574969B1 (fr) * 1992-06-19 2002-03-20 Agfa-Gevaert Méthode et dispositif de réduction du bruit
FI93068C (fi) * 1992-07-27 1995-02-10 Nokia Mobile Phones Ltd Kytkentä häiriöiden haittavaikutusten pienentämiseksi sovitettua suodatinta käyttävissä vastaanottimissa
FI91579C (fi) * 1992-08-20 1994-07-11 Nokia Mobile Phones Ltd Dekoodaus käyttäen lineaarista metriciä ja häiriön estimointia
DE69329670T2 (de) * 1992-09-14 2001-03-15 Thomson Multimedia Sa Verfahren und Gerät zur Rauschminderung
US5446501A (en) * 1992-10-22 1995-08-29 Accom, Incorporated Three-dimensional median and recursive filtering apparatus and method for video image enhancement
US5361105A (en) * 1993-03-05 1994-11-01 Matsushita Electric Corporation Of America Noise reduction system using multi-frame motion estimation, outlier rejection and trajectory correction
US5412436A (en) * 1993-04-22 1995-05-02 Thomson Consumer Electronics, Inc. Motion adaptive video processing system
KR970009469B1 (ko) * 1993-06-05 1997-06-13 삼성전자 주식회사 더블스무딩(Double Smoothing) 기능을 갖는 비월/순차주사변환장치 및 그 방법
DE4405723A1 (de) * 1994-02-23 1995-08-24 Daimler Benz Ag Verfahren zur Geräuschreduktion eines gestörten Sprachsignals
JP2949186B2 (ja) * 1994-03-18 1999-09-13 富士通株式会社 画像処理方法及び画像処理装置
US5798846A (en) * 1996-05-28 1998-08-25 Hewlett-Packard Company Apparatus and method for selectively processing a scanned image
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
US5995163A (en) * 1996-09-30 1999-11-30 Photobit Corporation Median filter with embedded analog to digital converter
DE19743124B4 (de) * 1997-09-30 2004-08-26 Harman Becker Automotive Systems (Xsys Division) Gmbh Verfahren zur Bestimmung der Qualität eines Videosignals und/oder eines Fernsehbildes
DE19743125A1 (de) * 1997-09-30 1999-04-08 Ise Interactive Systems Entwic Verfahren zur Bewertung der Qualität eines Fernsehbildes
JPH11113030A (ja) 1997-10-07 1999-04-23 Pioneer Electron Corp コンポジットビデオ信号のsn比測定方法及びシステム
US6058405A (en) * 1997-11-06 2000-05-02 Motorola Inc. SIMD computation of rank based filters for M×N grids
SE515674C2 (sv) 1997-12-05 2001-09-24 Ericsson Telefon Ab L M Apparat och metod för brusreducering
US6654510B1 (en) 1998-10-14 2003-11-25 Konica Corporation Image processing method and image processing system
JP4147647B2 (ja) 1998-11-09 2008-09-10 ソニー株式会社 データ処理装置およびデータ処理方法、並びに記録媒体
JP4344964B2 (ja) * 1999-06-01 2009-10-14 ソニー株式会社 画像処理装置および画像処理方法
US6760737B2 (en) * 2000-03-27 2004-07-06 Lucent Technologies Inc. Spatial median filter
US7133568B2 (en) 2000-08-04 2006-11-07 Nikitin Alexei V Method and apparatus for analysis of variables
DE10138353A1 (de) * 2001-08-04 2002-02-28 Grundig Ag Verfahren und Vorrichtung bei einer pulsbreitengesteuerten Bildanzeigevorrichtung
US7092573B2 (en) * 2001-12-10 2006-08-15 Eastman Kodak Company Method and system for selectively applying enhancement to an image
US6934421B2 (en) * 2002-03-20 2005-08-23 Eastman Kodak Company Calculating noise from multiple digital images having a common noise source
US6958783B2 (en) * 2002-07-16 2005-10-25 Broadcom Corporation Adaptive non-linear noise reduction techniques
US7050501B2 (en) * 2002-07-16 2006-05-23 Broadcom Corporation Digital noise reduction techniques
US7617270B2 (en) * 2002-10-07 2009-11-10 Nikitin Alexei V Method and apparatus for adaptive real-time signal conditioning, processing, analysis, quantification, comparison, and control
US7107306B2 (en) * 2002-10-07 2006-09-12 Nikitin Alexei V Method and apparatus for adaptive real-time signal conditioning, processing, analysis, quantification, comparision, and control
KR100483697B1 (ko) * 2002-12-06 2005-04-18 엘지전자 주식회사 비디오 데이터 처리장치 및 방법
JP4144563B2 (ja) * 2003-07-11 2008-09-03 ティアック株式会社 ディスク装置
KR100522607B1 (ko) * 2003-07-15 2005-10-19 삼성전자주식회사 노이즈 상태를 고려한 적응적 비디오 신호 처리 장치 및방법
US7477319B2 (en) * 2005-06-17 2009-01-13 Lsi Corporation Systems and methods for deinterlacing video signals
US20150173756A1 (en) 2013-12-23 2015-06-25 Ethicon Endo-Surgery, Inc. Surgical cutting and stapling methods
CN108416740B (zh) * 2018-01-22 2021-09-14 大连大学 一种消除椒盐噪声的迭代自适应中值滤波方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4177430A (en) * 1978-03-06 1979-12-04 Rockwell International Corporation Adaptive noise cancelling receiver
US4376955A (en) * 1980-02-28 1983-03-15 Rca Corporation Two dimensional adaptive dropout compensator and chroma inverter
CA1157525A (fr) * 1980-03-08 1983-11-22 Majid Sepehri Recepteur pour systeme de radiocommunication en diversite de frequence comprenant un circuit d'evaluation d'interference et un filtre d'adaptation
US4413282A (en) * 1981-04-06 1983-11-01 Rca Corporation Television ghost signal detection during the video information interval
DE3114275A1 (de) * 1981-04-09 1982-11-04 Robert Bosch Gmbh, 7000 Stuttgart Verfahren und schaltungsanordnung zum verdecken von fehlern in einem digitalen videosignal
US4517600A (en) * 1982-05-17 1985-05-14 Rca Corporation Error concealment system using dropout severity information for selection of concealment method
US4562470A (en) * 1983-04-07 1985-12-31 Centre De Recherche Industrielle Du Quebec Electrical noise reducer for T.V. signals
JPS6066544A (ja) * 1983-09-21 1985-04-16 Toshihiko Namekawa インパルス性雑音に対する最適受信機
US4682230A (en) 1986-03-21 1987-07-21 Rca Corporation Adaptive median filter system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
COMMUNICATIONS OF THE ASSOCIATION OF COMPUTING MACHINERY, vol. 27, no. 8, août 1984, pages 807-818, New York, US; D.R.K. BROWNRIGG: "The weighted median filter" *
IEEE TRANSACTIONS ON ACOUSTICS, SPEECH AND SIGNAL PROCESSING, vol. ASSP-32, no. 3, juin 1984, pages 571-576, IEEE, New York, US; C.A. POMALAZA-RAEZ et al.: "An adaptative, nonlinear edge-preserving filter" *
PATTERN RECOGNITION LETTERS, vol. 2, no. 2, décembre 1983, pages 89-93, Elsevier Science Publishers B.V. (North-Holland), Amsterdam, NL; H.H.S. IP et al.: "Impulse noise cleaning by iterative threshold median filtering" *

Also Published As

Publication number Publication date
DE3709068C2 (de) 1995-10-12
HK146894A (en) 1994-12-30
FR2596221B1 (fr) 1990-12-21
JPH0783236B2 (ja) 1995-09-06
GB2188511A (en) 1987-09-30
US4682230A (en) 1987-07-21
SG25092G (en) 1992-05-15
GB8706706D0 (en) 1987-04-23
DE3709068A1 (de) 1987-09-24
JPS62231509A (ja) 1987-10-12
GB2188511B (en) 1989-11-29

Similar Documents

Publication Publication Date Title
FR2596221A1 (fr) Systeme de filtre median adaptif
FR2594281A1 (fr) Appareil d'estimation du bruit dans des signaux ayant des intervalles redondants
US6151415A (en) Auto-focusing algorithm using discrete wavelet transform
JP4601623B2 (ja) ノイズ除去装置
US20070025482A1 (en) Flexible sampling-rate encoder
FR2579395A1 (fr) Appareil pour controler de maniere adaptive un filtre recursif de signaux video
FR2766946A1 (fr) Procede et dispositif de pretraitement pour estimation de mouvement
FR2549316A1 (fr) Dispositif de creusement de signaux numeriques avec niveau reglable du seuil de creusement
FR2807902A1 (fr) Procede de detection de changement de plan dans une succession d'images video
FR2508261A1 (fr) Circuit de reduction de bruit pour signaux video
FR2557410A1 (fr) Circuit de traitement de signaux
FR2702914A1 (fr) Dispositif de codage de suites d'images constituées d'images de nature film et d'images de nature vidéo, et dispositif de décodage correspondant.
Bosco et al. Fast method for noise level estimation and integrated noise reduction
EP0342756A1 (fr) Système d'émission et réception de télévision à haute définition à débit de données réduit
EP1297710B1 (fr) Dispositif et procede de codage video avec filtrage recursif compense en mouvement
EP0506535A1 (fr) Procédé et système de traitement des pré-échos d'un signal audio-numérique codé par transformée fréquentielle
EP0400520B1 (fr) Dispositif de détection de perte de signal de réception pour récepteur de signaux numériques
FR2860941A1 (fr) Dispositif et procede d'estimation de bruit d'un signal video
EP0123573B1 (fr) Procédé de codage adaptatif, et de décodage, d'une image de télévision, et dispositifs pour la mise en oeuvre de ce procédé
EP1223554A1 (fr) Méthode de réduction de bruit
EP0690623B1 (fr) Procédé et dispositif d'insertion de données asynchrones sur un signal numérique
FR2474722A1 (fr) Dispositif de determination numerique de correlation d'images
CA3105372C (fr) Traitement d'un bruit impulsionnel dans une sequence video
EP1297494B1 (fr) Procede et dispositif de traitement de sequences d'images avec masquage
FR2829345A1 (fr) Dispositif et procede de codage d'images video par compensation de mouvement

Legal Events

Date Code Title Description
TP Transmission of property