processing:complex-multi-phase
This is an old revision of the document!
Table of Contents
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.1559894560.txt.gz · Last modified: by matthias
