function voxelList = readRoiImage()

Vm = spm_vol(spm_select([1 Inf],'image','Select ROI image'));
nRoiFiles = size(Vm,2);

% for iRoiFile = 1:nRoiFiles
%     x = []; y = []; z = [];
%     [x1 y1] = ndgrid(1:V(k).dim(1),1:V(k).dim(2));
%     for p = 1:V(k).dim(3) % resample mask Vm(iRoiFile) in space of V(k)
%         B = spm_matrix([0 0 -p 0 0 0 1 1 1]);
%         M = inv(B*inv(V(k).mat)*Vm(iRoiFile).mat);
%         msk = find(spm_slice_vol(Vm(iRoiFile),M,V(k).dim(1:2),0));
%         if ~isempty(msk)
%             z1 = p*ones(size(msk(:)));
%             x = [x; x1(msk(:))];
%             y = [y; y1(msk(:))];
%             z = [z; z1];
%         end;
%     end;
% 
% end
voxelList = Vm;

end