PLOT_HARDI

 plot ngdir original directions in which the HADC was simulated and 900 interpolated directions
 
 Input: 
     1. points (ngdir directions)
     2. ADC_allcmpts_alldir
     
 Output:
     1 figure with title of "ADC in ngdir directions"

0001 function PLOT_HARDI(points,ADC_allcmpts_alldir)
0002 
0003 % plot ngdir original directions in which the HADC was simulated and 900 interpolated directions
0004 %
0005 % Input:
0006 %     1. points (ngdir directions)
0007 %     2. ADC_allcmpts_alldir
0008 %
0009 % Output:
0010 %     1 figure with title of "ADC in ngdir directions"
0011 
0012 ngdir = size(points,1);
0013 
0014 [sph_pts,C_sph] = spheresurface_regularpoints(1,900);   
0015 YY = spherical_harmonics(points(:,1),points(:,2),points(:,3),ones(size(points,1),1));
0016 sh_coeff = YY\ADC_allcmpts_alldir;
0017     
0018 YY_sph = spherical_harmonics(sph_pts(:,1),sph_pts(:,2),sph_pts(:,3),ones(size(sph_pts,1),1));
0019 ADC_interp = YY_sph*sh_coeff;
0020 ADC_alldir = ADC_interp.*sph_pts;
0021         
0022 figure; hold on;
0023 ADC_Ellipsoid = trisurf(C_sph,ADC_alldir(:,1),ADC_alldir(:,2),ADC_alldir(:,3),ADC_interp); view(3); colorbar; axis equal;
0024 ADC_Ellipsoid.EdgeColor = 'none';
0025 xlabel('x'); ylabel('y'); zlabel('z');
0026 title(['ADC in ',num2str(ngdir),' directions']);
0027 plot3(ADC_allcmpts_alldir.*points(:,1),ADC_allcmpts_alldir.*points(:,2),ADC_allcmpts_alldir.*points(:,3),'k.',...
0028 'markersize',30);