Mathematica中文论坛-非官方

标题: 豆粑粑 matlab 4D slice transparent [打印本页]

---

作者: meatball1982    时间: 2017-9-20 09:11
标题: 豆粑粑 matlab 4D slice transparent

1. clear all
   
2. clc
   
3. close all
   
4. 
   
5. %% outline
   
6. % plot scatter points, surf, countour3D
   
7. 
   
8. %% main
   
9. % load data ---------------------------------------------------------------
   
10. % dat=load('./cv_357.txt');
    
11. % n=length(dat(:,1));
    
12. %
    
13. % x=dat(:,1);
    
14. % y=dat(:,2);
    
15. % z=dat(:,3);
    
16. %
    
17. % % grid size ---------------------------------------------------------------
    
18. % n_cv=30;
    
19. % cv_g=linspace(-pi,pi,n_cv);
    
20. %
    
21. % cv_int=0.5; % 2 calculate point denisty, half of cubic size
    
22. % cv_x=cv_g;
    
23. % cv_y=cv_g;
    
24. % cv_z=cv_g;
    
25. %
    
26. % [cv_X,cv_Y,cv_Z]=meshgrid(cv_g,cv_g,cv_g);
    
27. %
    
28. % % points denisty ----------------------------------------------------------
    
29. % for i=1:n_cv
    
30. %     ind_x= (x>=cv_x(i)-cv_int) & (x<cv_x(i)+cv_int);
    
31. %     for j=1:n_cv
    
32. %         ind_y= (y>=cv_y(j)-cv_int) & (y<cv_y(j)+cv_int);
    
33. %         for k=1:n_cv
    
34. %             ind_z= (z>=cv_z(k)-cv_int) & (z<cv_z(k)+cv_int);
    
35. %             p_n = sum(ind_x & ind_y & ind_z);
    
36. %             cv_V(j,i,k)= p_n;
    
37. %         end
    
38. %     end
    
39. % end
    
40. %
    
41. % save Mat_cv.mat
    
42. 
    
43. load Mat_cv.mat
    
44. cv_V=smooth3(cv_V,'box',3);
    
45. 
    
46. % potential( u can mod the value with your style  ---------------------
    
47. cv_V=-log(cv_V./sum(cv_V(:)));
    
48. tm = cv_V(cv_V~=Inf);
    
49. cv_V(cv_V==Inf)=max(tm);
    
50. cv_V=cv_V-min(cv_V(:));
    
51. cv_V(cv_V>3.7)=NaN;
    
52. 
    
53. % slice u want 2 plot ---------------------------------------
    
54. xslice =  cv_x(1:1:n_cv);
    
55. yslice =  cv_y(1:1:n_cv);
    
56. zslice =  cv_z(1:1:n_cv);
    
57. 
    
58. h=slice(cv_X,cv_Y,cv_Z,cv_V,xslice,yslice,zslice);
    
59. for i = 1:3*n_cv
    
60. set(h(i),'EdgeColor','none','FaceColor','interp',...
    
61.    'FaceAlpha',0.01+0.008*(n_cv-mod(i,n_cv))) ;
    
62. end
    
63. 
    
64. % set only one alpha value ----------------------------------
    
65. % set(h,'EdgeColor','none','FaceColor','interp',...
    
66. %    'FaceAlpha',0.2) ;
    
67. 
    
68. % choose the colormap u want ---------------------------------
    
69. colormap(flipud(hot))
    
70. % colormap((hot))
    
71. 
    
72. % just figure setting -----------------------------------------
    
73. axis tight equal
    
74. set(gca,'GridLineStyle',':')
    
75. box on
    
76. ax=gca;
    
77. ax.BoxStyle = 'full';
    
78. grid on
    
79. axis([-pi pi -pi pi -pi pi])
    
80. view(20,20)
    
81. 
    
82. h = gcf;
    
83. fig_na=['./file_3d_grid_hot'];
    
84. fun_work_li_035_myfig_out(h,fig_na,3)
    
85. 
    
86. %% logs
    
87. % mod : 2017年 09月 20日 星期三 09:15:39 CST
    
88. % if u have question about this code. contact me
    
89. % meatball1982@163.com
    
90. %
    

复制代码

---

欢迎光临 Mathematica中文论坛-非官方 (http://ilovemathematica.com/)

Powered by Discuz! X3.2