Plotting the intensity of one pole/reflection as a function of the polar angle

[Martin Lentz]

• Please be as detailed as possible and answer the following questions.
  
  1. What input kind data do you use?
  Diffraction data and EPSC texture output files
  2. How do you import them?
  I'm using the import scripts
  3. What do you want to compute from your data?

4. How to you want to plot or export you result?

I would like to calculate and plot the intensity/density of one pole (e.g. {0002}) as a function of the polar angle

integrating over the azimuth angle
5. What MTEX version do you use?
3.3.1
If possible attach simplified m-files and data files.

[Ralf Hielscher]

Taking means over circles of the pole figure is not supported by a single command. However, assuming you have given an odf you can do the following:

define a regular grid where you are going to evaluate the pole figure

r = S2Grid('regular')

evaluate the odf at this grid

pf = pdf(odf,Miller(0,0,1),r);

rehape it to have the azimuth angles as columns and the polar angles as rows

pf = reshape(pf,size(r)).' ;

take the mean over the azimuth angles

mean_pf = mean(pf);

plot the result

plot(mean_pf)