matlab怎么读取Img格式的文件,60行报错

猪猪虾

1. 
   
2. function [t, files] = get_totals_revised_by_long(files, thr, msk)
   
3. %get_totals - Returns image totals (sum over all voxels), in ml
   
4. %  t = get_totals
   
5. %  [t files] = get_totals(files, thr, msk)
   
6. % GUI file-selection is used if files not specified as argument (or empty).
   
7. %
   
8. % If thr is given, this will be treated as an absolute threshold
   
9. % (i.e. values below this will be zeroed, hence the total will better match
   
10. % the GM analysed in
    
11. 
    
12. % the voxelwise stats, with the same threshold masking).
    
13. %
    
14. % Similarly, if msk is specified this image will be used as an explicit
    
15. % mask (i.e. only non-zero mask voxels will be included).
    
16. % GUI file-selection is used if msk is given as empty string ('').
    
17. % [Currently, masking assumes that msk matches the voxel dimensions of each
    
18. % image, and that therefore, all images have the same dimensions.]
    
19. 
    
20. % check spm version:
    
21. if exist('spm_select','file') % should be true for spm5
    
22.     spm5 = 1;
    
23.     select = @(msg) spm_select(inf, 'image', msg);
    
24. elseif exist('spm_get','file') % should be true for spm2
    
25.     spm5 = 0;
    
26.     select = @(msg) spm_get(inf, 'img', msg);
    
27. else
    
28.     error('Failed to locate spm_get or spm_select; please add SPM to Matlab path')
    
29. end
    
30. 
    
31. if ( ~exist('files', 'var') || isempty(files) )
    
32.     files = select('choose images');
    
33. end
    
34. if ( ~exist('thr', 'var') || isempty(thr) )
    
35.     thr = -inf; % default to include everything (except NaNs)
    
36. end
    
37. if ~exist('msk', 'var')
    
38.     msk = 1; % default to include everything
    
39. end
    
40. if isempty(msk)
    
41.     msk = select('Choose mask image');
    
42. end
    
43. if ischar(msk)
    
44.     msk = spm_vol(msk);
    
45. end
    
46. if isstruct(msk)
    
47.     msk = spm_read_vols(msk);
    
48. end
    
49. msk = msk ~= 0;
    
50. vols = spm_vol(files);
    
51. N = length(vols);
    
52. 
    
53. the_number_of_mask=1;  %%
    
54. HIV_1234vol=zeros(N, the_number_of_mask);%%dingyimubiao
    
55. 
    
56. t = zeros(N,1);
    
57. for n = 1:N
    
58.     hwait=waitbar(n/N,'请等待>>>>>>>>');
    
59.     for nn=1:the_number_of_mask
    
60.         Imask=load_nii(strcat('F:\data\t\',num2str(nn),'myMask.img'));   %% add the mask files
    
61.         Imask=double(Imask.img);
    
62.         msk=Imask;
    
63.         vsz = abs(det(vols(n).mat));
    
64.         img = spm_read_vols(vols(n));
    
65.         img = img .* msk;
    
66.         t(n) = sum(img(img > thr)) * vsz / 1000; % vsz in mm^3 (= 0.001 ml)
    
67.         HIV_1234vol(n,nn)=t(n);  %%fuzhi
    
68.     end
    
69.     close(hwait);
    
70. end
    
71. save('HIV_1234vol.mat','HIV_1234vol')
    
72. 
    
73. 
    

复制代码