Hi Peyman,
yes, you can calculate a kam using the function KAM. I'm not entirely sure what you mean with "for different orientations" since [0001] etc are crystal directions, however assuming you have an orientation as a reference, you can select the ebsd of the grains
which have a meanOrientation within a certain threshold, e.g. 10degrees to this reference orientation and calculate the kam for this subset.
Assuming your reference orientation would be o and grains would be of a single phase,conceptually this would be:
kam = KAM(ebsd(grains(angle(o,grains.meanOrientation) < 10*degree)))
Of course you could similarly select a subset of grains e.g. with just a crystal direction pointing towards some specimen direction, e.g. 0001 close to parallel to z would be:
kam = KAM(ebsd(grains(angle(zvector, grains.meanOrientation.*Miller(0,0,0,1,cs)) < 10*degree)))
I hope this helps,
Cheers, Rüdiger