How to handle the grainid beyond the number of grains

[kun-Jiang]

Hello everyone,

I'm currently working on calculating the dislocation density around the crack path, and I've been following a tutorial to achieve this. However, I encountered an error when running the following command:
Command:
ipfKey.oriRef = grains(ebsd('Ni (Nickel)').grainId).meanOrientation;
The error message states:
Out of range; index must be a positive integer or logical.
I checked the code and find that the maximum id obtained from ebsd('Ni (Nickel)').grainId is 1503, but there is only 477 grains. This inconsistency has left me puzzled. So how did this problem arise?

Following is my code and structure:

[zmichels]

Does it work if you add 'id' in the command?

Like:

ipfKey.oriRef = grains('id', ebsd('Ni (Nickel)').grainId).meanOrientation

[kun-Jiang]

Thanks, zmichels, but it still doesn't work. This is the error message
No data with the specified ids in the data set

[zmichels]

What about if you add (after filtering out some grains and before the problematic line):

ebsd = ebsd(grains);

[kun-Jiang]

Thanks, zmichels. That works well. So what is the meaning of this command?

[zmichels]

@kun-Jiang I admit I'm not entirely sure why it's necessary, but somehow the id affiliations are confused when there are EBSD points for grains that have been removed from the grainset. Whenever I filter out grains I do this to avoid similar issues you had. I didn't notice the first time looking through the code tho.