豆粑粑 3D sphere画scatter3散点图

meatball1982

别人的问题，
要画散点图。
只用scatter3，无法调alpha

有一个scatter3sph的函数，思路差不多，一个一个的画球体，可以设置color和size，还有transparent.
比较了两种方式。
都还不错。

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. load data.mat
   
6. x=XYZ2(:,1);
   
7. y=XYZ2(:,2);
   
8. z=XYZ2(:,3);
   
9. v=sqrt(sum((XYZ2(:,[4:6]).^2)'))';
   
10. ind=v<20;
    
11. x(ind)=[];
    
12. y(ind)=[];
    
13. z(ind)=[];
    
14. v(ind)=[];
    
15. n=length(v);
    
16. 
    
17. tm=flipud(jet(140));
    
18. col_tm=tm(5:end-6,:);
    
19. 
    
20. [ba,ind]=sort(v);
    
21. for i=1:n
    
22.     col_mm(i,:)=col_tm(floor(v(i))-19,:);
    
23. end
    
24. 
    
25. % sbp_width=0.8;
    
26. % sbp_heig=0.8;
    
27. % n_row = 2;
    
28. % n_col = 2;
    
29. % [out_pos]=fun_mm_subplot_pos(n_row,n_col,sbp_width,sbp_heig);
    
30. % h=figure(1)
    
31. % % set(h, 'Position', [100, 100, 1100, 1000]);
    
32. % set(h, 'Position', [100, 100, 800, 800]);
    
33. % % %potential
    
34. % % ax=axes('position',out_pos(1,:));
    
35. % % scatter3sph(x,y,z,'color',col_mm,'size',log(log(log(v)))/100,'Trans',0.6)
    
36. % scatter3sph(x,y,z,'color',col_mm,'size',log(log(v))/300,'Trans',0.6)
    
37. % grid on
    
38. % axis tight
    
39. % xlabel('x');
    
40. % ylabel('y');
    
41. % zlabel('z');
    
42. % view(30,20)
    
43. 
    
44. % ax=axes('position',out_pos(2,:));
    
45. % scatter3sph(x,y,z,'color',col_mm,'size',log(log(log(v)))/100,'Trans',0.9)
    
46. % grid on
    
47. % xlabel('x');
    
48. % ylabel('y');
    
49. % zlabel('z');
    
50. % view(0,90)
    
51. %
    
52. % ax=axes('position',out_pos(3,:));
    
53. % scatter3sph(x,y,z,'color',col_mm,'size',log(log(log(v)))/100,'Trans',0.9)
    
54. % grid on
    
55. % xlabel('x');
    
56. % ylabel('y');
    
57. % zlabel('z');
    
58. % view(0,0)
    
59. %
    
60. % ax=axes('position',out_pos(4,:));
    
61. % scatter3sph(x,y,z,'color',col_mm,'size',log(log(log(v)))/100,'Trans',0.9)
    
62. % grid on
    
63. % xlabel('x');
    
64. % ylabel('y');
    
65. % zlabel('z');
    
66. % view(90,0)
    
67. 
    
68. % h=gcf;
    
69. % fig_na=['./fig_sphere_3D_single'];
    
70. % fun_work_li_035_myfig_out(h,fig_na,3)
    
71. 
    
72. 
    
73. sbp_width=0.8;
    
74. sbp_heig=0.8;
    
75. n_row = 2;
    
76. n_col = 2;
    
77. [out_pos]=fun_mm_subplot_pos(n_row,n_col,sbp_width,sbp_heig);
    
78. h=figure(1)
    
79. % set(h, 'Position', [100, 100, 1100, 1000]);
    
80. set(h, 'Position', [100, 100, 800,800]);
    
81. 
    
82. 
    
83. maxax= max([range(x), range(y), range(z)]);
    
84. ratios= [range(x)/maxax, range(y)/maxax, range(z)/maxax];
    
85. 
    
86. % %potential
    
87. % ax=axes('position',out_pos(1,:));
    
88. colormap(col_tm)
    
89. h=scatter3(x,y,z,(v)*3,v,...
    
90.     'filled','markeredgecolor','c')
    
91. view(30,20)
    
92. daspect([ratios(1), ratios(2), ratios(3)]);
    
93. xlabel('x');
    
94. ylabel('y');
    
95. zlabel('z');
    
96. axis tight
    
97. %
    
98. % % ax=axes('position',out_pos(2,:));
    
99. % % colormap(col_tm)
    
100. % % h=scatter3(x,y,z,(v)*3,v,...
     
101. % %     'filled','markeredgecolor','c')
     
102. % % view(0,90)
     
103. % % daspect([ratios(1), ratios(2), ratios(3)]);
     
104. % % xlabel('x');
     
105. % % ylabel('y');
     
106. % % zlabel('z');
     
107. % %
     
108. % %
     
109. % % ax=axes('position',out_pos(3,:));
     
110. % % colormap(col_tm)
     
111. % % h=scatter3(x,y,z,(v)*3,v,...
     
112. % %     'filled','markeredgecolor','c')
     
113. % % view(0,0)
     
114. % % daspect([ratios(1), ratios(2), ratios(3)]);
     
115. % % xlabel('x');
     
116. % % ylabel('y');
     
117. % % zlabel('z');
     
118. % %
     
119. % %
     
120. % % ax=axes('position',out_pos(4,:));
     
121. % % colormap(col_tm)
     
122. % % h=scatter3(x,y,z,(v)*3,v,...
     
123. % %     'filled','markeredgecolor','c')
     
124. % % view(90,0)
     
125. % % daspect([ratios(1), ratios(2), ratios(3)]);
     
126. % % xlabel('x');
     
127. % % ylabel('y');
     
128. % % zlabel('z');
     
129. % %
     
130. % %
     
131. h=gcf;
     
132. fig_na=['./fig_circle_3D_sing'];
     
133. fun_work_li_035_myfig_out(h,fig_na,3)
     
134. % % subplot(1,2,2)
     
135. % colormap(col_tm)
     
136. % h=scatter3(x,y,z,(v)*3,v,...
     
137. %     'filled','markeredgecolor','c')
     
138. % maxax= max([range(x), range(y), range(z)]);
     
139. % ratios= [range(x)/maxax, range(y)/maxax, range(z)/maxax];
     
140. % view(0,90)
     
141. % daspect([ratios(1), ratios(2), ratios(3)]);
     
142. 
     
143. 
     
144. % scatter3sph(x(ind),y(ind),z(ind),'color',col_mm,'size',log(v(ind))/1000)
     

复制代码

1. function scatter3sph(X,Y,Z,varargin)
   
2. %SCATTER3SPH (X,Y,Z) Makes a 3d scatter plot with 3D spheres
   
3. %        SCATTER3SPH is like scatter3 only drawing spheres instead
   
4. %                of flat circles, at coordinates specified by vectors X, Y, Z. All three
   
5. %                vectors have to be of the same length.
   
6. %        SCATTER3SPH(X,Y,Z) draws the spheres with the default size and color.
   
7. %        SCATTER3SPH(X,Y,Z,'size',S) draws the spheres with sizes S. If length(S)= 1
   
8. %                the same size is used for all spheres.
   
9. %        SCATTER3SPH(X,Y,Z,'color',C) draws the spheres with colors speciffied in a
   
10. %                N-by-3 matrix C as RGB values.
    
11. %        SCATTER3SPH(X,Y,Z,'transp',T) applies transparency level 'T' to the spheres
    
12. %                T= 0 => transparent, T= 1 => opaque.
    
13. %        Parameter names can be abreviated to 3 letters. For example: 'siz' or
    
14. %                'col'. Case is irrelevant.
    
15. %
    
16. % Example
    
17. % %Coordinates
    
18. %  X= 100*rand(9,1); Y= 100*rand(9,1); Z= 100*rand(9,1);
    
19. %
    
20. % %Colors: 3 blue, 3 red and 3 green
    
21. % C= ones(3,1)*[0 0 1];
    
22. % C= [C;ones(3,1)*[1 0 0]];
    
23. % C= [C;ones(3,1)*[0 1 0]];
    
24. %
    
25. % %Sizes
    
26. % S= 5+10*rand(9,1);
    
27. %
    
28. % scatter3sph(X,Y,Z,'size',S,'color',C,'trans',0.3);
    
29. % axis vis3d
    
30. 
    
31. 
    
32. %-- Some checking...
    
33. if nargin < 3 error('Need at least three arguments'); return; end
    
34. if mean([length(X),length(Y),length(Z)]) ~= length(X) error ('Imput vectors X, Y, Z are of different lengths'); return; end
    
35. 
    
36. %-- Defaults
    
37. C= ones(length(X),1)*[0 0 1];
    
38. S= 0.1*max([X;Y;Z])*ones(length(X),1);
    
39. nfacets= 15;
    
40. transp= 0.5;
    
41. 
    
42. 
    
43. %-- Extract optional arguments
    
44. for j= 1:2:length(varargin)
    
45.         string= lower(varargin{j});
    
46.         switch string(1:min(3,length(string)))
    
47.                 case 'siz'
    
48.                         S= varargin{j+1};
    
49.                         if length(S) == 1
    
50.                                 S= ones(length(X),1)*S;
    
51.                         elseif length(S) < length(X)
    
52.                                 error('The vector of sizes must be of the same length as coordinate vectors (or 1)');
    
53.                                 return
    
54.                         end
    
55. 
    
56.                 case 'col'
    
57.                         C= varargin{j+1};
    
58.                         if size(C,2) < 3        error('Colors matrix must have 3 columns'); return; end
    
59.                         if size(C,1) == 1
    
60.                                 C= ones(length(X),1)*C(1:3);
    
61.                         elseif size(C,1) < length(X)
    
62.                                 error('Colors matrix must have the same number of rows as length of coordinate vectors (or 1)');
    
63.                                 return
    
64.                         end
    
65. 
    
66.                         case 'fac'
    
67.                                 nfacets= varargin{j+1};
    
68. 
    
69.                         case 'tra'
    
70.                                 transp= varargin{j+1};
    
71. 
    
72.                 otherwise
    
73.                         error('Unknown parameter name. Allowed names: ''size'', ''color'', ''facets'', ''transparency'' ');
    
74.         end
    
75. end
    
76. 
    
77. 
    
78. %-- Sphere facets
    
79. [sx,sy,sz]= sphere(nfacets);
    
80. 
    
81. 
    
82. %--- Correct potential distortion
    
83. maxax= max([range(X), range(Y), range(Z)]);
    
84. ratios= [range(X)/maxax, range(Y)/maxax, range(Z)/maxax];
    
85. sx= sx*ratios(1);
    
86. sy= sy*ratios(2);
    
87. sz= sz*ratios(3);
    
88. 
    
89. 
    
90. %-- Plot spheres
    
91. hold on
    
92. for j= 1:length(X)
    
93.         surf(sx*S(j)+X(j), sy*S(j)+Y(j), sz*S(j)+Z(j),...
    
94.                 'LineStyle','none',...
    
95.                 'FaceColor',C(j,:),...
    
96.                 'FaceAlpha',transp);
    
97. end
    
98. 
    
99. daspect([ratios(1), ratios(2), ratios(3)]);
    
100. light('Style','infinit','Color',[1  1  1]);
     
101. % light('Position',[1 1 1],'Style','infinit','Color',[1 1 1]);
     
102. % lighting gouraud;
     
103. % view(30,30)
     

复制代码