spiread.m

matlab处理图像的一些基本方法。其中有一部分mex程序需要安装编译

function [img, ftime]=spiread(fname)
%SPIREAD Read the image from Spider format.
%
%       IMAGE = SPIREAD(FILENAME)
%
%       IMAGE is the image matrix in nm.
%       FILENAME the spider file name.
%
%       On error, an error message is displayed.
%
%       Claudio 5 Apr, 1995.

%
% Copyright (c) 1995 by Claudio Rivetti and Mark Young
% [e-mail],    [e-mail]
%

if nargin < 1
  error('Too few input arguments.');
end

OFFSET=2048;


[fp, errmsg]=fopen(fname, 'r');

if fp==-1
  error(errmsg);
end

s=fread(fp, 2, 'float');
SIZE=[s(2),s(2)];

fseek(fp, 21*4, 'bof');
OFFSET=fread(fp, 1, 'float');

fseek(fp, 211*4, 'bof');
d=fread(fp,12, 'char');
d=setstr(d');

t=fread(fp,8, 'char');
t=setstr(t');
ftime=[t, '  ' d];

fseek(fp, OFFSET, 'bof');
img=fread(fp, SIZE, 'float');
img = fliplr(img);
img = rot90(img, 1);

fclose(fp);

return;