IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

MATLAB Discussion :

Ameliorer la nettete d'une image avec un filtre


Sujet :

MATLAB

  1. #1
    Membre du Club
    Homme Profil pro
    Étudiant
    Inscrit en
    Mars 2018
    Messages
    171
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Maine et Loire (Pays de la Loire)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Mars 2018
    Messages : 171
    Points : 55
    Points
    55
    Par défaut Ameliorer la nettete d'une image avec un filtre
    Bonjour,

    J'aimerais creer un script matlab pour ameliorer la nettete d'une image en me basant sur les travaux de cet article: https://ieeexplore.ieee.org/document/6271882

    L'auteur explique comment utiliser un filtre a derivee fractionnelle combine a l'effet de bandes de mach pour ameliorer la nettete de l'image.

    Le filtre est le suivant:
    Nom : fonction.jpg
Affichages : 480
Taille : 31,2 Ko
    Avec a(k) egal a:
    Nom : a(k).jpg
Affichages : 571
Taille : 61,6 Ko

    Afind'essayer d'appliquer ce filtre a une image, j'aimerais coder la formule presentee ci-dessus sur matlab avec le code suivant:

    Code MATLAB : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    z=0:1:15;
    L=4;
    v=0.3;
    e=filtre(z,L,v);
    figure
    plot(z,e);
     
    function [A] = filtre(z,L,v)
        k = 0:1:L;
     
        S = sum(((-1).^k)*nchoosek(v,k)*(z.^(-k)));
    end

    Cependant, lorsque j'essaye de compiler je recois l'erreur suivante:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    >> script
    Error using nchoosek (line 29)
    The second input has to be a non-negative integer.
     
    Error in script>filtre (line 29)
        S = sum(((-1).^k)*nchoosek(v,k)*(z.^(-k)));
     
    Error in script (line 22)
    e=filtre(z,L,v);
    Je ne comprends pas d'ou vient cette erreur etant donne que mon parametre k n'est jamais negatif.
    Si vous avez des avis sur mon code et des idees sur comment corriger l'erreur obtenue je serais ravi de lire vos reponses.
    Merci!

  2. #2
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    k doit être un entier, pas un vecteur d'entiers.
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  3. #3
    Membre du Club
    Homme Profil pro
    Étudiant
    Inscrit en
    Mars 2018
    Messages
    171
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Maine et Loire (Pays de la Loire)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Mars 2018
    Messages : 171
    Points : 55
    Points
    55
    Par défaut
    Bonjour,

    Merci pour ta reponse.

    J'ai enleve le vecteur k et a la place j'ai place ma somme dans une boucle for:
    Code MATLAB : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    function [A] = filtre(z,L,v)
     
        for k =0:L,
         S = sum(((-1).^k).*nchoosek(v,k).*(z.^(-k)));
        end
    end

    J'ai desormais l'erreur suivante:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    >> script
    Error using  .* 
    Matrix dimensions must agree.
     
    Error in script>filtre (line 29)
         S = sum(((-1).^k).*nchoosek(v,k).*(z.^(-k)));
     
    Error in script (line 22)
    e=filtre(z,L,v);
    J'ai essaye d'isoler chaque element de la somme et l'erreur se trouve dans cette partie du code:
    Code matlab : Sélectionner tout - Visualiser dans une fenêtre à part
    nchoosek(v,k).*(z.^(-k))

    Apparement je ne peux pas multiplier ces deux elements car les dimensions des deux matrices ne correspondent pas, mais je ne vois pas vraiment en quoi ce sont des matrices ?

  4. #4
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Le premier argument de nchoosek devrait être un entier, non ?

    Car :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    >> nchoosek(0.3,0)
     
    ans =
     
       Empty matrix: 1-by-0
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  5. #5
    Membre du Club
    Homme Profil pro
    Étudiant
    Inscrit en
    Mars 2018
    Messages
    171
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Maine et Loire (Pays de la Loire)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Mars 2018
    Messages : 171
    Points : 55
    Points
    55
    Par défaut
    Justement l'objectif de l'utilisation de la derivee fractionnaire est de pouvoir obtenir des filtres d'ordre v avec v un reel (et non un entier).

    En regardant la fonction nchoosek, il faut effectivement mettre des entiers non negatifs en parametre. Sais-tu s'il existe un equivalent de nchoosek avec des parametres pouvant etre reels?
    Ou sinon comment coder une telle fonction?

    J'ai l'impression qu'il y a quelques pistes par ici mais la reponse semble dire qu'il n'est pas possible de calculer de coeficient binominaux ayant des parametres autre que entier: https://fr.mathworks.com/matlabcentr...nd-k-arguments

Discussions similaires

  1. Recuperation de la taille d'une image avec la propriété SRC
    Par CDRIK dans le forum Général JavaScript
    Réponses: 6
    Dernier message: 08/10/2004, 04h05
  2. [Image]Créer une image avec JAVA 1.1
    Par burno dans le forum 2D
    Réponses: 4
    Dernier message: 11/08/2004, 10h19
  3. Afficher une image avec swing
    Par arnonote dans le forum AWT/Swing
    Réponses: 11
    Dernier message: 07/08/2004, 13h10
  4. [Débutante] Création d'une image avec un composant
    Par gwendo dans le forum AWT/Swing
    Réponses: 9
    Dernier message: 09/07/2004, 10h58
  5. PL/SQL - Comment afficher une image avec HTP ?
    Par patmaba dans le forum PL/SQL
    Réponses: 2
    Dernier message: 08/07/2004, 10h28

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo