--
If you want to reduce the number of emails you get through this forum login to https://groups.google.com/forum/?fromgroups=#!forum/mtexmail, click "My membership" and select "Don't send me email updates". You can still get emails on selected topics by staring them.
---
You received this message because you are subscribed to the Google Groups "MTEX" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mtexmail+u...@googlegroups.com.
Visit this group at http://groups.google.com/group/mtexmail.
For more options, visit https://groups.google.com/d/optout.
plotIPDF(ori,zvector,'contourf')%% Import Script for ODF Data%% This script was automatically created by the import wizard. You should% run the whoole script or parts of it in order to import your data. There% is no problem in making any changes to this script.
%% Specify Crystal and Specimen Symmetries
% crystal symmetryCS = crystalSymmetry('m-3m', [3.615 3.615 3.615], 'mineral', 'Copper', 'color', 'light blue');
% specimen symmetrySS = specimenSymmetry('1');
% plotting conventionsetMTEXpref('xAxisDirection','east');setMTEXpref('zAxisDirection','outOfPlane');
%% Specify File Names
% path to files%pname = 'C:\Users\ItaloP\Documents\AADottorato\Lavoro-Andrew\VPSCsimulations\Verifica_texture_independence_with_hardening_parameters\tensile_test\param_1';
% which files to be importedfname = ['TEX_PH1.OUT'];
%% Import the Data
% % specify kernel% psi = deLaValeePoussinKernel('halfwidth',13*degree);% % % create an EBSD variable containing the data% odf = loadODF(fname,CS,SS,'density','kernel',psi,'resolution',5*degree,...% 'interface','VPSC');
ebsd = loadEBSD_generic(fname,'ColumnNames',{'Euler1','Euler2','Euler3'},'HEADER',4);
%% plot inverse pole figures %Specify Miller Indice
m1 = Miller(0,0,1,CS);m2 = Miller(1,0,1,CS);m3 = Miller(1,1,1,CS);
%Specify vector
v1 = vector3d(1,0,0);v2 = vector3d(0,1,0);v3 = vector3d(0,0,1);
%% plot IPF from EBSD discrete
figure(01)plotIPDF(ebsd.orientations,v1)annotate([Miller(m1),Miller(m2),Miller(m3)],'all','labeled')saveas(figure(01),'IPF from EBSD discrete.jpeg')
%% plot IPF from EBSD contourfigure(02)plotIPDF(ebsd.orientations,v1,'contourf')colorbarannotate([Miller(m1),Miller(m2),Miller(m3)],'all','labeled')saveas(figure(02),'IPF from EBSD contour.jpeg')
%% Import the Data
% specify kernelpsi = deLaValeePoussinKernel('halfwidth',13*degree);
% create an EBSD variable containing the dataodf = loadODF(fname,CS,SS,'density','kernel',psi,'resolution',5*degree,... 'interface','VPSC');
%% plot inverse pole figures from ODF
figure(03)plotIPDF(odf, v3,'contour','antipodal','resolution',5*degree) %'contour',0:0.2:3.2, 'minmax', 'off')annotate([Miller(m1),Miller(m2),Miller(m3)],'all','labeled')colorbarsaveas(figure(03),'IPF from ODF.jpeg')Hi,
carefully check your reference direction for IPDF:
In your script it is v1 for figure 1 and 2, but v3 for figure 3.
I was getting crazy for this stupid mistake. I am really sorry to have wasted your time.
One last question if I can, now the qualitative results are the same for both the IPF from the EBSD and ODF. The problem still remain for the quantitative results, which are different. Is it a good method, in your opinion, to play with the halfwidth and make the value of density for the IPF from ODF equal to the one of the IPF from EBSD?
Is it a good method, in your opinion, to play with the halfwidth and make the value of density for the IPF from ODF equal to the one of the IPF from EBSD?
plotIPDF(ori,zvector,'contourf','halfwidth',10*degree)