FR2596221A1 - Systeme de filtre median adaptif - Google Patents
Systeme de filtre median adaptif Download PDFInfo
- 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
Links
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 33
- 230000004044 response Effects 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 15
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000013479 data entry Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000000034 method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/14—Picture signal circuitry for video frequency region
- H04N5/21—Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- 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/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/20024—Filtering details
- G06T2207/20032—Median 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)
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.
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)
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)
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 |
-
1986
- 1986-03-21 US US06/842,644 patent/US4682230A/en not_active Expired - Lifetime
-
1987
- 1987-03-19 DE DE3709068A patent/DE3709068C2/de not_active Expired - Lifetime
- 1987-03-20 JP JP62064479A patent/JPH0783236B2/ja not_active Expired - Fee Related
- 1987-03-20 FR FR8703916A patent/FR2596221B1/fr not_active Expired - Lifetime
- 1987-03-20 GB GB8706706A patent/GB2188511B/en not_active Expired
-
1992
- 1992-03-06 SG SG250/92A patent/SG25092G/en unknown
-
1994
- 1994-12-22 HK HK146894A patent/HK146894A/xx not_active IP Right Cessation
Non-Patent Citations (3)
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 |