User Tools

Site Tools


processing:complex-multi-phase

This is an old revision of the document!


In the works: indexing a complex dataset with many phases

If your dataset is complex, here is a list of tricks to get it indexed. This page is being improved over time.

Trick 1: plot intensity vs. 2theta from peak histogram in ImageD11:

from ImageD11 import columnfile
c = columnfile.columnfile('peaks_t100.flt') % if 'expection: problem interpreting your colfile' check the extention of your file. Should be .flt
c.titles % to see what variables you have in your columnfile
c.parameters.loadparameters(“CeO2_parameters.prm”)
c.titles
c.updateGeometry()
import pylab
pylab.figure()
pylab.show()
import matplotlib
matplotlib.use(“GTK3Agg”)
from pylab import *
show()
from ImageD11.columnfile import *
c = columnfile(“peaks_t100.flt”)
c.parameters.loadparameters(“CeO2_parameters.prm”)
c.updateGeometry()
tth=arange(0,15,.01)
plot(np.histogram(c.tth,bins=tth)[0],“-”)
show()
plot(tth[1:],np.histogram(c.tth,bins=tth)[0],“-”) % give a better peaks separation
show()
plot(tth[1:],np.histogram(c.tth,bins=tth,weights=c.sum_intensity)[0],“-”) % take in account the intensity of the peaks

% open a python console
python

import matplotlib
matplotlib.use(“GTK3Agg”)
from pylab import *
from ImageD11.columnfile import *
c = columnfile('peaks_t100.flt')
c.parameters.loadparameters(“CeO2_parameters.prm”)
c.updateGeometry()
tth = arange(0,15,.01)
plot(tth[1:],np.histogram(c.tth,bins=tth)[0],“-”)
show()

Trick 2: work with large grain first

processing/complex-multi-phase.1559894251.txt.gz · Last modified: 2019/06/07 07:57 by matthias