function [DIM,VOX,SCALE,TYPE,OFFSET,ORIGIN,DESCRIP] = biac_hread(P)
% reads a header
% FORMAT [DIM VOX SCALE TYPE OFFSET ORIGIN DESCRIP] = biac_hread(P);
%
% P       - filename 	     (e.g spm or spm.img)
% DIM     - image size       [i j k [l]] (voxels)
% VOX     - voxel size       [x y z [t]] (mm [secs])
% SCALE   - scale factor
% TYPE    - datatype (integer - see spm_type)
% OFFSET  - offset (bytes)
% ORIGIN  - origin [i j k]
% DESCRIP - description string
%___________________________________________________________________________
%
% spm_hread reads variables into working memory from a SPM/ANALYZE
% compatible header file.  If the header does not exist global defaults
% are used.  The 'originator' field of the ANALYZE format has been
% changed to ORIGIN in the SPM version of the header.  funused1 of the
% ANALYZE format is used for SCALE
%
% see also dbh.h (ANALYZE) spm_hwrite.m and spm_type.m
%
%__________________________________________________________________________
% Edited by Josh Bizzell, 2003-April-04, for the BIAC
% @(#)spm_hread.m	2.7 99/10/29

P = deblank(P);
P2 = [P(1:end-3),'hdr'];
if (strncmp(P(end-2:end),'img',3) & exist(P2,'file')) | ...
    strncmp(P(end-2:end),'hdr',3)
  [DIM,VOX,SCALE,TYPE,OFFSET,ORIGIN,DESCRIP] = spm_hread(P);
else
  V = readmrtest(P,'=>INFOONLY'); 
  readmrtest(V,'=>CLEANUP');
  if ~isempty(V)
    DIM = [getfield(V.info.dimensions,{1},'size'),...
        getfield(V.info.dimensions,{2},'size'),...
        getfield(V.info.dimensions,{3},'size')];
    VOX = [getfield(V.info.dimensions,{1},'spacing'),...
        getfield(V.info.dimensions,{2},'spacing'),...
        getfield(V.info.dimensions,{3},'spacing')];
    SCALE = 1;
    TYPE = spm_type(getfield(V.info,'elemtype'));
    OFFSET = 0;
    % Maybe set to actual origin
    ORIGIN = [0 0 0];
    DESCRIP = ['BIAC - ',P];
  else
    error(sprintf('Error reading %s.',P));
    global DIM VOX SCALE TYPE OFFSET ORIGIN
    DESCRIP = ['defaults'];
  end
end