Différences

Ci-dessous, les différences entre deux révisions de la page.

--- wiki:projet:trait_image [2017/03/29 17:30]
http_fablab.sorbonne-universites.fr_wiki_doku.php_id_wiki_projet_laser
+++ wiki:projet:trait_image [2017/04/24 09:54] (Version actuelle)
http_fablab.sorbonne-universites.fr_wiki_doku.php_id_wiki_projet_laser [COMMENTAIRES]
@@ Ligne 7: / Ligne 7: @@
 <code>
+close all
 %-------------------------------------------------------------------------%
 %                   Traitement de l'image                                 %
@@ Ligne 12: / Ligne 13: @@
 %   Ce premier fichier est le pogramme de traitement. Il lit une image,   %
 %   puis appelle differentes fonctions afin de la traiter.                %
+clear figure
-Im = imread('nom_fichier.jpg'); %Im : variable contenant 3 matrices qui
+Im = imread('faisc_31_03_a.jpg'); %Im : variable contenant 3 matrices qui
                                 %contiennent les intensites de rouge, vert
                                 %pour chaque pixel de l'image de base
@@ Ligne 39: / Ligne 40: @@
 image(imf2) %affichage de l'image imf2
-[imech, ech] = echelle(Im); %appelle de la fonction echelle. imech est une
+ime = imread('cb_31_03_a.jpg');
+[imech, ech] = echelle(ime); %appelle de la fonction echelle. imech est une
                             %image et ech un flottant (plus d'explication
                             %dans le fichier affinage.m)
@@ Ligne 46: / Ligne 48: @@
 image(imech) %affichage de l'image imech
-res = pos_deriv(ech, zx, 0, 0.1); %appelle de la fonction pos_deriv. res
+res = pos_deriv(ech, zx, 0.018, 0.066); %appelle de la fonction pos_deriv. res
                                   %est une matrice nx5 (plus d'explication
                                   %dans le fichier affinage.m)
@@ Ligne 53: / Ligne 55: @@
 plot(res(:,1), res(:,2)) %graphe du faisceau
-save('nom_fichier.txt', 'res', '-ascii'); %ecrit les donnes de res dans un
+title('Height of the beam as a function of x')
+xlabel('x (m)')
+ylabel('Height (m)')
+figure
+plot(res(:,2),res(:,4), '.')
+title('Refractive index as a function of height')
+xlabel('Height (m)')
+ylabel('Refractive index')
+%errorbar(res(:,4),res(:,5))
+% figure
+% ft = fit(res(:,1), res(:,2), 'poly9')
+% plot(ft)
+% save('fit_31_03_b.txt', 'ft');
+save('31_03_a_fit.txt', 'res', '-ascii'); %ecrit les donnes de res dans un
                                           %fichier .txt
 %-------------------------------------------------------------------------%
@@ Ligne 314: / Ligne 333: @@
     zx(2,i) = zx(2,i)*ech;
 end
 %-------------------------------------------------------------------------%
@@ Ligne 329: / Ligne 349: @@
         tab_donnees(j,1)=(i-1)*ech; %Remplissage des x
         tab_donnees(j,2)=zx(1,i); %Remplissage des z
-        tab_donnees(j,3)=(zx(1,i+10)-zx(1,i))/(10*ech); %derive spatiale
+        tab_donnees(j,3)=devfite((i-1)*ech); %derive spatiale du fit
-        tab_donnees(j,4)=sqrt(C0*(tab_donnees(j,3)^2+1)); %indice
+        %Ce fit a ete obtenu a l'aide des donnees contenu dans zx.
-        tab_donnees(j,5)=1/2*(zx(2,i)/zx(1,i)); %incertidue de l'indice
+        %La ligne pour l'obtenir est : ft = fit(zx(2,:),zx(1,:), 'poly9')
+        tab_donnees(j,4)=sqrt(C0*(devfite((i-1)*ech)^2+1)); %indice
+        %tab_donnees(j,5)=1/2*(zx2(2,i)/zx2(1,i)); %incertidue de l'indice
         j = j+1;
     end
@@ Ligne 340: / Ligne 362: @@
 </code>
+===== GRAPHES =====
+{{ :wiki:projet:traitement.jpg?nolink&800 |}}
+===== COMMENTAIRES =====
+<WRAP center round box 100%>
+La première figure représente l’image de départ, soit le faisceau laser dans la cuve.\\
+La deuxième est l’image que sort la fonction faisceau, le faisceau est sélectionné grâce à son intensité de rouge, puis on place les pixels valides sur une image initialement noire.\\
+La troisième est l'image que sort la fonction affinage, le faisceau est réduit à un pixel d’épaisseur. (remarque: cette fonction sort aussi les coordonnées de ces pixels ainsi que l’épaisseur de chaque abscisse)\\
+La quatrième figure est un graphe représentant l'altitude du laser en fonction de la position en abscisse dans la cuve.\\
+Enfin la dernière image est le graphe de l’indice en fonction de l'altitude.\\
+Lors de l'affinage du faisceau, on obtient une altitude approximative en fonction de sa position en abscisse en pixels, après une conversion en mètre, on obtient une altitude. Celle ci est quasiment des morceaux de constante, rendant imprécis le calcul de l'indice qui utilise la dérivée spatiale de l'altitude. On va donc réaliser un "fit" de cette altitude à l'aide d'une approximation d'un polynôme d'ordre 9 (Un tel ordre est sûrement trop grand, mais aux vus de nos approximations et erreurs parfois grossières, cela nous permet de nous exonérer d'une source d'erreur en plus). On obtient alors un graphe pour l'indice qui ne semble pas incohérent.
+</WRAP>

FablabSU

Outils pour utilisateurs

Outils du site

Différences

Outils de la page