FR2555003A1 - Procede et systeme d'amelioration du contraste d'une image - Google Patents

Procede et systeme d'amelioration du contraste d'une image Download PDF

Info

Publication number
FR2555003A1
FR2555003A1 FR8417206A FR8417206A FR2555003A1 FR 2555003 A1 FR2555003 A1 FR 2555003A1 FR 8417206 A FR8417206 A FR 8417206A FR 8417206 A FR8417206 A FR 8417206A FR 2555003 A1 FR2555003 A1 FR 2555003A1
Authority
FR
France
Prior art keywords
contrast
image
function
objects
modulation transfer
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
FR8417206A
Other languages
English (en)
Other versions
FR2555003B1 (fr
Inventor
Yair Shimoni
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.)
Elscint Ltd
Original Assignee
Elscint Ltd
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 Elscint Ltd filed Critical Elscint Ltd
Publication of FR2555003A1 publication Critical patent/FR2555003A1/fr
Application granted granted Critical
Publication of FR2555003B1 publication Critical patent/FR2555003B1/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/30Transforming light or analogous information into electric information
    • H04N5/32Transforming X-rays
    • H04N5/321Transforming X-rays with video transmission of fluoroscopic images
    • H04N5/325Image enhancement, e.g. by subtraction techniques using polyenergetic X-rays
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05GX-RAY TECHNIQUE
    • H05G1/00X-ray apparatus involving X-ray tubes; Circuits therefor
    • H05G1/08Electrical details
    • H05G1/60Circuit arrangements for obtaining a series of X-ray photographs or for X-ray cinematography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • G06T2207/10121Fluoroscopy
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Closed-Circuit Television Systems (AREA)
  • Controls And Circuits For Display Device (AREA)

Abstract

L'INVENTION CONCERNE L'IMAGERIE MEDICALE. UN SYSTEME DESTINE A CORRIGER LA PERTE DE CONTRASTE PRODUITE PAR DES EFFETS DE FONCTION DE TRANSFERT DE MODULATION, POUR AMELIORER LES IMAGES DE PETITS OBJETS SANS AFFECTER DEFAVORABLEMENT LES IMAGES DE PLUS GRANDS OBJETS, COMPREND NOTAMMENT DES MOYENS 13 DESTINES A ETABLIR DES FONDS HOMOGENES POUR LES OBJETS; DES MOYENS DESTINES A DETERMINER LE CONTRASTE DE L'IMAGE COMPTE TENU DES EFFETS DE FONCTION DE TRANSFERT DE MODULATION DU SYSTEME GLOBAL; ET DES MOYENS 17, 18 DESTINES A MANIPULER LE CONTRASTE DETERMINE POUR FAIRE EN SORTE QU'IL APPROCHE UN CONTRASTE IDEAL NON AFFECTE PAR LA FONCTION DE TRANSFERT DE MODULATION. APPLICATION A L'ANGIOGRAPHIE NUMERIQUE PAR SOUSTRACTION.

Description

2,5 5-4- 3
La présente invention concerne l'amélioration
d'images par ordinateur, et elle porte plus particulière-
mernt sur des procédés et des systèmes utilisant une mani-
pulation d'échelle dans la formation d'image par soustrac-
tion, pour corriger une dégradation du contraste. Lorsqu'un rayonnement traverse un petit trou ou une fente étroite, ou est atténué par une petite tige ou une barre étroite, l'image résultante "s'étale"-et est plus grande ou plus large que l'objet d'origine, et le contraste entre l'image et le fond est réduit. Ceci est le résultat combiné de la diffraction et de la diffusion de l'énergie, d'erreurs dans les moyens de détection, de bruit électronique et d'erreurs de numérisation. La forme fonctionnelle exacte de "l'étalement" est ce qu'on appelle
la "fonction de transfert de modulation" ou FTM. Les sys-
tèmes d'imagerie ont des FTM globales qui sont le résul-
tat de la combinaison des FTM de chacun des mécanismes d'étalement spécifiques qui interviennent dans l'image.En résumé, les images de trous et de fentes ou de tiges et de
barres sont "étalées" et ont un contraste réduit en compa-
raison des objets réels.
La FTM définit habituellement une augmentation
fixe des dimensions linéaires et une diminution corrélati-
ve du contraste d'un objet quelconque. Ainsi, lorsque la fente ou la barre ou le trou ou la tige est suffisamment grand, la taille et le contraste moyens sont en pratique inchangés par rapport à ceux de l'objet réel. Lorsque les objets sont petits, le contraste moyen dans l'image est proportionnel à la fois au contraste et à la largeur des
objets.
Il y a eu de nombreuses tentatives visant à cor-
riger l'étalement lié à la fonction de transfert de modula-
tion. A titre d'exemple, on a utilisé des filtres de "réta-
blissement" complexes qui tentent de redonner leur netteté
aux images étalées. Les filtres de rétablissement qu'on uti-
lise à l'heure actuelle correspondent à un traitement qui prend du temps et ils sont 'bruyants". Ainsi, ils renforcent le bruit déjà présent dans l'image et tendent ainsi à nuire à l'efficacité du système et à la qualité de l'image.
Dans de nombreuses études utilisant des tech-
niques d'amélioration d'image, le faible agrandissement de l'objet ne constitue pas un problème aussi important
que la dégradation du contraste qui est produite par l'éta-
lement dû à la FTM.
A titre d'exemple, en radioscopie numérique et dans des études similaires dans lesquelles on étudie le système vasculaire, le fait que les vaisseaux sanguins
soient légèrement agrandis a peu d'influence sur l'aptitu-
de à détecter ces vaisseaux. Au contraire, si le contraste
est réduit, il devient impossible de détecter ou de distin-
guer les plus petits vaisseaux sanguins.
Une solution au problème de l'étalement dû à
la FTM consiste à utiliser une "fenêtre" différente (étale-
ment de l'échelle de gris) pour observer les plus petits vaisseaux sanguins. Cette "solution" conduit à une échelle plus contrastée, qui fait perdre des détails essentiels dans de plus grands vaisseaux. De plus, lorsqu'on utilise cette "solution", on doit trouver manuellement la fenêtre
"correcte", ce qui est un processus qui prend du temps.
Une autre solution proposée au problème de l'étalement dû
à la FTM consiste à utiliser une échelle non linéaire.
C'est par exemple le cas pour la "double fenêtre", compre-
nant deux parties linéaires séparées qui s'appliquent res-
pectivement aux gros vaisseaux et aux petits vaisseaux. On a utilisé d'autres échelles non linéaires (par exemple des
échelles exponentielles et logarithmiques), mais en l'ab-
sence de connaissances théoriques sous-jacentes, elles ten-
dent à distordre le contraste relatif, à faire perdre des
détails et à rendre impossiblesdes calculs densitométriques.
Il existe donc un besoin réel et ressenti de-
puis longtemps concernant la correction de la dégradation
du contraste d'images sous l'effet de l'étalement des ob-
jets dans les images, à cause de la FTM.
L'invention présente donc des procédés préfé- rés destinés à être utilisés dans la formation d'image
par soustraction pour corriger la détérioration du con-
traste d'objetsqui est due à la FTM. Ces procédés amélio-
rent les images de petits objets sans affecter défavorable-
ment les images de grands objets, et la correction est ef-
fective même lorsqu'on n'utilise qu'une seule fenêtre pour
observer à la fois de grands objets et de petits objets.
Selon un aspect général de l'invention, ces procédés comprennent les opérations suivantes: - on établit des fonds homogènes pour les objets dont on désire former les images, - on détermine le contraste réel de l'image, tel qu'il est obtenu avec les effets de fonction de transfert de modulation, et - on manipule l'échelle du contraste réel déterminé de
l'image pour faire en sorte que le contraste réel dé-
terminé de l'image s'approche d'un contraste idéal non
affecté par la dégradation due à la fonction de trans-
fert de modulation.
Selon une caractéristique de l'invention, on
manipule l'échelle en utilisant une table préparée au préa-
lable, qui a été établie à partir de la manipulation d'une
fonction de contraste idéale.
On peut effectuer l'opération de manipulation de l'échelle de nombreuses manières entrant dans le cadre de l'invention. A titre d'exemple, l'une des manières consiste à manipuler l'échelle en utilisant une table préparée au
préalable, établie à partir d'une comparaison entre une fonc-
tion de contraste idéale et une fonction de constraste mesu-
rée. Les mesures exigées pour établir la table sont
habituellement fastidieuses et longues. Il est donc sou-
haitable d'utiliser des procédés qui exigent la quantité minimale de mesures. Pour réduire la quantité de mesures, on peut déterminer la pente de la fonction à partir de con- sidérations théoriques, si on fait certaines hypothèses, comme par exemple celle consistant à considérer que l'image ne contient que des vaisseaux dans le plan étudié et que la densité de la matière contrastante est constante. Si on
fait des hypothèses appropriées, les seules mesures néces-
saires pour obtenir la forme de la fonction de réponse sont des mesures de FTM qu'on utilise également pour la normalisation et le cadrage. On peut aisément obtenir la pente de la valeur mesurée en fonction de la valeur vraie
à partir des mesures de FTM.
Un autre procédé de manipulation comprend les opérations suivantes:
- en utilisant une FTM approchée avec une forme échantil-
lon, qui est une fonction triangulaire ayant la même largeur moyenne que la FTM exacte connue, on détermine approximativement la perte de contraste qui est due à
l'étalement résultant de la FTM pour des tailles d'ob-
jet données. La fonction triangulaire donne les rela-
tions suivantes pour le contraste en fonction de la taille: C = A(d2-d 3/3) lorsque d. 1 C = A(d -d2/3) lorsque d > 1 tandis que la fonction de contraste idéale est: C = A*d lorsque d > 1 2C = A*d3 lorsque d 1 On a ici: A =per d = R/r avec les notations suivantes: est une constante de proportionalité; est la densité de la matière de contraste; r est la largeur moyenne (largeur à la moitié de la hauteur maximale) de la FTM; et R est la dimension d'un côté du vaisseau sanguin
(supposé être un cylindre de section carrée).
L'invention sera mieux comprise à la lecture de
la description qui va suivre de modes de réalisation et en
se référant aux dessins annexés sur lesquels: La figure 1 est un schéma synoptique montrant un
système d'imagerie destiné à utiliser les procédés de cor-
rection de dégradation du contraste conformes à l'invention; La figure 2a est une représentation graphique d'une FTM caractéristique; Les figures 2b et 2c sont des représentations graphiques de différentes versions de courbes de réponse basées sur la FTM; La figure 3 est une représentation graphique montrant les résultats de l'utilisation d'un procédé de correction de la dégradation du contraste; et
La figure 4 est une autre représentation graphi-
que montrant les résultats de l'utilisation d'un autre pro-
cédé de correction de la dégradation du contraste.
Comme le montre la figure 1, le procédé convient
parfaitement à des systèmes d'imagerie du type à soustrac-
tion comme ceux qu'on utilise en radioscopie numérique. Le
système 11 de la figure 1 est un système de radioscopie nu-
mérique caractéristique. Une unité d'acquisition d'image 12 effectue l'acquisition d'images qui sont prises avant et
après l'ajout de la matière de contraste. L'unité de sous-
traction 13 soustrait les images acquises avant l'ajout de la matière de contraste, par rapport aux images acquises après l'ajout de la matière de contraste, pour donner des
images d'objets ayant des fonds homogènes. Les données d'ima-
6 2555003
ge qui sont émises par l'unité de soustraction 13 sont de
préférence, mais pas nécessairement, sous forme numérique.
On notera que les images considérées ici sont fondamenta-
lement des images de petits objets, par exemple des vais-
seaux sanguins, ayant un effet de FTM qui réduit le con-
traste dans l'image finale.
Il existe des moyens prévus pour corriger la ré-
duction du contraste qui se produit normalement lorsqu'on forme l'image de vaisseaux sanguins. Plus précisément, il existe des moyens destinés à manipuler l'échelle de gris de l'image sur le moniteur 14 pour améliorer le contraste
en rétablissant le contraste perdu à cause de la caractéris-
tique de fonction de transfert de modulation. L'image appa-
raît sur le moniteur lorsque le convertisseur numérique-
analogique 16 convertit en signaux vidéo analogiques les données qui proviennent de l'unité de soustraction 13. Le système est sous la commande d'un dispositif de traitement
17 qui entre autres, manipule l'échelle de gris pour défi-
nir le contraste. Comme le montre la figure 1, le système utilise de préférence dans ce but une unité arithmétique
18 et une table 19, sous la commande du dispositif de trai-
tement. Les figures 3 et 4 montrent différents types de
manipulations d'échelle de gris qui sont envisagées. Le ty-
pe effectivement utilisé dépend de l'équipement disponible et des résultats désirés. A titre d'exemple, lorsque le but
est d'obtenir des images montrant de petits vaisseaux san-
guins, on utilise une manipulation d'échelle de gris diffé-
rente de celle qu'on utilise lorsque le but est de former
l'image de gros vaisseaux sanguins. Dans toutes les manipu-
lations de l'échelle de gris, il est important que l'objet soit observé sur un fond relativement homogène. Toutes les
manipulations décrites ci-après améliorent l'image de vais-
seaux de la taille particulière désirée, sans affecter défa-
vorablement le contraste et donc la visibilité des vaisseaux
ayant des tailles différentes.
Pour expliquer les principes sur lesquels repo-
se l'invention, la figure 2a représente graphiquement
l'étalement d'une source ponctuelle qui est dû à la FTM.
Les axes représentent l'intensité en fonction de la distan-
ce à partir de la position de l'image linéaire de la sour-
ce ponctuelle. L'image idéale est concentrée en un point (p) ou en une ligne L sans dimensions. L'image réelle est
étalée par la FTM. La courbe représentée est fondamenta-
lement une fonction symétrique à un seul maximum, par exem-
ple une fonction gaussienne. On peut utiliser dans le ca-
dre de l'invention de nombreux types de fonctions, en plus
de fonctions gaussiennes, comme des fonctions triangulai-
res, rectangulaires, à échelons et de type mesa.
La figure 2b montre la fonction de réponse du contraste, c'est-à-dire le contraste mesuré divisé par le contraste vrai en fonction de la dimension linéaire R du
côté d'un vaisseau sanguin, avec l'hypothèse selon laquel-
le le vaisseau sanguin est un tube creux ayant une section transversale rectangulaire et avec un côté parallèle à la direction d'observation. Lorsque la dimensicn du côté du vaisseau sanguin augmente, la valeur mesurée s'approche de la valeur vraie. Pour de petits vaisseaux, la fonction de réponse tend asymptotiquement vers une fonction linéaire
du côté.
La figure 2c montre le même contraste mesuré (réel) en fonction du contraste vrai. Cette représentation permet de déterminer plus aisément les corrections qu'il est nécessaire d'apporter au contraste à n'importe quel point
désiré.
Dans une procédure préférée de l'invention, on
suppose que la fonction de transfert de modulation est trian-
gulaire. Il entre cependant dans le cadre de l'invention de considérer pour cette fonction une forme gaussienne ou toute
autre forme.
Si on suppose que les vaisseaux sanguins sont
des tubes creux carrés, avec leurs côtés toujours paral-
lèles à la direction d'observation, plutôt que des tubes
cylindriques, on peut montrer que le contraste vrai en an-
giographie numérique par soustraction est: CCC OR, en dé- signant par 4 la densité de la matière de contraste, par R les dimensions du côté du tuyau et par une constante
de proportionnalité. La FTM fait que la fonction de répon-
se est approximativement: R Contraste mesuré/contraste vrai = R-r R+I r en désignant par r la largeur moyenne de la FTM, et par Y
une constante de proportionnalité qui dépend dans une cer-
taine mesure de la forme de la FTM.
Dans de petits vaisseaux sanguins, le rapport entre le contraste mesuré et contraste vrai est donc égal
à r (comme le montre la figure 2b)2et le contraste me-
suré est approximativement égal à é ir. Ceci indique que l'image de petits vaisseaux sanguins présente des niveaux d'énergie réduits. En utilisant la notation: A =Pr d= R/rr on a dans le cas idéal, sans effets de FTM: C = A*d pour tout R, mais en réalité, à cause des effets de FTM et avec une FTM
triangulaire, ce sont les équations suivantes qui s'appli-
quent: C = A (d-d /3) lorsque d >m 1 C = A (d2/2 - d3/3) lorsque d. 1 La figure 2c montre la fonction idéale et la
fonction réelle. Si on connaît le terme r de la FTM, l'uni-
té arithmétique peut calculer la valeur théorique idéale
que doit avoir le contraste pour n'importe quelle valeur me-
surée. Cependant, la correction de chaque valeur d'élément
d'image dans une grande image peut prendre trop de temps.
Un procédé préféré consiste donc à calculer une table à con-
sulter dont l'information d'entrée est une valeur mesurée quel-
conque et dont l'information de sortie correspondante est
la valeur idéale approchée. De tels mécanismes de consul-
tation de table sont bien connus dans les techniques de
traitement d'images.
Si on ne connaît pas le paramètre r, on peut supposer que sa valeur est égale à r et l'utiliser dans a le procédé et les moyens de correction. L'utilisation de ra au lieu de r change les valeurs de contraste obtenues qui deviennent: (r-r a) C' =A.d 1(-11/3 (R-r /3) j lorsque d 1 et a r 2 (r-r) C' =A*d (1 -) lorsque d C 1 r (r -R/3)
Ainsi, lorsque ra r, le résultat obtenu est proche du ré-
sultat idéal. Parmi les approximations qui sont utiles pour
améliorer l'image de petits vaisseaux figure celle qui uti-
lise la relation: C" = A.Cmesuré Le résultat obtenu par cette approximation peut s'écrire au moyen de variables A et d de la façon suivante: C" = A \-d (l-d/6) lorsque d > 1 et C" = A.d (l-d/6) lorsque d C 1 ce qui est une bonne approximation pour de petits vaisseaux
(d 41).
Cette relation est représentée sur la figure 3.
Une autre approximation utile pour de petits vaisseaux sur une certaine "fenêtre"'seulement de la plage complète de valeurs correspond à ce qu'on peut appeler
"l'1élévation carrée de l'inverse".
On peut choisir une "fenêtre intéressante" et transformer la courbe contenue dans la fenêtre pour qu'elle soit proche de la courbe linéaire idéale. On supposera que
la fenêtre se trouve entre la valeur de seuil T et la va-
2555003
leur supérieure U. On obtient une échelle "étirée" par la
transformation faisant passer T à O et U à une limite nu-
mérique supérieure W, de façon que:
C' = W (C-T)/(U-T)
Une échelle inversée est définie par
C" = W -C' - W(U-C)/(U-T)
dans laquelle T est remplacé par W et U est remplacé par O.
L'opération "d'inversion" est sa propre opération complé-
mentaire:
C"' = W - C" = C'
Une fenêtre spéciale commence à T = O. On a alors: C' = W/U -C (C correspond à la plage O-U et C' à la plage O-W) L'échelle "inversée" est: C" = W(U-C)/U (toujours dans les mêmes plages) A ce point, on utilise une échelle quadratique: C"' =]/Wc"2 =W(U-C)2/U2 (dans les mêmes plages) En inversant maintenant à nouveau, on obtient: CiV = W-W(U-C)/U (dans les mêmes plages)
et dans le cas spécial dans lequel il n'y a pas d'étire-
ment: W-U Cv = U-(U-C)2/U = 2C-C2/U (C et Cv sont dans la plage O-U) Si le contraste mesuré est C, on peut décrire le contraste de la manière suivante, en utilisant A et d: cVy 2Ad2 _ 2Ad3/3 lorsque d <1 et Cv 2 2A CV 2Ad - (2A/3-A /U) d lorsque d >1 ce qui correspond à une fonction qui démarre avec une loi quadratique (parabole) pour dû 0, puis franchit la valeur
idéale à environ d*0,64 et dû1l,5 (voir la figure 4).
Cette approximation est donc valide pour une "fenêtrd'plus large que la précédente. Cependant, on ne doit pas utiliser cette échelle pour C >U, du fait qu'elle donne une solution à deux valeurs. Elle est très utile pour
des images réelles lorsqu'on prend soin de choisir U supé-
rieur à toute valeur de C dans l'image. Elle est particu-
lièrement utile du fait que des images obtenues par ana-
lyse numérique de signal sont en réalité souvent formées avec une échelle inversée, de façon que la matière de contraste apparaisse en noir. I1 va de soi que de nombreuses modifications peuvent être apportées aux dispositifs et aux procédés particuliers décrits et représentés, sans sortir du cadre
de l'invention. On notera en particulier qu'il est possi-
ble de faire d'autres approximations en demeurant dans le
cadre de l'invention.

Claims (18)

REVENDICATIONS -
1. Procédé de correction de la perte de con-
traste produite par des effets de fonction de transfert de modulation dans la formation d'une image, ce procédé améliorant les images de petits objets sans affecter dé-
favorablement les images de plus grands objets, et la cor-
rection étant effective même lorsqu'on n'utilise qu'une seule fenêtre pour observer à la fois des objets grands et petits, caractérisé en ce que: on établit des fonds homogènes pour des objets dont on désire former l'image; on détermine le contraste de l'image qui est obtenu avec
des effets de fonction de transfert de modulation du sys-
tème global; et on manipule l'échelle du contraste déter-
miné de l'image pour faire en sorte que le contraste dé-
terminé s'approche d'un contraste idéal non affecté par la réduction du contraste due à la fonction de transfert
de modulation.
2. Procédé selon la revendication 1, caractéri-
sé en ce que l'opération de manipulation de l'échelle com-
prend l'opération qui consiste à utiliser une table à con-
sulter (19) préparée au préalable.
3. Procédé selon la revendication 2, caractéri-
sé en ce que l'opération de manipulation comprend la pré-
paration de la table à consulter (19), par la comparaison d'une fonction de contraste idéale avec une fonction de
contraste mesurée.
4. Procédé selon la revendication 1, caractéri-
sé en ce que l'opération de manipulation comprend l'utili-
sation d'une largeur moyenne connue pour la fonction de
transfert de modulation.
5. Procédé selon la revendication 1, caractéri-
sé en ce que l'opération de manipulation comprend l'utilisa-
tion d'une largeur moyenne approchée pour la fonction de
transfert de modulation.
6. Procédé selon l'une quelconque des revendica-
tions 1 à 5, caractérisé en ce qu'il comprend l'utilisa-
tion d'une fonction de transfert de modulation triangu-
laire prise par hypothèse.
7. Procédé selon l'une quelconque des revendi-
cations 1 à 6, caractérisé en ce qu'il comprend l'utilisa- tion de l'hypothèse selon laquelle des objets longs et minces dans l'image tels que des veines ont des sections
transversales carrées.
8. Procédé selon la revendication 1, caracté-
risé en ce que l'opération consistant à établir des fonds
homogènes comprend l'utilisation de procédés angiographi-
ques soustractifs numériques.
9. Procédé selon la revendication 1, caracté-
risé en ce que l'opération de manipulation comprend l'élé-
vation au carré de l'inverse.
10. Système destiné à corriger la perte de con-
traste produite par des effets de fonction de transfert de modulation dans la formation d'une image, ce système
améliorant les images de petits objets sans affecter défa-
vorablement les images de plus grands objets, et la correc-
tion étant effective même lorsqu'on n'utilise qu'une seule fenêtre pour observer à la fois des objets grands et petits, caractérisé en ce qu'il comprend: des moyens (13) destinés
à établir des fonds homogènes pour des objets dont on dé-
sire former l'image; des moyens destinés à déterminer le contraste de l'image qui est obtenu avec des effets de fonction de transfert de modulation du système global; et
des moyens (17, 18) destinés à manipuler l'échelle du con-
traste déterminé de l'image pour faire en sorte que le contraste déterminé s'approche d'un contraste idéal non affecté par la réduction de contraste due à la fonction de
transfert de modulation.
11. Système selon la revendication 10, caracté-
risé en ce que les moyens de manipulation (17, 18) compren-
nent des moyens destinés à utiliser une table à consulter
(19) qui a été préparée.
12. Système selon la revendication 11, carac-
térisé en ce que les moyens de manipulation (17, 18) com-
prennent des moyens destinés à préparer la table à consul-
ter (19), et les moyens de préparation de la table à con- sulter comprennent des moyens destinés à comparer une
fonction de contraste idéale avec une fonction de contras-
te mesurée.
13. Système selon la revendication 10, caracté-
risé en ce que les moyens de manipulation (17, 18) com-
prennent des moyens destinés à déterminer et à utiliser une largeur moyenne connue pour la fonction de transfert
de modulation.
14. Système selon la revendication 10, caracté-
risé en ce que les moyens de manipulation (17, 18) com-
prennent des moyens destinés à utiliser une largeur moyen-
ne approchée pour la fonction de transfert de modulation.
15. Système selon l'une quelconque des revendi-
cations 10 à 14, caractérisé en ce qu'il comprend des mo-
yens (17, 18) destinés à utiliser une fonction de trans-
fert de modulation triangulaire qui résulte d'une hypothè-
se.
16. Système selon l'une quelconque des revendi-
cations 10 à 15, caractérisé en ce que des objets longs et minces dans l'image, tels que des veines, sont supposés
avoir des sections transversales rectangulaires.
17. Système selon la revendication 10, caractéri-
sé en ce que les moyens destinés à établir des fonds homo-
gènes comprennent des moyens angiographiques soustractifs
numériques (13).
18. Système selon la revendication 10, caractéri-
sé en ce que les moyens de manipulation (17, 18) comprennent
des moyens destinés à élever au carré l'inverse.
FR848417206A 1983-11-13 1984-11-12 Procede et systeme d'amelioration du contraste d'une image Expired - Lifetime FR2555003B1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
IL70214A IL70214A (en) 1983-11-13 1983-11-13 Image contrast enhancement arrangement

Publications (2)

Publication Number Publication Date
FR2555003A1 true FR2555003A1 (fr) 1985-05-17
FR2555003B1 FR2555003B1 (fr) 1991-11-29

Family

ID=11054673

Family Applications (1)

Application Number Title Priority Date Filing Date
FR848417206A Expired - Lifetime FR2555003B1 (fr) 1983-11-13 1984-11-12 Procede et systeme d'amelioration du contraste d'une image

Country Status (6)

Country Link
US (1) US4907288A (fr)
JP (1) JP2575612B2 (fr)
DE (1) DE3441162C2 (fr)
FR (1) FR2555003B1 (fr)
IL (1) IL70214A (fr)
NL (1) NL193182C (fr)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0235042A1 (fr) * 1986-02-28 1987-09-02 General Electric Cgr S.A. Système d'imagerie numérique à réglage de l'échelle des gris, notamment pour la visualisation des vaisseaux sanguins
WO1989012278A1 (fr) * 1988-06-08 1989-12-14 Quantel Limited Procede et appareil d'amelioration d'images
EP0349075A1 (fr) * 1988-07-01 1990-01-03 Laboratoires D'electronique Philips Méthode numérique de modification d'image et système de mise en oeuvre
FR2648302A1 (fr) * 1989-06-09 1990-12-14 Peugeot Procede de traitement d'un signal composite video d'une scene a examiner
EP0418921A2 (fr) * 1989-09-22 1991-03-27 Fuji Photo Film Co., Ltd. Méthode et appareil de traitement de soustraction d'énergie, méthode et appareil de traitement de superposition et appareil de lecture d'image de radiation
FR2666426A1 (fr) * 1990-08-31 1992-03-06 Gen Electric Cgr Procede de correction des mesures de densite optique effectuees sur un film radiographique.
FR2670979A1 (fr) * 1990-12-21 1992-06-26 Philips Electronique Lab Procede de segmentation binaire locale d'images numerisees, par seuillage d'histogramme.
EP0756246A1 (fr) * 1995-07-25 1997-01-29 Kabushiki Kaisha Toyota Chuo Kenkyusho Appareil pour le traitement d'images

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01248187A (ja) * 1988-03-30 1989-10-03 Toshiba Corp ディスプレイシステム
NL8900217A (nl) * 1989-01-30 1990-08-16 Philips Nv Werkwijze voor het bepalen van een modulatie overdrachtfunctie van een digitaal beeldvormend systeem.
CA2028092C (fr) * 1989-10-23 1996-08-06 Yuji Tanaka Appareil de traitement d'images binaires
NL9001265A (nl) * 1990-06-05 1992-01-02 Philips Nv Werkwijze voor bepaling van een modulatie-overdrachtfunctie van een digitaal beeldvormend systeem.
US5125046A (en) * 1990-07-26 1992-06-23 Ronald Siwoff Digitally enhanced imager for the visually impaired
US5359675A (en) * 1990-07-26 1994-10-25 Ronald Siwoff Video spectacles
US5271064A (en) * 1991-06-14 1993-12-14 University Of Cincinnati Apparatus and method for smoothing regions and enhancing edges in gray scale images
US6466689B1 (en) * 1991-11-22 2002-10-15 Arch Development Corp. Method and system for digital radiography
US5297550A (en) * 1992-08-06 1994-03-29 Picker International, Inc. Background darkening of magnetic resonance angiographic images
US5297551A (en) * 1992-08-06 1994-03-29 Picker International, Inc. Weighted ray projection imaging for MR angiography
GB9420654D0 (en) * 1994-10-13 1994-11-30 Kodak Ltd A method for improving the contrast in photographic film materials
JPH09261490A (ja) * 1996-03-22 1997-10-03 Minolta Co Ltd 画像形成装置
US6864888B1 (en) * 1999-02-25 2005-03-08 Lockheed Martin Corporation Variable acuity rendering for a graphic image processing system
GB9904692D0 (en) * 1999-03-01 1999-04-21 Isis Innovation X-ray image processing
GB0028491D0 (en) * 2000-11-22 2001-01-10 Isis Innovation Detection of features in images
JP2002271611A (ja) * 2001-03-14 2002-09-20 Fujitsu Ltd 画像処理装置
US20030160800A1 (en) 2002-02-22 2003-08-28 Agfa-Gevaert Multiscale gradation processing method
TWI237497B (en) * 2004-05-07 2005-08-01 Quanta Comp Inc High speed and high accuracy image adjust apparatus and method
US20060056722A1 (en) * 2004-09-14 2006-03-16 Nathan Moroney Edge preserving method and apparatus for image processing
DE102004049227B4 (de) * 2004-10-08 2007-03-01 Yxlon International Security Gmbh Verfahren zur Bestimmung der Lageänderung eines Objekts in einem Gepäckstück
JP2009211313A (ja) * 2008-03-03 2009-09-17 Fujitsu Ltd 画像処理装置、画像処理方法および画像処理プログラム

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0101336A2 (fr) * 1982-06-21 1984-02-22 Thomson-Csf Appareil de génération d'un signal vidéo traité, à contraste amélioré pour radiologie, et procédé pour traiter un tel signal

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4193089A (en) * 1976-03-05 1980-03-11 Picker Corporation Television radiation imaging system and method
JPS5914790B2 (ja) * 1976-08-16 1984-04-06 富士電機株式会社 光学的識別装置
US4156914A (en) * 1977-08-18 1979-05-29 Baird Corporation Computer image display system and processor therefor
US4317179A (en) * 1978-12-26 1982-02-23 Fuji Photo Film Co., Ltd. Method and apparatus for processing a radiographic image
US4315318A (en) * 1978-12-26 1982-02-09 Fuji Photo Film Co., Ltd. Method and apparatus for processing a radiation image
JPS5691735A (en) * 1979-12-25 1981-07-24 Fuji Photo Film Co Ltd Method and apparatus for treating xxray image
SE426281B (sv) * 1981-05-05 1982-12-20 Ericsson Telefon Ab L M Forfarande for visning pa en bildskerm av utsignaler fran en stralningskenslig detektor
US4394688A (en) * 1981-08-25 1983-07-19 Hamamatsu Systems, Inc. Video system having an adjustable digital gamma correction for contrast enhancement
JPS5834677A (ja) * 1981-08-25 1983-03-01 Matsushita Graphic Commun Syst Inc 画情報読取装置のmtf補正方式
US4517607A (en) * 1981-11-09 1985-05-14 Ricoh Company, Ltd. Method of and apparatus for compensating image in image reproduction system
US4438495A (en) * 1981-11-13 1984-03-20 General Electric Company Tomography window-level gamma functions
JPS5928144A (ja) * 1982-08-09 1984-02-14 Fuji Photo Film Co Ltd 放射線画像再生装置
US4559557A (en) * 1984-06-01 1985-12-17 General Electric Company Region-of-interest digital subtraction angiography

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0101336A2 (fr) * 1982-06-21 1984-02-22 Thomson-Csf Appareil de génération d'un signal vidéo traité, à contraste amélioré pour radiologie, et procédé pour traiter un tel signal

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, tome NS-29, no. 3, juin 1982, pages 1176-1180, New York, US; W.R. BRODY "Digital subtraction angiography" *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0235042A1 (fr) * 1986-02-28 1987-09-02 General Electric Cgr S.A. Système d'imagerie numérique à réglage de l'échelle des gris, notamment pour la visualisation des vaisseaux sanguins
FR2595153A1 (fr) * 1986-02-28 1987-09-04 Thomson Cgr Systeme d'imagerie numerique a reglage de l'echelle des gris, notamment pour la visualisation des vaisseaux sanguins
US4827492A (en) * 1986-02-28 1989-05-02 Thomson-Cgr Digital imaging system with grey scale setting, especially for the display of blood vessels
WO1989012278A1 (fr) * 1988-06-08 1989-12-14 Quantel Limited Procede et appareil d'amelioration d'images
EP0349075A1 (fr) * 1988-07-01 1990-01-03 Laboratoires D'electronique Philips Méthode numérique de modification d'image et système de mise en oeuvre
FR2633746A1 (fr) * 1988-07-01 1990-01-05 Labo Electronique Physique Methode numerique de modification d'image et systeme de mise en oeuvre
US4999711A (en) * 1988-07-01 1991-03-12 U.S. Philips Corp. Digital method of modifying an image contrast so as to retain imaging of small objects
FR2648302A1 (fr) * 1989-06-09 1990-12-14 Peugeot Procede de traitement d'un signal composite video d'une scene a examiner
EP0418921A2 (fr) * 1989-09-22 1991-03-27 Fuji Photo Film Co., Ltd. Méthode et appareil de traitement de soustraction d'énergie, méthode et appareil de traitement de superposition et appareil de lecture d'image de radiation
EP0418921A3 (en) * 1989-09-22 1991-12-18 Fuji Photo Film Co., Ltd. Energy subtraction processing method and apparatus, superposition processing method and apparatus, and radiation image read-out apparatus
US5315507A (en) * 1989-09-22 1994-05-24 Fuji Photo Film Co., Ltd. Energy subtraction processing method and apparatus, superposition processing method and apparatus, and radiation image read-out apparatus
FR2666426A1 (fr) * 1990-08-31 1992-03-06 Gen Electric Cgr Procede de correction des mesures de densite optique effectuees sur un film radiographique.
EP0478406A1 (fr) * 1990-08-31 1992-04-01 General Electric Cgr S.A. Procédé de correction des mesures de densité optique effectuées sur un film radiographique
US5297036A (en) * 1990-08-31 1994-03-22 General Electric Cgr S.A. Method for the correction of the measurements of optical density made on a radiographic film
FR2670979A1 (fr) * 1990-12-21 1992-06-26 Philips Electronique Lab Procede de segmentation binaire locale d'images numerisees, par seuillage d'histogramme.
EP0493855A1 (fr) * 1990-12-21 1992-07-08 Laboratoires D'electronique Philips S.A.S. Procédé de segmentation binaire locale d'images numérisées, par seuillage d'histogramme
US5832111A (en) * 1990-12-21 1998-11-03 U.S. Philips Corporation Method and device for binary segmentation of a signal representing a digital image with a threshold obtained using a histogram
EP0756246A1 (fr) * 1995-07-25 1997-01-29 Kabushiki Kaisha Toyota Chuo Kenkyusho Appareil pour le traitement d'images
US6160576A (en) * 1995-07-25 2000-12-12 Kabushiki Kaisha Toyota Chuo Kenkyusho Image compensating device based on age-related characteristics

Also Published As

Publication number Publication date
US4907288A (en) 1990-03-06
FR2555003B1 (fr) 1991-11-29
DE3441162C2 (de) 1994-01-13
NL193182C (nl) 1999-02-02
IL70214A (en) 1987-10-20
JP2575612B2 (ja) 1997-01-29
JPS60169980A (ja) 1985-09-03
NL193182B (nl) 1998-10-01
IL70214A0 (en) 1984-02-29
DE3441162A1 (de) 1985-07-25
NL8403468A (nl) 1985-06-03

Similar Documents

Publication Publication Date Title
FR2555003A1 (fr) Procede et systeme d&#39;amelioration du contraste d&#39;une image
Soler et al. Fourier depth of field
Battiato et al. High dynamic range imaging for digital still camera: an overview
US20180067197A1 (en) Method and system for real-time motion artifact handling and noise removal for tof sensor images
Huck et al. Visual communication: an information theory approach
WO2005122086A2 (fr) Reduction de bruit d’image numerique
EP0840252A1 (fr) Procédé de traitement d&#39;image numérique pour l&#39;extraction automatique d&#39;objets en forme de rubans
EP0114559B1 (fr) Procédé de correction permanente de couleur applicable à la prise de vues vidéo et système mettant en oeuvre ce procédé
FR2885717A1 (fr) Procede d&#39;estimation de correction d&#39;une image et dispositif d&#39;angiographie le mettant en oeuvre
FR2829268A1 (fr) Procede de traitement d&#39;images pour angiographie soustractiv e numerisee
EP1298592A1 (fr) Stabilisation des images d&#39;une scène, correction des offsets de niveaux de gris, détection d&#39;objets mobiles et harmonisation de deux appareils de prise de vues fondés sur la stabilisation des images
CN110930311B (zh) 一种提高红外图像与可见光图像融合的信噪比方法与装置
Fonseca et al. Combined interpolation—restoration of Landsat images through FIR filter design techniques
FR2831753A1 (fr) Procede d&#39;adaptation de contraste de plusieurs images sur une image de reference commune
FR2633746A1 (fr) Methode numerique de modification d&#39;image et systeme de mise en oeuvre
Li et al. Image deblurring for satellite imagery using small-support-regularized deconvolution
CN109410180B (zh) 衰减系数的确定方法、装置、计算机设备和存储介质
CN112203068A (zh) 单像素成像方法、系统、装置及介质
EP0782327A1 (fr) Procédé de modification de la résolution d&#39;une image numérisée
FR2547966A1 (fr) Circuit numeriseur adaptatif pour un systeme de traitement d&#39;information
Du Research on image de‐disturbing algorithm based on dark channel prior and anisotropic Gaussian filtering
Anusha Automated Speckle Noise Suppression from Sentinel-1A Synthetic Aperture Radar Imagery Using Adaptive Filtering Techniques
Xu et al. Underwater Fourier single pixel imaging based on water degradation function compensation method
Liu et al. Low-light Image Enhancement Using Chain-consistent Adversarial Networks
Chopra et al. Mathematical formulation for optimization of the application of image processing technique for pollution control