mtex-toolbox / mtex

MTEX is a free Matlab toolbox for quantitative texture analysis. Homepage:
http://mtex-toolbox.github.io/
GNU General Public License v2.0
281 stars 185 forks source link

Very large data files ODF and pole figures and memory problems #149

Closed DavidMainprice closed 8 years ago

DavidMainprice commented 8 years ago

I am doing some standard calculations of ODFs and Pole figures. However the EBSD object has over 17 million indexed points and ebsd_subregion has over 7 million index points.

So when use the script below I get the the following error

Error using subsref Index exceeds matrix dimensions.

Error in rotation/subsref (line 6) r.a = subsref(r.a,s(1));

Error in symmetry/subsref (line 8)

So we have run out of space. I using Mac, but I could try running it on PC as memory management is different between the two systems, Linux I do not know too well.

Have someone any suggestions for any operating system about increasing the memory for such applications. Does anyone know the current limit for indexed points in MTEX applications.

all the best David

%% Estimate ODF of selected sub-area Ice Fabric Analyser measurements %*** % Calculate ODF with Standard de Vallee Poussin kernel with 10 degree HW odf_Ice = calcODF(ebsd_subregion('Ice Ih').orientations,'halfwidth',10_degree) % plot pole figure % make list of pole figures, this case only (0001) PFs_Ice = Miller(0,0,0,1,CS('Ice Ih'),'hkil'); % Plot Ice Pole Figure (Lower hemisphere) close all figure plotpdf(odf_Ice,PFs_Ice,'smooth','antipodal','lower','resolution',1_degree); mtexColorbar drawNow(gcm,'figSize','small') % file name saveFigure('/MatLab_Programs/PFs_LH_Ice.png');

DavidMainprice commented 8 years ago

With correct syntax it works even 17 million index points - sorry my error

all the best David

Corrected syntax below

% %% Estimate ODF of selected sub-area Ice Fabric Analyser measurements % % Calculate ODF with Standard de Vallee Poussin kernel with 10 degree HW odf_Ice = calcODF(ebsd_subregion('Ice Ih').orientations,'halfwidth',10_degree) % plot pole figure % make list of pole figures, this case only (0001) PFs_Ice = Miller(0,0,0,1,ebsd_subregion('Ice Ih').CS,'hkil'); % Plot Ice Pole Figure (Lower hemisphere) close all figure plotPDF(odf_Ice,PFs_Ice,'smooth','antipodal','lower','resolution',1_degree); mtexColorbar drawNow(gcm,'figSize','small') % file name saveFigure('/MatLab_Programs/PFs_LH_Ice.png'); %