I see that VO.fname is empty. Perhaps specify a file name and try again.
Best regards,
-John
On 26 January 2012 09:52, Bonobo 2000 <[log in to unmask]> wrote:
> Hi,
>
> I am trying to save a volume using
>
> spm_write_vol( VO , Volume);
>
> where VO is
>
> VO(1) = struct(...
> 'fname', [],...
> 'dim', vars.VolDimen,...
> 'mat', vars.Mat,...
> 'pinfo', [1 0 0]',...
> 'descrip', 'SearchLightMod: ',...
> 'Sphere', vars.Sphere);
>
> and Volume is my data.
>
> This procedure works fairly good, but occasionally I am getting the
> following error:
>
> There was a problem writing to the header of
> "XXXXXXXXXXXXXXXX"
> ok =
> 0
> ??? Error using ==> nifti.create>create_each at 29
> Unable to write header for
>
>
> I tracked the error upto the function ok = write_hdr_raw(fname,hdr,be)
> more precisely here:
>
> for i=1:length(org)
> if isfield(hdr,org(i).label),
> dat = hdr.(org(i).label);
> if length(dat) ~= org(i).len,
> if length(dat)< org(i).len,
> dat = [dat(:) ; zeros(org(i).len-length(dat),1)];
> else
> dat = dat(1:org(i).len);
> end;
> end;
> else
> dat = org(i).def;
> end;
> % fprintf('%s=\n',org(i).label)
> % disp(dat)
> len = fwrite(fp,dat,org(i).dtype.prec);
> if len ~= org(i).len,
> ok = false;
> end;
> end;
> fclose(fp);
> if ~ok,
> fprintf('There was a problem writing to the header of\n');
> fprintf('"%s"\n', fname);
> end;
>
>
> But I have no idea what is happening here, may be somebody could help me out
> there.
>
> Best,
> Jane
|