loadPoleFigure_geesthacht.m

function pf = loadPoleFigure_geesthacht(fname,varargin)
% load ptx file
%
% Input
%  fname - file name
%
% Output
%  pf    - vector of @PoleFigure
%
% See also
% ImportPoleFigureData loadPoleFigure

% read data

try
  fid = efopen(fname);
  
  d = textscan(fid,'%n %n %*n %n %n %n %*n %*n %*n %*n %*n:%*n',...
    'commentStyle','$','HeaderLines',1);
  
  fclose(fid);
  
  % identify data
  assert(all(d{1}.'==(1:length(d{1}))));
  d = {d{2:5}};
  assert(~isempty(d{1})); % all items have the same length! (equal(cellfun('length',d),1))
  assert(all(d{2}>=0 & d{2}<=90 & d{3}>=-370 & d{3}<=370));
  
  pos = 1;
  r = {};
  h = {};
  
  intensities = {};
  while pos <= length(d{1})
    
    npos = find(d{1}(pos)==d{1},1,'last');
    r{end+1} = vector3d.byPolar(pi/2-d{2}(pos:npos)*degree,d{3}(pos:npos)*degree); %#ok<AGROW>
    assert(length(r{1})>=5);
    intensities{end+1} = d{4}(pos:npos); %#ok<AGROW>
    h{end+1} = Miller(1,0,0,crystalSymmetry('m-3m')); %#ok<AGROW>
        
    pos = npos+1;
    
  end
  
  hh = getClass(varargin,'Miller');
  if length(hh) == length(h)
    if iscell(hh)
      h = hh; 
    else
      h = vec2cell(hh); 
    end
  end
  
  pf = PoleFigure(h,r,intensities,varargin{:});
  
catch
  interfaceError(fname);
end