mtexdata forsterite
region = [5 2 10 5]*10^3;
close all
plot(ebsd)
rectangle('position',region,'edgecolor','r','linewidth',2)
condition = inpolygon(ebsd,region);
ebsd = ebsd(condition)
figure
plot(ebsd)
If it doesn't work, check also number 2 here: https://groups.google.com/forum/#!topic/mtexmail/adn2k7Vpt70
Best,
Filippe
ebsd=ebsd(condition)
ebsd
mtexdata forsterite
figure; plot(ebsd)
x1=5*10^3;
y1=2*10^3;
x2=15*10^3;
y2=7*10^3;
ebsd=ebsd(ebsd.prop.x>x1 & ebsd.prop.y>y1 & ebsd.prop.x<x2 & ebsd.prop.y<y2)
figure; plot(ebsd)